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Definition 

Entrainment velocity = v = — \ 

p dx 

Fluctuating component of vertical velocity 

Cartesian coordinate in streamwise direction; 
origin at leading edge of plate 

Downstream distance between leading edge of 
plate and turbulence manipulator 

Cartesian coordinate taken normal to flat 
plate ; origin at plate surface- 
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Symbol 


Definition 
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Wall shear stress 
Kinematic viscosity 

Non-dimensional damping coefficien 


ABSTRACT 


The potential of passive "manipulators" for altering 
the large-scale turbulent structures in boundary layers 
was investigated. Utilizing sr.oke-wire visualization and 
multi-sensor probes, the experiment verified that the 
outer scales, defined by the intermittent excursions of 
potential fluid into the boundary layer could be 
suppressed by simple arrangements of parallel-plates, A 
parametric investigation of different types of 
parallel-plate devices verified that optimum designs 
existed for the different configurations of stacked and 
tandem plate arrangements. A final optimized design 
consisting of two single plates residing at a height of 
0.8- in the boundary layer and spaced 8c apart in the 
streamwise direction was chosen for the final set af 
measurements to document the characteristics of the 
"manipulated" layers. 

As a result of suppressing the outer scales in the 
turbulent layers, a decrease in the streamwise growth of 
the boundary layer thickness was achieved and was coupled 
with a 30 percent decrease in the local wall-friction 
coefficient. After accounting for the drag on the 
manipulator plates, the net drag reduction reached a 
value of 20 percent within 55 boundary layer thicknesses 
downstream of the device. No evidence fcr the 
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1 reoccurrence of the outer scales was present at this 

streamwise distance thereby suggesting that -further 
reductions in the net drag are attainable. The 
manipulated layers also exhibited a 17 percent thickening 
of the wall sublayer. In turn, the frequency of the wall 
"burst" process was reduced by a comparable amount. 

Utilizing new techniques, many of them introduced 
in chis investigation, the objective detection and 
description of the two-dimensional velocity field near 
the wall was achieved. These included match filtering, 
two-dimensional velocity reconstruction as well as 
digital image processing and pattern recognition. The 
results of these established that the frequency of 
occurrence of the wall events is simultaneously dependent 
on the two parameters, 2 nd The correlation of 

the inner layer Strouhal frequency for these events with 
the two parameters resulted in a simple relation which 
correlates wall-layer event (or "burst") results of this 
investigation as well as those in the literature, 
covering a range of Reynolds numbers from 600 to 38000. 

As a result of being able to independently control the 
• , inner and outer boundary layer characteristics with these 

manipulators, a new view of these layers emerged. It now 

I 

I appears that the large outer scales in laboratory 

\ turbulent boundary layers are remnants of the transition 

f process, i.e., slowly decaying turbulent spots that are 

i 

B 

i xxiii 

m 

■ 

i 

! 


ATXX 


^ 'aouaxnqjn:} aqq jo ,But6b, 

laqmnu spiou^an sq^ uaqseq ajojajaqq pue aauaxnqinq 
JO oqj qJT*' mutt syqq aAomaj qoagja 

uj saoqcxndTueui aqx 'JaAex Aacpunoq aqq UT pappaqoa 



CHAPTER I 
INTRODUCTION 


The use of screens, grids and honeycombs for 
-managing- turbulence has been an ongoing interest at 
I.I.T. for some time. The application of such devices 
for controlling free stream turbulence in wind tunnels is 
best summarized in the publications by Loehrke and Nagib 
(1972 and 1976 ) and Tan-atichat et al. ( 1982 ). Through 
careful scaling of mesh sizes and in-depth appreciation 
of the various suppression and generation mechanisms that 
control the turbulent flow, these techniques have been 
quite successful in this application. 

Recently, a growing energy consciousness has 
generatea interest in developing techniques for reducing 
viscous drag on aerodynamic bodies. Recent uork 
conducted by NASA's Langley Research Center by D. 

Bushnell and J. Hefner and their colleagues has focused 
on this subject. Their work is best summarised in the 

. reference by Bushnell ( 1978 ). In a majority of 

applications, this would involve removing or altering, in 
some way, the turbulent boundary layers that es^st on 
their surfaces. It was the latter of these two 
alternatives that was of interest to us and which lends 
itself as a natural estension to turbulence management 
techniques. To be most effective, however, it is 
essential that these -turbulence manipulators- be 
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tailored to the particular mechanisms governing the 
flowfield of interest. Devices such as screens and 
honeycombs are not necessarily the suitable ones here. 

This warranted a review of the current experimental and 
theoretical information on wall-bounded- turbulent shear 
flows . 

Relevant Literature 

In’ the last 20 years, the research cn turbulent 
boundary layers has transformed the picture given by 
earlier long-time average statistics. Once perceived as 
quasi-steady turbulent eddies transported by a mean 
shear, such flows now appear, as a result of space-time 
velocity correlations and short-time conditioned 
m.easu r ements , to be highly unsteady coherent motions of 
strongly interacting scales. This revolution in thinking 
was prim.arily set forth by Townsend (1956) and Grant 
(195B) who, while being primarily interested in turbulent 
wake flows, viewed the turbulent flowfield as one where 
turbulent fluid is moved about by the motions of a system 
of large coherent eddies whose dimensions scaled with the 
width of the flow. The motions of these large scales 
carried eddies of smaller comparative size but which 
contained most of the turbulent energy. Townsend in 
particular inferred from velocity correlation 
measurements, the presence of characteristic eddies of 
finite length with a structure inclined in the flow 
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direction. He sumised that these large eddies in 
turbulent shear flows form coherent and identifiable 
groups and control the overall rate of spreading by 
contorting the layer between the turbulent and the 
non-turbulent fluid and by affecting the bulk convection 
of turbulence energy from regions of maximum production. 
Later Townsend (1970) described these structures as an 
inclined pair of double— roller eddies, a picture which was 

consistent with the eddy pattern perceived by Grant 
(1958), 

At about the same time, Kline and others at Stanford 
University were studying turbulent boundary layers with 
the then new hydrogen bubble visualization technique. 

The culmination of that work was the 1967 paper by Kline 
et al. which revealed the presence of surprisingly 
well-organized motions in the near-wall region of a 
turbulent boundary layer. These motions lead to the 
formation of low-speed wall "streaks" which eventually 
interacted with the outer portions of the flow through a 
process "of gradual lift-up, sudden oscillation, bursting 
and ejection". This sequence of events later became 
known as "bursting". It was felt that this process 
played an important and dominant role in the production 
of new turbulence and the transport of turbulence away 
fro- the wall. 

The turbulent boundary layer is uniquely different 
from any of the other turbulent flcwfields in that its 





mean velocity profile can be divided into four main 
•parts. The first of these is the wall region in which 
the flew processes are dominated by the presence of the 
wall. Within the wall region in the near vicinity of the 
wall is the viscous sublayer (0 j< y+ £ 7) . In this 
region^ viscosity plays a dominant role, however 
turbulent fluctuations are still large relative to the 
local mean velocity. Scaling analysis in that region 
yields the result that u+ = y+ in the sublayer and 
subsequently is also used to define it. Further from the 
wall, 50 j< y+ £ 200 , there is a region of nearly 
homogeneous turbulent flow which is often referred to as 
the constant stress layer. In this region the velocity 
is found to vary with height in a logarithmic manner 
defined by u"^ = ^ In y^ + B, The viscous sublayer 
and constant stress layer are smoothly matched by the 
buffer layer which can take a number of implicit forms to 
define the velocity distribution (Whitfield, 1979; 
Spalding, 1981) . The buffer layer typically exhibits the 
greatest amount of turbulence production. Measurements by 
Klebanoff (1954) showed that th^ peak production occurred 
in that region just at the edge of the sublayer. Beyond 
the constant stress layer, the . 'i^.lence is bound jd by 
the potential flew outside the dary layer. Coles 

(1956) designated this as the wake region because of its 
similarity to flew characteristics exhibited in a wake 
flew. He further developed an empirical formulation of 


the mean velocity distribution in that region which was 
found to work well for a variety of boundary layer flows 
with pressure gradients. At the outer edge of the wake 
region^ there exists a thin interfacial region separating 
the turbulent and the non-turbulent flow. This interface 
was termed the superlayer by Corrsin and Kistler (1955). 

It is across this region that vorticity is acquired by 
the potential flow and eventually is incorporated into 
the wake flow. 

Each of these regions has a characteristic length 
scale. In the viscous sublayer that length is v/u^. The 
wake region has a characteristic length on the order of 
the boundary layer thickness. Kovas?.nay (lObT) estimates 
the thickness of the superlayer to bo 10\/V^^ where is 
the entrainment velocity. 

The fact that the maximum turbulence production 
coincided with the region of large coherent activityr 
discovered by Kline ct al., prompted more investigations 
of this near wal 1 phenomena . For examplCr Corino and 
Brcdkey (1969) visualized the near wall motions in a 
turbulent pipe flew. In their paper they reported that 
the most important and distinguishing characteristic of 
the wall region was the intermittent ejection of discrete 
fluid elements away from the wall. They catagorized the 
process as a local deceleration of the fluid near the 
wall followed by the appearance of faster moving fluid 
which approached from the outer parts upstream resulting 
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in eruptions of low velocity fluid away from the wall. 

The ejection of low velocity fluid was always found to be 
followed by a sweep of upstream fluid from the previous 
ejection event. They also investigated the Reynolds 
number dependence of these events and found that their 
frequency of occurrence increased like the 1.7 power of 
Re . Corino and Brodkey (1969) further estimated that 

X 

the wall ejections accounted for 70 percent of the uv 
Reynolds stress in the boundary layer. This amount was 
substantiated by the visual study made by Kim et al. 

(1971) . The "burst" process examined by that reference 
may be somewhat altered in boundary layer flows over 
large roughness elements as pointed out by Grass (1971) . 
Using a uv quadrant discrimination VJillmarth and Lu 
(1972) determined that 60 percent of the uv Reynolds 
stress is produced when the sublayer velocity is less 
than the mean and that the largest contribution occurs 
when u<0 and u>0. 

A growing issue at that timie was how the frequency 
of occurrence of these events scaled with the boundary 
layer characteristic lengths and velocities. In the 
initial studies by Kline et al.(1967) the "bursting" was 
associated with the wall streaks and, therefore, thought 
to scale with the inner variables and v. The results 
of Rao et al. (1971), which utilized narrow band-pass 
filtering of a differentiated velocity tim.e series to 
detect the intermittent turbulent fluctuations associated 
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with these events, indicated that their frequency scaled 
with the outer variables, and 6. Kim et al. (1971) 

tried to reconcile the difference in their proposed 
mechanism by suggesting that the low-speed streak lifting 
was triggered by large* scale disturbances. There was, 
however, no explanation for the clear dependence of the 
burst period on the wall friction velocity in that 
reference. Through some dimensional arguments and 

r • 

assumptions about the boundary layer profile shape, 

Laufer and Narayanan (1971) suggested that the 
nondimensional burst frequency Uo^T/o = 5, was 
approxim.ately 5, which is in near agreement with the 
value of 6 obtained by Rao et al.(1971). They further 
pointed out that the average wavelength of the 
turbulent-laminar interface is approximately 
suggesting an interaction between this outer flovr 

structure and the inner layer regular events. • 

In a recent review article by Cantwell (1981) he 
states that "it appears fairly well established that the 
mean time between bursts scales with outer variables". 

However, he prefaces his statement by saying that the 
interaction process that creates conditions under which 
the bursts occur is far from understood. Blackwelder and 
Haritonidis (1981), however, in their recent paper in 
which they utilize a locally high variance as a burst 
discriminator (see also Blackwelder and Kaplan, 1976) , 
claim that the burscing period scales with inner 
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variables which is in support of the Kline et al. (1967) 
concept based on the lift up of the low-speed wall 
streaks. A more detailed investigation of this aspect 
was performed by Blackwelder and Eckelmann (1979), where 
they identified the low-speed streaks as the regions 

between pairs of counter-rotating streamwise vortices. 

« 

This verified an earlier result obtained by Bakewell and 
Lumley (1967) who performed an eigenfunction 
decomposition of space-time correlations of streamwise 
velocity fluctuations in the wall region and obtained a 
dominant eigenfunction structure consisting of 
counter-rotating eddy pairs of elongated stream\’ise 
extent . 

In order to form a conceptualized view which would 
lead us to distinguish between inner and outer variable 
dependence of the wall "bursting" process, we must review 
the available data on the organized motions in the outer 
part of the turbulent boundary layer. Since the 
measurements of Corrsin and Kistler (1955) of 
intermittency in the outer parts of the boundary layer, 
it became clear that the turbulent/non-turbulent 
interface extends well into the boundary layer. The 
stuoies of intermittency by Fiedler and Head (1966) gave 
visual evidence of this. This region is of interest 
because it is the interface for energy transfer from the 
external flow into the turbulence producing region. 

Kaplan and Laufer (1969) utilized a rake of hot-wires 
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spanning the whole boundary layer and came to the 
conclusion that the non-turbulent intermittent regions 
did not appear to be surrounded by turbulent fluid but 
were connected to the free-stream fluid. The interface 
also seemed to be highly corrugated. Kovasznay (1970 ) r 
utilizing conditional sampling technigues, defined 
-fronts" and "backs" to these turbulent "bulges". He 
found that the "front" velocity was approximately 5 
percent faster than that for the "back". The average 
velocity of the turbulent bulges was of the order of 0.93 
which suggested that the outer flow rides over these 
Etructi^res. From velocity correlations Blackwelder and . 
Kovasznay (1972) later defined a mean motion inside a 
turbulent "bulge" which was rotational with a mean 
vorticity component in the same direction as the mean 
strain rate. Furthermore, they gave these features a 
travel life-time, estimated from a 1/e correlation 
decay, to be from 5 to 10 boundary layer thicknesses. 
From their measurements they estimated that these large 
eddies contribute as much as 80 percent of the Reynolds 
stress for y>0.2d. The conditioned measurements of 
Antonia (1972) were in substantial agreement with them. 
In addition he determined that the maximum uv 
distribution within the turbulent bulges represents 
approximately 45 percent of the wall shear value. This 
provided support to the Bradshaw et al. (196/) theory 
that the strength of the large scale eddies is closely 


related to the wall shear stress, t.^. 

On the trailing edge or "back" of the large scale 
eddies, Falco (1977) noticed predooinant small-scale 
eddies, with sizes on the order of 100 to 2C0 v/u_, which 
he associated with Reynolds-stress-oroducing motions. 

Head and Bandyopadhyay (1951) present a picture of the 
larce eddy structures as beinc an ar.alca-ation of 
elongated hairpin vortices. It was suggested by them 
that Falco' s (1977) small scale eddies are actually 
tips of these hairpin structures. 

In spite of the well established existence of these 
large scale motions and their obvious importance to the 
entrainment of potential fluid, turbulence energy 
transfer and wall friction, the mechanisms for their 
generation are still in question. Zilbermann et al. 

(1977) and v:ygnanski (1975) generated turbulent spots 
(Emmons, 1951) in a transitional boundary layer flow and 
followed their development downstream. They concluded 
that the turbulent spots were related to the large scale 
motion in a turbulent boundary layer. Schraub and Kline 
(1965) suggested that these bulges are the same as or at 
least formed from the ejected wall streaks. This was 
similarly conjectured in the results of Kovasznay et 
al.(1970). Later Offen and Kline (1975) suggested that 
the large-scale outer bulges night be a consequence of 
vortex pairing between the vertices associated with two 
or more wall bursts. They prefaced this concept by 
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stating that it was related to data taken at fairly low 
Reynolds numbers. However, they suggested that the basic 
structure should be the same at high and low Reynolds 
numbers, Taneda (1981) proposes that the large scale 
bulges might be a preferred eigenmode of the outer flow. 

The suggestion by Offen and Kline (1975) tnat the 
large outer scales are a result of the pairing of old 
bursts is only part of their hypothesis. They further 
say that the circulatory flow or vortex in a "burst" 
impose an adverse pressure gradient at the wall as it i^ 
convected downstream, thereby resulting in t..e — ft~up 
of a newly formed wall streak. This defined a 
regenerative process for the wall bursts. Coles (1978) 
discussed and evaluated a large volume of boundary layer 
and ducted flow data based on a Gdrtler~like instability 
mechanism. This mechanism would lead to the streaks and 
may be triggered by the flow curvature associated with 
the large scales. 

Theoretical models for a regenerative bursting 
process based on the instability of the sublayer were 
proposed by Hanratty (1956), Einstein and Li (1957) and 
Black (1966) . According to theil: analysis the bursting 
frequency should scale with the inner variables v/u^. 

This view cannot, however, be easily reconciled with tbe 
findings of Rao et al. (1971) which indicate that the 
bursting period shows Reynolds number independence when 
scaled with outer variables. This suggests tr*e influe. -ce 
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of the outer large-scale eddy motions. The dependence of 
the semi-recular wall events on the passage of an outer 
large-scale eddy was also supported by Piaturi and 
Brodkey (1978) , based on their interpretation of 
visualization results, Doligalski et al. (198C) proposed 
a model for the lift up of low-speed fluid from the wall 
which was based on the interaction with an eddy structure 
which originated in the outer part of the flow. 

Blackwelder and Woo (1974) investigated the bursting 
process while applying a pressure perturbation on thf flow 
of a type exhibited by a passing large scale eddy. Those 
results indicated no changes in the wall bursts as a 
result of the pressure excitation, Chen and 
Blackwelder's (1978) temperature measurements above a 
heated wall highlight an internal shear layer by which 
the large scale outer structure may be related to the 
wall bursting. Brown and Thomas (1977) suggest that an 
organized structure, inclined in the flow direction at an 
angle of 18 degrees and is on the order of 2 i long, 
relates the outer flow to the wall bursts. This 
structure, they believe, may be connected with the 
interaction of turbulent spots. 

If we summarize these results into one picture for 
the process governing the turbulent boundary layer, we 
can point to the existence of two types of coherent 
motions, one associated with the wall layer and the other 
with the outer flow. Each of these has been observed to 
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contribute a large percentage of the Reynolds stress in 
the boundary layer. The wall layer is associated with 
the generation of low-speed streaks which lift away from 
the wall, interact with faster moving fluid and 
eventually break-up. This process has been associated 
with the generation of as much as 70 percent of the net 
turbulence production. Away from the w’all, the outer 
turbulent bulges maintain a net circulation which plays 
an important role in the entrainment of potential fluid 
and the extraction of turbulence energy from the mean 
flow. The strength of these large scale eddies has been 
found to be related to the wall shear stress, although 
the mechanism, for their generation is not clear. Through 
a conceptualized view it has been thought that these two 
coherent motions should interact, however, the method by 
which this occurs is again not clear. Of primary 
importance to this is the manner in which these events 
scale. From an intuitive standpoint, it would seem that 
the period of these wall layer events or "bursts" should 
scale with wall units. The data in this regard seem to 
be contradictory, although it is safe to say that no 
clear trend has been put forth. There is also little 
agreement in the literature as to the source of the 
large-scale outer eddies. It is undeniable, however, 
that these features of the flow are important in 
sustaining the generation of new turbulence and in the 
growth of the boundary layer. 
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Objectives 

In light of the fundamental questions that remain 
unanswered about the coherent processes which have bean 
observed to play a major role in the growth and evolution 
of a turbulent boundary layer, we undertook an experiment 
which would separate the two classes of coherent features 
and point to the mechanisms of their interrelations. 
Loehrke and Nagib (1976) and Tan-atichat et al. (1982) 
had a great deal of success with this approach. This 
involved altering the outer part of the turbulent 
boundary layer by suppressing the large-scale 
intermittent turbulent bulges. Based on some of the 
conjectures previously cited, the absence of this feature 
in the outer part of the flow would either interrupt or 
impede the "bursting" process. On the »ther hand, 
another view may suggest that these wall layer events 
would continue undisturbed and their outer interaction 
would generate new large scale bulges, or some new 
equilibrium would be established and remain throughout 
the lifetime of the boundary layer. Also of extreme 
importance from a technological point of view is the 
prospect of reducing skin— friction drag beneath turbulent 
boundary layers by removing the intermittent excursion of 
higher velocity potential fluid near the wall. We 
proposed to use two-dimensional "boundary layer 
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manipulators" formed from thin, stacked parallel plates 
which would be suspended in the boundary layer parallel 
to the wall. This. concept to reduce drag was first 
proposed by H. M. Nagib in 1977 (see Bushnell, 1978) . 

These plates would present a minimum obstruction to the 
flow and provide us with the greatest amount of design 
freedom for matching the device to the scales of 
interest. Other devices for altering the large scale 
turbulence in a boundary layer are summarized by Hefner 
et al. (1980). A parametric study was first planned in 
order to pinpoint the important criteria for controlling 
the large scale motions. 

To aid as tools for understanding this complex 
process we initiated our research intent on utilizing 
advanced flow visualization and velocity measurement 
techniques. These included digital processing techniques 
which could provide consistent detection and statistical 
analysis of these coherent events. The techniques would 
include one- and two-dimensional analysis of single 
velocity sensors and rakes of sensors, along with the 
simultaneous acquisition of digitized images of the 
smoke-wire visualized flowfiela- Signal processing 
techniques such as one-dimensional match filtering and 
two-dimensional reconstruction^ and match filtering would 
be employed to isolate coherent features in the velocity 
time series and digital images. 

The results of this study are aimed at documenting a 
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method by which net drag reduction could be obtained in 
fully turbulent boundary layer flows, Ke further 
intended to document the coherent features in the 
turbulent boundary layer and to try to answer sone of the 
outstanding questions about their interaction. 

Many thanks are due to John L. Way for his ready 
advice on the signal processing aspects as well as the 
expert work he has done in building the data proceSiSing 
system. The expert craftsmanship of Edward D. Nieman 
exhibited in the modification of the wind tunnel and in 
the other fabrications needed throughout the experiment 
is also acknowledged. The careful work of Jeffrey A. 

Corke in helpina to fabricate the analog preprocessing 
circuitry is appreciated. Special thanks to Amy K. Lee 
and Dennis J. Koga for helping to prepare the final figures. 
Finally, special appreciation is due to Roberta J. Corke and 
Thelma T/ Grymes for their help in typina the jneinuscript . 

It is obvious that such as this is built on and 

shaped by the efforts and accomplishments of others who 
have preceded us. To them we extend a special thanks- The 
authors would like co especially acknowledge the efforts 
and stimulating discussions with Dennis Bushnell and Jerry 
Hefner, the technical monitor of the grant, of the NASA 
Langley Research Center. 
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CHAPTER II 


EXPERIMENTAL FACILITIES AND TECHNIQUES 
Facilities 

The experiment was performed in the high-speed test 
section of the I.I.T. Environmental Wind Tunnel. The 
tunnel runs in a closed return configuration with the 
return leg being used as a low-speed test-section for 
atmospheric boundary layer simulations. In the 
high-speed section, honeycombs and screens upstream of a 
4:1 contraction give a turbulence intensity level of 0.03 
percent over a range of velocities from 1 to 30 m/s. The 
rectangular test section is 0.61 m wide and 0.91 m high 
with plexiglass back and front walls to allow for easy 
flow visualization. In addition, the front walls are 
hinged in two sections to allow easy access to the test 
section. 

Wind Tunnel Modifications . During the experiment, 
between the preliminary and final data runs, 
modifications were made in the high-speed leg of the wind 
tunnel to accommodate the addition of 1.83 m. to the 
streamwise dimension of the test— section. To accomplish 
this, a 2.44 m. long section from the 4.88 n.. long 
diffuser, downstream of the test section, was removed. 

The existing 3.05 m. long instrumented test section was 
moved back toward the remaining diffuser and a 1.83 m. 
long parallel-wall section made the connection to the 
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contraction. A re:?.ovablc 0.5 by 1.3 in. plexiglass 
window was incorporated into the new section and all of 
the other walls were painted black to aid in flow 
visualization. 

A 0.58 m. long square diffuser section was made to 
connect the test section to the remaining diffuser 
section. This "quick" diffuser was bolted to the test 
section. A 2.54 cm. cap was left between the two 
diffuser sections to allow air to enter the tunnel and 
maintain its internal pressure above the room atmospheric 
pressure. A schematic of the wind tunnel and the new 
4.88 m. long high-speed test section is shown in Figures 
1 and 2 . 

Without the boundary layer plat? inside the test ^ 
section, measurements were m.ade to determine if these 
modifications .lad any detrimental erfects on the tunnel 
performance, using a digital, real-time spectrum, 
analyzer, special attention was given to the free-strean 
turbulence level and especially its low frequency 
content, where we felt the effects of any diffuser 
separation would appear. These measurements indicated 
that the tunnel modifications had no adverse effect on 
its performance and actually appear to have improved the 
tunnel efficiency as indicated by a lower operating air 
temperature rise and aerodynamic noise levels. 

The reason for these improvements may be in the 
(jj^ffuser design. Before the modification, air entered 
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the tunnel through a 0.6 cm. slot in the wall perimeter 
just downstream of the test section, it now appears that 
the air was jetting sufficiently enough to have caused 
the flow in the long diffuser to separate. Since this 
flow leads directly into the fan assembly, 
non-uniformities in the flow could result in inefficient 
fan operation. In the present case, the combination of 
the "quick" diffuser and the four times wider slot 
results in lower velocity air to enter in a direction 
more tangent to the "slow" diffuser wall. 

Traversing Mechanism . The downstream 3 m. length of 
the tunnel test-section is equipped with a motorized 
traversing mechanism. Using this, two probes can be 
moved to any position in a plane 0.6 m. high and 3 m. 
long respectively along the test-section vertical 
ccnterplanc. The two-dimensional probe positions are 
monitored through the linear outputs of two 10 turn poten- 
tiometers acting as voltage dividers. The potentiometers 
are mounted on the movable beds carrying the probes and 
their translational motion is transTerred through rack 
and pinion gears. This ensures that the outputs are 
proportional only to the relative locations of the probes 
and not sensitive to the tolerances in the traversing 
drive mechanism. In the vertical direction, a probe can 
be repeatably positioned with 0.03 mm. u^ccuracy. 
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Boundary Layer Plate . A 5.16 in. long, 0.57 m. wide 
flat plate served as the boundary layer plate for the 
experiment. To be able to place the plate inside the 
test-section, it was broken up into three parts with 
lengths 1.68, 1.85 and 1,63 m. respectively, from the 
leading edge. This arrangement is shown in the schematic 
of the plate in Figure 3. 

The first plate was fabricated from a 1.6 cm. thick 
aluminum plate. The leading edge of the plate is 
machineo to a razor edge. In addition, the noise is 
sloped downward so as to fix the stagnation point on the 
measurement side of the plate. The other two plates are 
made from 1.3 cm. thick clear plexiglass which were 
mounted to a metal frame made from 2.5 cm. square 
aluminum channel to maintain a flat surface. Each plate 
section was suspended on four independently adjustable 
legs. The vertical placement of the plates partitioned 
the test section so that 0.61 ir.. by 0.61 m. souare 
cross-section was formed above the plate. • The seams 
between the three plates were glued together with epoxy 
and drawn tight by threaded rods attached to each end of 
adjoining plates. To guarantee a smooth surface across 
the seams, they were covered with a plastic filling 
compound and wet-sanded to provide a smooth, flush 
surface. 

V^ith the boundary layer plate inside the test 
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section, the flow at its leading edge was observed using 
smoke-wire visualization. As a result of the greater 
blockage below the plate due to the support legs, the 
stagnation point was observed to occur on the underside 
of the plate which resulted in flow separation on the top 
measurement side. In order to equalize the blockage, a 
number of perforated plates were tested at the exit of 
the test-section across the part of the cross-section 
above the plate. A 40 percent solidity plate was found 
to give the best results. With this plate, a slightly 
higher pressure was maintained on the measurement side of 
the test-section so that the stagnation streamline 
impinged on the downward sloping leading edge. This 
slight pressure difference was further utilized by 
-leaving a .16 cm. gap along the intersection of the 
boundary layer plate and side walls. This resulted in a 
very slight flow from this intersection to the underside 
of the plate. As observed by injecting smoke into the 
corner, the influence of this downflow only extended out 
to approximately 1 cm. from the side walls. Since the 
. corners are likely sites for the generation of 
longitudinal vorticity from the merging of the plate and 
side-wall boundary layers, we wanted to draw off as much 
of the boundary layer there as possible, without 
interrupting the two-dimensionality of the test boundary 
layer . 
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Surface Plug . A 6.03 cm. circular hole was machined 
through one of tht plexiglass plates to accomragdate a 
removable plexiglass flush-surface plug. The 
instrumented plug shown in the photograph of Figure 4 was 
made up of a larger plug which fit snugly the hole in the 
piste and a smaller f 1.27 cm. diameter removable plug 
placed off-center, C.16 cm. from the edge of the larger 
^olug. An 0-ring filled a groove cut into the 
circumference of each plug to ensure an air-tight fit 
when they were in place. 

The smaller plug was designed to hold the surface 
wire. In order to minimize heat conductivity through its 
surface, it was constructed of wood. Two gold plated 
terminals were imbedded in its surface, 1.3 cm. apart, 
and connected to copper leads which extended out through 
the underside of the plug. The surface wire, designed to 
measure wall shear velocity, consisted of a 10 micron 
diameter tungsten wire which was copper plated on both 
ends and soldered to the terminals. Care was taken here 
to ensure that the wire was taught and that it rested on 
the surface. The unplated portion of the wire was 
approximately 1 mm, long. This results in a maximum 

imens i ona 1 sensor length, 1 = iu_/v , for the range of 

shear velocities encountered approximately 20. This 
value is approximately half of the width of the 
'*3ll“Streak structures (Blackwelder and Eckelmann, 1979), 

The larger plug was designed to hold a rake pf 10 
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hot-wires so that their sensors were aligned in a 
vertical plane on. top of the shear wire. The rake was 
made by soldering the ends of jeweler's broaches to a 
double-sided, copper clad p.c. board. Separate 
conduction lines were etched into both sides of the board 
by exposing the photo-sensitized surfaces to ultraviolet 
light through a photo-negative mask of the network. The 
** broaches have a minimum length of 3.1 cm. and their 

diameters tapered from. 0,3 mm. to 0,08 nun, at the sensor. 

The sensor spacing was selected for specific 
non-dimensional heights in the low Reynolds number 
boundary layer case at that downstream station. Because 
of the high mean shear near the floor, the wires are more 
closely spaced there. The non-dimensional heights of the 
sensors vary slightly from the design values in the other 
boundary layer, because of the change in the wall shear 
velocity. However, there was always a minimum cf six 
wires below y+ of 100 in all cases. 

The exact alignment, of the broaches for soldering 
onto the p.c. board v.-as accomplished on a specially made 
jig. The jig held the sensor ends of the broaches at the 
precise spacing and angles so that the opposite ends 
could oe soldered to the p.c. board. After mounting the 
broaches, thin glass fibers were glued across broach 
pairs to reduce vibration which could result in 
strain-gaging of the velocity sensors. Ten micron, 
copper-plated tungsten wire was soldered onto the tips of 
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the broaches. As with the shear sensor, the unplateo 
portion of the wires nave a maximum sensor length, 1 ,of 
approximately 20 , 

Tha P-C, board carryinq the probes extended up 
through the plug in a 2.9 cm. long, 0.15 cm. wide slot 
cut in the side diagonally opposite to the shear-wire 
plug. With the sensors on the top of the plug, the 
bottom of the board was fixed to a micrometer mounted to 
the underside of the plug. When the plug was placed in 
the boundary layer plate in the test section, the 
micrometer was used to lower the wire rake into position 
near the floor. In all cases, the lowest wire was 
positioned to be at y+ = 7. A twenty two conductor 
ribbon cable and connecter carried the power to the 
eleven hot wires ir the plug. As with the boundary layer 
plate, the- seams around the plug were filled with plastic 
filler and wet-sanded smooth. 

Boundary Layer Kanipulator . In the course of this 
study, a number of boundary layer manipulators were 
studied. All of these were configured from thin, flat 
plates which \?ere suspended above the boundary layer 
plate and spanned the tunnel test section. The plates 
were used singly, vertically stacked or placed in tandem 
in the streamwise direction. Because of the number of 
param.eters involved in the configuration of these 
m^anipulators (namely, the chord length and thickness of 
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the plates, the vertical spacing between then, their 
angle of attack and the distance between tanden plates as 
well as the relationship of all of these to the boundary 
layer .characteristic length scales) a nethod was needed 
to -optinite- their design. This was acconplished by 
holding the plates of various dimensions in a^ low 
disturbance rig inside the tunnel test section. The tig 
consisted of two sets of vertical, slotted rods which 
were supported on either side of the test section in = 
tail on the test section ceiling and in the snsll gap 
between the boundary layer plate and the side wall. Each 
manipulator plate was. clamped in a streamlined holder 
which could slide in the slotted rods. The plates were 
positioned and tightened in place so that they were 
tensioned to remove any sagging or ripples. Various 
sites of manipulator plates were placed in the rig for 
different spacings in three different turbulerit boundary 
layers. Using a smoke-wire to visualize the flow, we 
adjusted the parameters for four different basic plate 
configurations to obtain a maximum visible reduction in 
the large-scale, intermittent structures. 

The final design for 'the manipulator used in the 
quantitative measurements was optimized in the manne* 
previously described. The placement along the boundary 
laver plate was chosen so that the minimum Reynolds 
number Re,., was sufficiently large (20001 so that the 
layer is fully turbulent past the transition stage and 



containing a full spectrum of scales. The streamwise 
chosen location also provided a sufficient number of 
boundary layer thicknesses (approximately 60) from the 
start of the manipulator along which we could observe any 
regeneration of the large scale structures. Also, since 
measurement station 1 of Figure 3 was the farthest 
upstream location that could be conveniently measured 
from the traversing mechanism, it was necessary that the 
manipulator be placed at a point where the approaching 
flow field was well docujnented. 

The final manipulator consisted of two, 0.2 mm. 
thick brass plates. The chord of each of the plates and 
its elevation above the floor was equal to 80 percent of 
the boundary layer thickness at the location of the 
upstream most plate, i.e., equal to 3.2 cm. In the high 
Reynolds number case the boundary layer thickness at that 
station did not change sufficiently enough to change any- 
of the design values. The two plates were spaced 25.4 
cm. apart or eight times the boundary layer thickness at 
the upstream plate. The manipulator plates were held at 
one end by a 3.2 cm. wide, 0.7 mm, thick steel bracket 
which wrapped around the end of the plexiglass boundary 
layer plate, and was secured to thd aluminum frame below. 
On that end the brass manipulator plates were bent over 
the brackets and fixed by three small machine screws. The 
opposite end of each manipulator plate was attached to a 
0.7 mm. thick steel bracket which pivoted on a small 
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threaded rod. The rod protruded through the back wall of 
the test section and a nut on the outside was used to 
draw up the slack on the plates and put them under 
tension. Special care was taken to make sure that the 
manipulator plates were not sagging, twisted or under an 
angle of attack. In addition they were observed under a 
strobe light to ensure that they were not vibrating 
under the operating conditions. A photograph of the two 
tandem manipulators mounted in the wind tunnel for the 
final set of measurements is shown in th"e top portion of 
Figure 4. 

Instrumentation 

A variety of analog and digital instrumentation was 
used in the investigation. In general, the analog 
instrumentation was used to pre-process the data before 
the digital acquisition by a PDP-11/10 mini-computer 
con*^ oiling an A/D converter and other peripherals. 

Pr ssing of the data was performed later on a variety 
of computer systems including a larger UNIVAC 11/08. The 
analog instrumentation and peripherals used in the data 
acquisition and processing will be described in the 
following . 

Multi-channel Anemometers . To control the rake of 
11 hot wires, a system of 12 constant temperature 
anemometers based on a variation of the design first 
introduced by Kreider (1973) was constructed. The 


circuit diagram for the anemometer is shown ir. Figure 5. ' 
The system represents a low cost alternative to 
commercially available CTA's while giving comparable 
performance to a DISA 55DC1 anemometer (Wle2ien,1981) . 

The anemometers were grouped into two sets of six 
channels each. Two 15 volt,il Ampere power supplies, 
whose circuit diagram is shown in Figure 6, provided 
ample power for the anemometer operation. The circuits 
utilized CA3140 hi-gh-precision operational amplifiers 
which are pin for pin replacements for more traditional 
741 op-amps, but which require only a few milliamps to 
power and provide a near ideal input impedance of 10^ 
ohms and frequency response bandwidth of iMHz. To guard 
against cross-talk between channels, all of the power 
inputs were AC coupled with 0.1 uF capacitors to ground. 
Each set of six anemcmeters also utilized a built' in 
square wave generator which was used to adjust the 
feedback in the servo loop. 

The operation of the anemometers involved measuring 
the resistance of the not— wire and adjusting the overheat 
resistance to equal times the overheat ratio times 20, 
where the 20 accounts for the 20:1 bridge ratio. With 
the switch set to the +15 volt side, the 10 ohm 
resistance is adjusted to provide an output voltage of 
approximately 10 volts. Kith the servo loop closed and 
the square wave supplied to the bridge, the 10 ohm 
resistance is adjusted to give a slightly underdamped 
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response. The square wave generator was turned off 
during the normal operation of the anemometer channels. 


Multi-channel Signal Conditioners . ^ Prior to 
digitization of the analog outputs, the signals were c.c. 
biased, amplified and lew-pass filtered by the circuit 
shown- in Ficure 7. The bias and amplification was 
necessary to minimize the quantization error in the 

analog to digital conversion. The low-pass filtering vas 

« 

necessary to remove the amplitude energy in the signals 
above the ::vquist frequency for spectral analysis. 

The bias and gain are accomplished in the first tvo 
stages of the op-amp circuit. The gain in each stage is 
set by plug-in feed-back resistors. In this experiment, 
the first stage operated as a unity gain amplifier so as 
not to saturate the op-amp, and a majority of the d.c. 
was subtracted from the signal. The second stage 
amplification was normally 10 or 15 and any remaining 

d.c. was subsequently subtracted from the signal. The 

. • 

exact amplification was determined by applying a 
precisely known d.c., voltage into the, channel and 
measuring the output value. The last two stages of the 
circuit are two second-order low-pass butterworth 
filters. In series they give the fourth order response 
shown in Figure 8 with an amplitude decay rate of 
60dE/decade for frequencies above the filter cutoff. The 
damping values for these filters were fixed by the 
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capacitors in the circuit to give the desired amplitude 
response with frequency . The cutoff frequency was* 
changed by switching in dual pairs of precision resistors 
to give nominal cutoff yalues of 2.5, 5.0, 8.0 and 10.0 K 
hertz. The resistance values corresponding to these 
frequency cutcff values are given in the table of Figure 
7 . Provisions were made in the design to allow for 
2 ®thing of other cutoff values through the equation 

f^ = 1/27 tRC^C 

Digital Acquisition and Processing System 

The digital system (D.A.P.S.) for acquiring analog 
data and performing digital processing consists of three 
systems: The Fast Acquisition System (F.A.S.), the 

Processing and Development System (P.D.S.) , and the 
Special Purpose Module (S.P.M.) . Each of these systems, 
schematically shown in Figures 9, 10, and 11 was used in 
the experiment to develop software, to perform the 
digital acquisition, to process the data or in many 
cases, to display or plot results. For some of the 
results, proce*Ssing was also performed on a UNIVAC IICS. 
The FAS and PDS systems are based around PDP-11/10 and 
PDP-11/04 mini-computers, respectively. Information 
could be passed between the three systems through 
external PDP-11 bus lines, but typically the SPM was 
connected to either the FAS or PDS at one time. All 
these systems are portable, thereby allowing them to be 
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to the laboratory during data acquisition. These 
systems will be described in the following with more 
detailed attention being given to those devices which are 
an integral part of this experiment. 

Fast _^quisition Sys^. The FAS, schematically 
shown in Figure 9, was used along with the Special 
Purpose Module to acquire the data during the experiment. 
It consists of a PDP-ll/lO processor with 28 K, 16 bit 
words of RAM memory. The analog signals are acquired by 
a micro-processor controlled, 16 channel, 100 kHz 
analog-to-digital converter which is designated as the 
MIP. This device features its own 4K words of dual 
ported memory. The memory is double buffered so that one 
of the 2K buffers can be written to tape by the PDP-11 
processor while the MIP microprocessor fills the other 
buffer, thereby allowing continuous, uninterrupted 
acquisition of data. The high speed digital magnetic 

tap, d.ive in ,he FAS allows fast, continuous storage of 
•3 mega-words of data., 

Spec^ Purpose Module . The S?M, schem'atically 
shown in Figure 10, is made up of a number of special 
purpose peripherals which are used both in data 
acquisition and processing. Data from and control of 
these peripherals are handled by one of the PDP-11 
processors through the bus. Digital acquisition of ' 
analog signals is done in this module by a 100 kHz, 16 
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channel, 12 bit, A/D converter. Unlike the 
micro-computer controlled MI?, acquisition here is 
controlled by the PDP-11 processor. Acquisition rates 
are set with the programmable clock in the SPH, The 
acquisition software is written in computationally 
efficient PDP-11 assembly language which handles the A/D 
control, clock interrupts and'flata conversion at a rate 
which makes full use of the 100 kHz A/D. As in the case 
of the. other assembly language software" Written- for this, 
experiment, the macro acquisition software was configured 
as subroutines which were called from less efficient but 
more I/O compatible Fortran programs. Computer control 
of analog devices is made through the 3 channels of 12 
bit D/A conversion in the SPM. Again the D/A converter 
is controlled through Fortran callable, assembly language 
subroutines. 

The other three peripherals in the irK, i.e., the 
digital camera, color display processor and array 
oroccssor, performed very important functions in this 
exoeriment. Because of their unique nature and 
imocrhtance, they will be given a more detailed discussion 
in the following. 

Di gital Camera . A model 659F optical data digitizer 
which' is interfaced to the PDp^ll bus Has used to 
digitize the smoke-wire visualized flowfields. The 
camera incorporates a Vidicon imaging tube which stores 



an electron charge pattern proportional to the luminous 

energy of an image focused on the tube. As contrasted 

with a scanning type video tube in which the image data 

is retrieved by sequential sweeps, the vidicon tube has 

4096 X 4096 randomly addressable picture elements 

(pixels) from which intensity information can be 

retrieved. As a result of some spatial averaging of 

adjacent pixels, tfee number of resolvable points for the 

tube is 512 x 512. The analog intensity information- is 

converted- to an 8 bit parallel binary intensity . word 

through an A/D converter in the camera. This gives 256 

resolvable levels of intensity. Since the PDP-11 word is 

two 8 bit bytes, the image data was packed into double 

bytes to reduce storage requirements. The vidicon tube 

-6 

is sensitive to light energy in the range of 10 to 10 
footcandle-seconds which gives it an equivalent ASA 
rating of approximately 3000. The exposure time is 
program controllable from 15 ms to 60 s. , 

The camera also has an analog video and position D/A 
output and a 30 frame/scc sweep generator so that the 
image can be viewed on a black and white monitor. This 
feature is useful for setting up the camera before 
digitizing runs. The camera is equipped with a Nikon 
lens mount so that a conplemcnt of conventional, lenses 
can be. used. The camera data and control lines are 
interfaced to the PDP-11 bus through a DR-llC camera 
interface and a merging logic unit. Camera control and 
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data transfer were handled by a Fortran callable, 
assembly language subroutine. The image data is stored 
on digital magnetic tape for processing later. 

Display Processor . An MDP/3 display processor is 

used to display the image data retrieved from the digital 

camera. The MDP is a micro-processor controlled image 

processor having 64K x 18 'bits Words of refresh memory 

with 16 bits for graphics display and 2 bits for 2 

* 

graphic overlays in each processor word, as well as iK by 
16 bit RAM memory and 4K by 1’6 bit EPROM memory. Image 
data ‘is stored in the graphics memory, assigned a gray 
level intensity value or a pseudo color through 

look— up tables and displayed on a color 
monitor. The monitor screen formats are 512 x 512 pixels 
X 4 bits to give 16 levels of intensity or colors and 256 
x 256 X 8 bits to give 255 levels of intensity or colors. 
In addition, special firmware allows rapid manipulation 
of displayed data such as level windowing, zoom and 
linear interpolation to 128 x 1^8 or 64 x 64 as well as 
animation. 

A library of Fortran callable, assembly language 
subroutines was written in order to load up the KDP 
processor registers and to execute the desired firmware 
programs. In addition, programs were written in MDP 
micro assembly language and loaded into the MDP memory to 
allow data to be written directly into and retrieved from 
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the MDP graphics memory. Thus we were able to utilize 
the 65K words of MDP memory as a temporary storage media 
or "scratch pad". This was especially important because 

the PDP-*ll/04 and 10 can support only 28K words of 

* ' • • 

iT\emory. 

Array Processor , An MSP-^3X array processor was asec 
to perforin the floating point arithmetic on the digital 
image and hot-wire signal arrays. The MSP is a 
block-floating dedicated array processor with 4K words of 
24 bit data memory which is controlled by a high speed 
microprocessor. As with the MDP, 4K x 16 bit EPROMS 
contain firmware to support the execution of a 
comprehensive array function library. These functions 
includ*^ real and complex number addition, multiplication, 
magnitudes and forward and i-nverse fast Fourier 
transforms. As an indication of the processing speed of 
the MSP, it can perform a 1024 complex FFT in 26 
milliseconds. This is extremely important since it 
requires 512 256-point one-dimensional row and column 
FFT's to transform a typical 256 x256 image into the 
spatial frequency domain for filtering and pattern 
recognition. As with the MDP, a library of Fortran 
callable, PDP-11 assembly language subroutines were 
written to transfer data into the MSP memory and to 
execute the array arithmetic. 
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Processing and Development System . The PDS system 
shown in Figure 11 was used to develop programs and 
process image and hot-wire data. This system is 

controlled by a PDP-11/04 mini-computer with 28 K words 

« 

cff“16 bit memory. Efficient use of the processor, hard 
and floppy disks, magnetic tape drive and input/output 
devices is made through the DEC RT-11 operating system. 
Under this system programs can be created, edited, 
compiled, linked with subroutine libraries and executed. 

Flow Visualization 

A "smoke-wire" was used to visualize the boundary 
layer flow at the lower freestream speed of approximately 
6 m/s in this investigation. The smoke-wire, described 
by Corke, et al, (1977), consists of %0.1 mm. diameter 
stainless steel wire which is coated with oil by letting 
a drop run down the wire. As a result of a capillary 
instability, the oil forms droplets along the wire so 
that when the wire is resistively heated the vaporized 
droplets generate discrete streaklines of smoke. The 
heating of the "smoke-wire" and triggering of a camera 
and strobe light are controlled by a synchronization 
circuit. Digital logic in the circuit allows 
visualization records to be taken which are conditioned 
on an "event" in the flow. 

In this experiment a smoke wire, oriented normal to 
the boundary layer plate, was placed midway between 
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stations 6 and 7. The upper support arm was placed well 
into the free-stream so as not to disturb any part of the 
flow in the boundary layer. The stainless steel 
"smoke-wire" was connected to the upper support arm and 
passed through a- small hole drilled through the boundary 

layer plate to its underside. That end of the wire was 

* 

clamped and earth grounded through a connection to the 
metal* frame of the wind tunnel. When the smoke wire was 
not in the tunnel, the small hole in the boundary layer 
plate was filled and wet-sanded smooth so that it would 
not c’reate any disturbance to the flow. 

The streamwise position of the smoke-wire allowed us 
to visualize approximately the last 15 boundary layer 
thicknesses of the flew. Also, the smoke wire was 
sufficiently upstream of the plug station so as not to 
influence the hot-wire rake by its thermal wake. 
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CHAPTER III 
DATA ACQUISITION 


The digital data acquisition of the analog signals 
for this experiment can be broken up into t*he relatively 
straightforward process used in the calibration of the 
hot-wires and time-mean measurements in the boundary 
layer, and into the much more complicated acquisition for 
the instantaneous unsteady measurements and digital image 
retrieval. All of the data acquisition software was 
written for this experim.ent in a manner to make the most 
use of the computer core available and to optimize the 
use of the various peripherals described in the previous 
chapter. These experimental setups and methods will be 
described in the following sections of this chapter. 


Hot-wi re Calibration 

The -calibration of the sir.g-le hot-wire used in the 
mean measurement runs and the rake of 11 wires used in 
the unsteady measuremer.t was done in the wind tunnel 
test-sectipn. A schematic of the hot-wire rake is shown 
in the corner of Figure 3. Sven though the tunnel 
heating problem was improved by the test section 
modifications, the test section air temperature would 
still increase by approximately 4^C for the low Reynolds 
number runs and approximately 6-8-C during the high 
Reynolds number runs. The absolute tunnel temperature 
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was r however, dependent on the laboratory room temperature 
which could vary, about 3 to 4°C from day to day. For 
this reason, all calibration and data runs required that 
the hot-wire outputs be compensated for low frequency 
temperature variations. 

f 

For the single-senSor measurements, the hot-wire, 
sensitive to the streamwise velocity, U, was operated in 
a constant temperature mode by a hISA 55D01 anemometer. 

In this case and in the multiple wire case, a 1.8 
overheat was used to provide the minimum amount of 
ambient tem.perature sensitivity of the output voltage 
proportional to velocity. Two separate calibration runs 
were made for the low and high Reynolds number cases so 
as to calibrate over each range of velocities with the 
two different mean operating temperatures. The overheat 
temperature was set after the tunnel temperature had 
reached its equilibrium value. Since the absolute 
temperature varied from day to day, the calibration 
temiperature as referenced to an average long term 
temperature from, which the compensated, true velocity 
could be determined. 

During the calibration run for a single hot-wire, 
three channels of analog signals were digitized. One was 
the analog output of a Digatec 5810 temperature 
transducer proportional to the wind-tunnel air 
temperature. This signal was pre-processed through the 
bias/filter/gain circuitry to obtain an output equal to 
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the temperature in degrees Celsius divided hy 10 . A 
n,inimum, pre-designed cutoff frequency of 2.5 KHz. was 
used to low-pass filter the signal prior to digitization. 
With the minimum resolution of 4 rov. for .-he A/D. 
converter, an accuracy of ,1°C was obtained. . 

The absolute test-section velocity was measured by a 
pitot-static probe connected to a Validyne DP103 pressure 
transducer. This transducer has a sensitivity of 2.6o 
vAg/m in pressure, thereby giving it a sensitivity of 

2.52 V /m/s in velocity. This signal was directly 
digitized to give a resolution of approximately .002 m/s 
over the range of operating velocities. 

The output of the anemometer was d.c. biased, 
amplified nominally by 15 timss and lo«-pass filterad at 
5 KHZ prior to bain, input to the A/D converter. The 
exact bias voltage and amplification factors were 
determined vith a test signal to account for tolerance 
limits of the op-amp -circuitry. These three sensor 
probes were placed at the same streamwise location along 
the mid-span of the tunnel test section for calibration. 
The three signals were simultaneously sampled and simple 
statistics were computed on line during the calibration 

run. 

With the tunnel 'pre-heated- to one of the operating 
temperatures, the fan, speed was controlled to give a 
range of steady velocity readings. At each of these 
setungs, acquisition was initiated and a sampling of 
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points was taken until the voltage values converged to 
steady means. This typically, involved 1024 points taken 
over 4 seconds. Along with the mean values, the maximum 
and minimum values in the sample were ouptut to give an 
indication of the sample distribution. Each of the 
’digitized voltages, proportional to temperature, was 
converted to temperature through the calibration and the 
average temperature was computed. Given the average 
temperaturef the digitized voltages proportional to the 
pitot-static pressure w^re converted to velocity and the 
average free-stream velocity was computed. 

From each of the voltages proportional to the 
streamwise velocity, the. gain and bias were removed to 
give the actual bridge voltage. The temperature 
corrected voltage was then computed from the relation 
given by Bearman (1971) . 


E(T) 

corr 


[ 1 



- 1/2 

(T-T^)^ E(T) 


where 

r = overheat ratio at thd reference temperature 
o 

= resistance sensitivity to temperature of 

tungsten (^C) ^ 

T = actual temperature 

T = reference ter.oerature 
o • ’ ' • * 


The value 


for 


a was 
c 


experimentally determined by 
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allowing the tunnel to heat up at a fixed velocity. A 
typical value was approximately .Of)34°C. The average 
corrected bridge voltage was computed for each velocity 
setting to generate voltage-velocity data pairs. 

Given the data pairs, a squared, third order 

. 

polynomial least-squares curve fit was performed to find 
the "best-fit" coefficients. This fit had been 
demonstrated by Drubka and Klezien (1979) to bettef . 
represent the output from a hot-wire than King's Law and 
to be computationally more efficient. A typical average 
error in the fit was 0.0C6 m/s. 

In the case of the rake of hct wires, the 10 wires 
counted to the p.c. board were held in the free-stream 
for calibration. The hot-wires were operated by 10 of 
the 12 CTA channels built for the experiment. The 10 
analog outputs were d.c. biased, amplified nominally tfy 
10 times and lew-pass filtered at 5 KHz. These 10 
channels were simultaneously sampled with the signals 
proportional to temperature a.'.d pressure as with the 
single wire calibration. The bridge voltage outputs were 
temperature compensated and 10 sets of velocity-voltage 
data pairs were fit to determine their calibration 
polynoninal-coeff icier, ts . 

Shear-wire Calibration . A single hot-wire mounted 
to the tunnel floor was calibrated to measure the 
streamwise component of the surface shear stress. The 



wire, which was operated in a constant temperature mode, 
was calibrated in a turbulent boundary layer at different 
Reynolds numbers. An estimate of the mean wall shear was 
inferred from a single wire placed in the si^layer above 
the shear wire, which measured the slope of the velocity 
distribution near the wall. Although there is ev^^ence 
from measurements in the thick sublayers of turbulent 
boundary layers in oil (Ekelmann,1974) that the sublayer 
velocity distribution is not linear, for sufficiently 
small non-dimensipnal heights, the error in assuming a 
' linear distribution is small. For this reason the 
calibration was performed in the manipulated boundary 
layer which, as a result of its lower wall-f r icticn 
velocity, has a thicker sublayer in which a probe can be 
placed at a smaller non-dimensional distance from the 
wall. The height of the probe was carefully measured 
using a cathatom.eter . The results of the ^calibration are 
shewn in Figure 12. This calibration was considered 
satisfactory for the present study, since the main 
purpose for using the shear wire was to obtain 
information on the fluctuating shear stresses at the 
wall. Hi^h reliability for absolute values of the shear 
stress was not a primary requirement for this method.’ .• 

The analysis of heat transfer from a heated strip it^ 
a linear velocity distribution which was rer.ovec my 
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(1972> , suggests a 1/3 power relation ^_wee.. 


the shear stress and heat less, or 
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Ag + B 


where B is related to the energy lost into the surface 
substrate on which the sensor is mounted. Per constant 
temperature operation, the i^R heating equals E^/R\ 

W , ^ ^ 

Which equals CE^ since R is constant. Therefore, 


T = CE^ + B 

CJ 


was used to normalize the axis of Figure 12. Sandborn 

(19B1) points out through a oerturbaticn analvsis that - 

'' ' w 

is a function not only of E but also of the nor.ents of 
fluctuations up to e^ . Therefore, the mean value of the 
anemometer output voltage is composed of a mean value 
from conditions when no fluctuations are present plus 
higher order mean terms due to fluctuations in the shear 
as in the case of the turbulent boundary layer. These 
conclusions point to the fact that the calibration must 
be done where the conditions, such as the turbulence 
levels and scales, are similar to those in which 
measurements are to be taken, 

# 

ft is obvious in Figure 12 that the calibration 
would not be best fit by a straight line, and as in the 
case of the velocity sensing wires, a squared third-order 
polynomial was found to provide the best fit tc the 
voltage-shear distribution. The average error in the 
calibration fit was 0.3 percent. 


44 


Acquisition for Mean Measurements 

The digital acquisition of data for tii»e-mea„ 
results was broken up into profile runs which utilired 
relatively slow -acquisition rates for the calculation of 
thf first and second nonents of the strear-ise velocity 
fluctuations, and into relatively high acquisition rates 
for the calculation of one-dinensional power spectral 
density, auto correlations and integral ti.e scales. 

Er^le Riips. The instrumentation schematic 
fer the acquisition of the boundary layer profile data is 
Shown in Figure 13. As in the case for the single-wire 
calibration, the voltages proportional to the tunnel 
temperature, free-stream velocity and instantaneous 
velocity from the hot-wire were digitally sampled. m 
addition to these, the voltage proportional to the 
vertical position of the hot-wire was also simultaneously 
sampled. This voltage was pre-processed to give an 
output of lv/12.5 mm. sc as to give a resolution of 0.05 
through the A/D conversion. 

Since the profile runs were intended to calculate 
only the mean and rms velocity values, it was necessary 
only to have enough samples .taken over a sufficiently 
long period of time, equivalent to the passage of a 
sofficient number of the 'integral scales of interest,, for 
the quantities to reach steady values. As contrasted 
vith the spectra, the mcnents are not aliased in any wav 



by the sampling rate- For these data runs, a sequence of 
four channels,, sampled 64 times a second for 30 seconds 
was found to be sufficient. 

Profile data was taken by starting at • e fi-st 
upstream station and digitally acquiri/:q ar . storing on 
digital magnetic tape the 5012 values comprising a tape 
record for each of 30 vertical heights from the surface 
of the boundary layer plate. Nearer to the plate, in 
that region of greater mean shear., the vertical sampling 
intervals were more closely spaced. In all but the f-’rst 
and last profiles, the stream.wise intervals between the 

m 

data stations was 0r3 shewn in Figjre 

To ensure that the voltage corresponding to the 
vertical position of the probe accurately indicated the 
relative distance frci?, the plate surface at each station^ 
the lowest point was always set to be at x=0.5 mm, and to 
correspond to zero volts output. This was dene by 
sighting the probe and its reflection in the plate in a 
graticuled cathetometer , By this techniquOf the probe 
censor coulc? be positioned to an accuracy of ,01 rnin. 

Spectra Runs . The data runs for spectra 
calculations were cone for half of the streair.wise 

* 

stations. In order to be able to make cemparisons 
between the regular and manipulated boundary layers at 
the different stations, the data were taken at the same 
self similar heights of y+ = 8, 15, 3.0, 50, 75, 100 , 200, 
400 and 80C. 
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The h-’ghest -frequency of interest in the velocity 
fluctuations was 5 KHz so that a sampling frequency of 10 
KHz was necessary to ensure that the spectra would not be 
aliased. To further ensure no aliasing the signal was 
low-pass fil'tered at 5 KHz. For each height at each 
station, 100 records of 2052 pts/record were acquired and 
stored on digital magnetic tape. Each record contained 
0.21 seconds of continuous data so that scales as large 
as 0 .7 m. {a minimum of 9 boundary layer thicknesses) 
being convected at 70 percent of the free-strear 
velocity could be resolved, in regards to the turbulence, 
the 100 records were statistically independent, however, 
tV.f- acquisition program was written so that the data be 
acquired in phase with the 60 Hz clock, in this way any 
60 cycle line noise and all of its harmonics could be 
removed if necessary as demonstrated by Way (1975) . 

Sim.ultaneous Hot-Wire R ake and Digital Image Acquisition 
The simultaneous acquisition of the rake of 10 
wires, the shear wire and the digitized image of the 
smoke-wire visualized flowfield was designed to look in 
detail at the relationship be^tween the large-scale outer 
structures and the wall "burst* mechanism in the 
turbulent boundary layers. The design of the probe 
concentrated as many as 7 sensors below y''’=100 which has 
been observed to bound the influence of the wall events. 
Although other sensors as high as y =800 would provide 
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some information about the large, scale turbulehce* the 
digital images, viewing the whole boundary layer over a 
number of boundary layer thicknesses, would be primarily 
used for that information. 

The method by which the data was acquired reflects 
the type of information we wanted to extract and the 
limitations of the equipment. For example, we wanted the 
sampling rate to be fast enough to give sufficient 
resolution of points during a "burst" event and to be 
continuous for a sufficient amount of time to ensure more 
than one event per sampling record. For a number of 

i s t i ca 1 ly independent samples, the second criterion 
was necessary to compute frequency of events. The first 
criterion was limited only by the speed of the A/D 
converter and controlling assembly language software. 

For 13 input channels a maximum sequence sampling rate of 
approximately 2000 Hz could be attained. The second 
criterion was limited by sampling rate and the amount of 
computer memories available. Since the 65 K bytes of 
image data could not be retrieved from the digital camera 
fast enough to keep pace with the hot-wire acquisition, 
only one image per sampling was possible. Therefore, the 
images had to be "conditioned" on the event. Since we 
wanted the images and hot-wire data to coincide, the 
detection ot the event had to also initiate the hot-wire 
acquisition. However, we were also interested in the 
pre-event history of the flow. The acquisition strategy 
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was designed to allow for all these requirements. 

The results of Rao et.al. (1971), Laufer and 

Narayanan (1971) , Blackwelder and Kaplan (1976) and 

others have suggested that in the range of Reynolds 

numbers of this investigation the non-dimensional period 

between "burst" events is TUc/6 =5. In the case of the 

low Reynolds number regular boundary layer, a sampling 

- » 

frequency of approximately 1800 Hz would resolve this 
period into 100 points. Approximately 400 words of 
memory were available to store each of the 11 channels of 
anemometer outputs^ temperature and pressure so that at 
this sampling rate an average of 4 events could be 
detected. The acquisition rate was adjusted to maintain 
this resolution over the range of boundary layer 
parameters. Table 1 lists the acquisition rates for all 
of the boundary layer cases. 

The analog and digital instrumentation setup for the 
simultaneous acquisition is shown in Figure 14. The 11 
channels of anemometer outputs are d.c, biased, amplified 
nominally by 10 and low-pass filtered at 5 KHz before 
being input to channels 1 - 11 of the A/D converter in * 
the support processing module.* The analog signal 
proportional to the. free-stream pitot-static pressure 
difference was input to channel 13 of the A/D. 

The analog signal output of channel 10 of the 
anemometer was used to detect the "burst" events. This 
signal, which was proportional to the instantaneous 
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Table 1. Sampling Rate Parameters for Simultaneous Acqui- 
sition . • - * • 



B. L. Case 

(m/s) 

C (cm) 

T(s) 

f (Hz) 

SSR 

PTS/ T 

Regular 

6,71 

• 

8.08 

.0590 

16.96 

1694 

lOQ 

Regular 

10.36 

7 . 4 2 

.0356 

28.10 

2000* - 

T 

/ _ 

Manipulated 

6.71 

7.16. 

.0536 

18.66 

1865 

IOC 

Manipulated 

10.36 

6.71 

.0323 

30.92 

2000* 

65 

* Upper 

acquisition 

, rate 

limit 

for 13 channels 



velocity at y-r = 7 in 

the 

boundary layer, would be 


examined for occurrences 

of 

high variance compared 

to the 

long-time background 

values 

. Although it has been 

common " 

for investigators to 

use 

a 

sensor at y-f = 15, the r 

esults 


of Blackwelder and Kaplan (1976) show that the velocity 
signature of these events remains coherent in the whole 
portion of the boundary layer below y+ =100. The signal 

4 

at y* = 15f being near the peak in the vertical 
distribution of turbulence intensity, has considerably 
more turbulence fluctuations thereby making the detection 
more difficult. The (Blackv:elder and Kaplan, 1976) 
technique in effect low-pass filters the velocity data to 
perform the detection^ However in this case, analog 
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filtering would be undesirable since it would introduce a 
phase shift in the detection signal. The signal from the 
y+ = 7 sensor allowed easy detection without low— pass 
filtering. In all of the boundary layer cases, with the 
changing free-stream velocities and the addition of the 
manipulator, the vertical position of the hot-wire rake 
was adjusted to compensate for the different values of 
the wall shear velocity so as to always locate the lowest 
wire (channel 510) at y+ = 7. 

The anemometer output of the detection wire was 
linearized with a DISA 55D10 lineariz'er and input to a 
DISA 55D35 rms meter to measure the long-time average rms 
and to output a signal proportional to the instantaneous 
variance. The analog sigr>al was fed to the variable 
threshold level detection circuit schsmatically shown in 
Figure 15a. The output from the I detection circuit 

consisted of a TTL logic status cc ^ onding. to 


CONDITION 


1 0: u'Iq < ku-2 

^ 2 2 

(1: U',, < ku-2 


where k is an adjustable threshold coefficient. The 
logic output Was connected to the condition logic in the 
smoke-wire control. 

The smoke-wire circuit is made up of timLtg chips 
which control the duration of heating of the wire and the 
delay needed for the smoke to be ccnvected into the field 
of view before triggering a camera and/or strcbe light. 
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In addition, digital logic in the circuit al!^ows 
photographs to be taken when the smoke is in view and 
when an external condition is satisfied. The conditioned 
photographs are necessary to do ensemble image averaging. 
The timing diagrcun for the smoke-wire control circuit is 
shown in part in Figure 16. The sequence is initiated by 
a positive pulse which can be manually applied or, as in 
this case, was software controlled through the output of 
channel 0 of the D/A converter. 

The experiment was run in a darkened lab so that the 
image was recorded at the instant of the strobe flash 
The flash is approximately 1 ms. in duration which was 
sufficiently short to "freeze" the motion of the smoke 
image. In order to be able to mark that instant on the 
sequence of data, the trigger output pulse was 
multiplexed onto the slowly varying signal proportional 
to the wind tunnel air temperature. Before this could be 
cone however, it was necessary to "shape" the pulse using, 
the circuit scherr.atically shown in Figure 15b. This 
circuit was designed to transform the syncing pulse frcm 
the strobe light into a -iCv amplitude, 10 ms. duration 
square wave pulse. The 10 ms. guaranteed a minimum of 20 
detection points at the A/D conversion rates used in the 
acquisition. The shaped pulse was added to the 
temperature signal through the second stage of the op-amp 
amplifier section of the bias/filter/gain circuit. The 
combined signal was fed to channel 12 of the A/B 
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converter . 


The. digital camera, which was mounted on a tripod, 
was connected t J the PDP-11 bus through the rerging Ire: 
unit and DRll-c interface. Under software ccr.trol, tie 
camera was operated in the 30 Hz scanning mode while tie 
D/A outputs of the camera were viewed on a Tettronix t!4 
black/white monitor which sat on top of the camera. 
During this operation, the camera performed liice a 
conventional television camera w’hich cade it easy to 
position and focus. For both of the low Reynolds numter 
boundary layer cases the camera was set up to record tvo 
views of the flow. One view (zoomed cut view? which 
encompassed approximately 8 boundary layer thicknesses 
utilized a 50 mm., f;1.4 lens. The other (zocmed in 
view) encompassed approximately 2 boundary thicknesses 
and utilized a 210 mm., f:3.5 macro/zeem lens. In both 
cases, the rake of hot vires was visible in the image. 
The camera was not moved throughout this phase of the 
acquisition. To be able to compare the i.T=ge results 
with future data runs, the ca.rera always viewer two LET 
lamps which were precisely spaced to provide a scale. 
During the retrieval and storage of the image fata, the 
digital numbers were copied into the KDP-3 disrlay ' 

processor graphics memory and displayed or. the Conrac ' 
5211C color display monitor. The MDP-3 graphics memory 
was also used for temporary storage of the digitized 13 
channels of data. 


The acquisition was directed by the PDP-11/04 
processor in the fast acquisition system as controlled by 
the acquisition software. The timing diagrars for the 
simultaneous acquisition is shown in Figure 16.. After 
the smoke-wire is coated with oil, the sequence is. 
initiated by a keyboard command. Prior to this, the 
digital camera is in a prime mode which erases any 
previous images. After initiation, the camera expose 
lines are enabled aw^ting a flash from the strobe to 
store an image. A start pulse is sent to the smoke — wire 
control circuit to begin the generation of the smoke 
streak-lines. The duration of the smoke and the delay 
time to allow for the smoke to be convected had been 
preset in the smoke-wire control, it is within this 
"window" of time that the smoke is in the field of view 
and during which we search for a condition. The "window" 
was typically 1-2 seconds wide. The condition is only 
valid if it occurs when the smoke is in view. 

since we do not know a priori where in tCiis vincow 
an eient w'ill occur , the acquisition of the 13 cbanmels 
began at the initiation of the sequence and the data were 
stored ift computer memory. For the sampling rates used, 
the amount of available data storage was too small to 
accommodate the time interval within our conditicnal 
window. Therefore, the software configured the space 
in.to a "circular" buffer so that old data were 

f 

overwritten. Between each sequence of 13 channels, and 
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before the next interrupt of the programmable clock, the 
data value of channel 12 of the A/D was checked to see if 
the condition, marked by a -lOV pulse, had occurred. If 
the condition was detected, 50 points, or one half of an 
average "burst" cycle, preceding the event were kept and 
the remaining 350 points following the event were 
acquired to fill the buffer. Also at the instant the 
condition occurred, the strobe flashed and the smoke 
image was stored on the digital camera. Tkg 13 channels 
of data were then copied into a file on digital magnetic 
tape. The digital image was retrieved from the camera 
and copied to the next file on the digital magnetic tape. 

The camera was subsequently put in the prime mode 

awaiting the next .sequence. 

Depending on the delay time allowed for the 

convection o'f the smoke, it was possible that 50 points 

were not acquired before an event was detected. In such 
a case, the data set was discarded. It was also possible 
that no events be detected inside the "window". In this 
case, the acquisition loop was halted by manually 
applyi''g a -10 v signal from the pulse shaping circuit 
into the channel 12 A/D input 'to terminate the 
acquisition cycle. Again the data set was discarded. 

The threshold for the detection was set using an HP 
digital Spectrum Analyzer. The output of the threshold 
detector circuit was connected to the external trigger 
input of the spectrum analyzer to conditionally sample 
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the events. The threshold was raised. until the educted 

ensemble average resesbled the signature for the "burst" 

process. The ratio of the instantaneous threshold 

variance level to the long-time average variance used in 

the regular and roanipclated boundary layers was 3.5. 

This number is not to be confused with the threshold used 

in the VITA technique since those- values are the local 
2 

average over 10Tu_/v points. 

For the high Reynolds number boundary layer cases , 
the velocity was too high to operate the smoke wire. In 
those cases, where the digital camera was not used, the 
acquisition was net initiated by the detection of a 
"burst" event. Rather, 256 statistically independent 
records of 400 points per channel for the 13 digitized 
analog signals were sequentially acquired and stored on 
digital magnetic tape. 

Acquisition runs were also made for use in cemputing 
the power spectrum, auto-correlation and probability 
density of the signals from the 10 rake sensors and shear 
wire. In this case, the output of each wire was 
independently acquired at a 10 KHz sampling rate and 
stored on digital magnetic tape. The methods and 
acquisition software were identical to those described 
earlier in this chapter. 
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CHAPTER IV 


DIGITAL PROCESSING OF TRANSDUCER OUTPUTS . 

This chapter is intended to describe the various 
types of digital processing performed on the hot-wire 
data. These include the methods for the boundary laver 
profile calculations, and for determining integral 
boundary layer properties, as well as the approach used 
in calculating the momentum balance for obtaining the 
wall friction of the boundary layer. Special 
one-dimensional time sequence processing is next 
examined. This includes one-dimensional power spectrum, 
auto-correlation, probability density and integral length 
scale calculations. Also included is one-dimensional 
match filtering for pattern recognition of the signature 
associated with the regular wall phenomena. The latter 
will provide information about the frequency of their 
occurrence, i.e., their distribution in time. The 
two-dimensional processing will demonstrate 
two-dimensional sampling and reconstruction techniques 
which allowed comparisons to be made between the hot-wire 
rake data and the digital images. Techniques for 
processing of the images based on smoke-wire 
visualization are described in the next chapter. 

Boundary Layer Profile Calculations 

The boundary a.ayer profile data were made up of 20 
discrete data points, spaced in the vertical directicn 
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normal to the flat plate, taken at 10 (or 9 in the 
Manipulated boundary layers! streamwise stations.* The 
data were stored in raw form on four channels of digital 
magnetic tape as time series of unlinearized anemometer 
output and values proportional to the vertical probe 
position, tunnel air temperature and free-stream velocity 
for each of the discrete profile points. The first step 
in forming the profile calculations was to put the 
data in a more manageable form. For each profile point, 
this involved computing the mean air temperature, 
free-stream velocity and vertical probe position through 
their respective calibration equations. The anemometer 
output was compensated for any drift in tunnel 
temperature from the calibration temperature, in the 
manner described in Chapter III, and the mean and rms 
velocities were computed. These five values were written 
on another digital magnetic tape where they would be used 
for further processing. Since this processing, as well 
as the other processing on the profile data and raw data 
calibration for the hot-wire rake, required many floating 
point operations and tape-to-tape data transfers, they 
were performed on a large main frame UNIVAC 11/08 
computer which handled these more efficiently than our 
FDP-11 systems. Kith the average profile values stored 
on digital tape, the integral characteristics of the 
boundary layers were calculated. Since we were 
particularly interested in the friction coefficient, a 
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number -of these integral quantities were calculated in 


the same program which computed C^. 

A number of empirical relations exist for 

determining the wall shear from velocity profile data in 

a turbulent boundary layer. One of the most commonly 

used involves fittinc the "log law" relation 

+ , + 
u = A log y + B 

* m 

Here A and B are- ^xperir.entally determinad^.and 

« 

Oniversally accepted constants. A number cf values for 
these coefficients have been reported, and there is some 
evidence that they may he weakly dependent on pressure 
gradient (Hinze, 1975). The values of 2. *539 for A and 
5,0 for B have been recommended by Coles (1969) and are 
most commonly used. 

In a zero pressure gradient turbulent boundary layer 
flov;, Coles’ "Law of the v:ake" (1960, 1969) relations 


— = i loa{ 1) ^ B - ^ 

k V ' k 


and 


5.U 
1 « 
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can be combined to remove the effect cf boundary layer 
thickness to obtain 




c.U - 65 k 
1 

\‘ 

i + n 


+ B + 


2T. 

k 


59 


OF POOR QUALITY 

Given 6;^ and U. from the velocity profile, and using the 
value for n given by Coles (1960), the wall shear 
velocity can be obtained. The value for k and B are 
those given by the log relation. The value of 7. has been 
observed (Hinze, 1976) to be strongly dependent on the 

streamwise pressure gradient. 

\The empirical relation of Ludwig and Tillmann (1945) 


-.268 


-.678 H 


= .246 Pe^ 


10 


12 


has also been used to infer the friction coefficient from 
the integral profile properties. This relation was 
derived by correlating the experimental results from 
boundary layer measurements of different investigators. 

Although these relations could be used to infer the 
friction velocity from the profiles of the boundary 
layers without the manipulator, there is no reason to 
expect that the empirical constants used in any of these 
formulations would take the same values in the 
m..nipulated boundary layer. Rather than using any of 
these, we undertook to measure the balance of momentum 
from which we could obtain the shear force on the plate. 
This was done for both the manipulated and 
n n-manipulated boundary layers with the empirical 
relations being used onl^ for comparison purposes. 

If we solve the momentum equation for a 
two-dimensional, incompressible turbulent boundary layer 
we can obtain the integral form 
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which was first formulated by yon Karroan in 1921 (see 
Hinze, 1975) . The contribution of turbulent fluctuations 
to the balance are trivial and neglected in this 
formulation. The terms 5^ and are the integral 
displacement and momentum thicknesses 


6 = \a - -U(y))dy 

and 

i“iX> ,1 - 

^ CC CD 

The computational method for solving this equation with 
the prescribed profile data is given in the flow-chart of 
Figure 17 and described in the following. 

For each particular boundary layer case, the mean 
velocity and turbulence intensity values corresponding to 
the vertical y positions at each streamwise station were 
read from digital magnetic tape into a computer array. 

For use in comparing values to those dete..inined by the 
momentum balance, the U,y data pajrs were least-squares 
curve fit to the log relation over the range from y+ = 50 
to 200 to determine fot the profiles. This was 

done in an iterative manner by initially guessing a value 
for u and computing’ the average deviation of the points 
from the equation 
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u - 2-44 log y"*" + 5.0 

A correction to the previous was computed and the 
solution continued to loop until a "best fit" was 
obtained. 

The free-stream velocity, , at each station was 

computed from the top portion of each profile. Since the 

value for Ua, appears inside the integral for 6 and 6 

1 2 

small errors in measuring that quantity are summed up in 
the area under the distribution. We therefore formulated 
an accurate general method for determining the 
free-stream velocity, which used the fact that at the top 
of the profile the slope is a constant. This involved 
computing the cumulative average velocity starting from 
the top roost point of the profile which is well into the 
free-stream. The contribution of each point to the 
cumulative average is determined and compared to a 
tolerance. If the contribution is less than the 
tolerance, the contribution of the next point down the 
profile is computed and the tolerance is reduced by an 
amount proportional to the number of points in the 
cumulative average. When a point's contribution exceeds 
the tolerance, the value of the cumulative averaoe taken 
two points previous to that one was assigrred the value of 
• The sensitivity of the final result to the choice of 
the two parameters, namely, the initial tolerance and 
number of points to be stepped back after the tolerance 
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had been exceeded was investigated and used to provide 
consistent results. The values of ll*, from each 
neasurement station constituted the streamwise 
distribution of Uco.(x) . 

For a particular boundary layer case, the 
acquisition of the profiles, starting at the upstream 
most station, required data runs which lasted 
approximately six hours. During those runs, a very slow 
free stream velocity drift was detected. This resulted 
in the free-stream velocity for the last profiles being 
slightly lower than that for the first profiles. To 
compensate for this, so that it would not be interpreted 
as pressure gradient effect, the free-stream velocity 
time dependence was computed from the digitized output of 
the free-stream reference pitot-static probe and used to 
correct the free-stream velocity values. This involved 
increasing the It, values by an amount equal to the 
difference between the free-stream velocity read by the 
pitot-probe at the beginning of the data run and the 
average value over the time when each set of profile 
points was taken. The velocity values of the profile 
points were suitably increased by an amount proportional 
to their percentage of the free-stream speed so as* to 
preserve the profile shape. 

• If one is to use the mcmentum balance for 
determining the wall-shear stress, it is very important 
that any time dependent changes in the tunnel operating 
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conditions be accounted for. A great deal of effort has 
been spent over the last ten years quantifying the wind 

tunnel used in this investigation. Additional hot-wire 

0 

diagnostics and flow visualiation experiments were 
conducted early in this study. The excellent agreement 
achieved in regular boundary layers between this approach 
and several well accepted formulas, as discussed in 
Chapter VII, was a prerequisite to proceeding with this 
method in the canipulated boundary layers. 

To obtain a functional representation of the U^^Cx) 
distribution, the values obtained from the above technique 
were fit with a least squares cubic spline fit with 
weighting and damping. The advantage of the cubic spline 
is that it fits a smooth function through the data points 
rather than fitting the points tc some pre-described 
function. In the limits, a zero damped spline function 
will pass through all of the data points, while an 
infinitely damped spline will fit a straight line through 
the data. In its application here, all of the points 
were equally weighted and the deimping was adjusted to 
obtain a smooth distribution through the data values. 

After smoothing, the velocity values of the profile 
points were raised or lowered by a proportional amount so 
that the free-stream value equaled the spline smoothed 
value at that stream.wise position. The streamwise 
distribution of the dU^/dx was obtained from the direct 
differentiation of the spline function. The local dU^/dx 
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values were further smoothed by a spline fit for use in 
the momentum computation. 

Since we placed less confidence on the accuracy of 
the points taken closest to the floor we developed an 
iterative method for computing the quantities needed in 
the solution of the momentum equation. Initially, for 
each velocity profile, the datum point closest to the 
surface of the flat plate was chosen as an "initial 
guess" for y corresponding to y+ = 10. The velofcity 
profile was subsequently taken to be linear ro that heich 
and to pass through the data points at all higher 
elevations. The profiles were then integrated to obtain 
5^(x), 5^ (x) ar.d the energy thickness 6^(x) . In the 

process of integration, the profiles were smoothed by 
fitting a parabola through successive groups of three 
data points. Two data points of each group were 
overlapping. The integral values of the parabolic 
profiles taken across the overlapping regions were 
subsequently averaged and summed to obtain the total 
integral value for the complete mean profile. Given 
values for , the Ludwig— Tillmann relation was 

computed at this point to be used later for comparison to 
the values obtained by the momentum balance. 

For use in the momentum equation, a functional 
relation for was obtained from a least squares 

cubic spline fit. Discrete values of d: 2 (x)/dx were 
obtained from direct dif ferentiaticn of the cubic spline 
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function 


• The d52(x)/dx values were further smoothed by 
fitting the values with a cubic spline. 

The values taken from the "initial guess" were 
substituted into the momentum equation to generate new 
values of The values of were used to compute 

new values of y at y+ = 10 and the corresponding velocity 
values at that height taken from the linear velocity 
distribution 

+ + 
u = y 

Using these values the profiles containing new bottom 

points were again integrated to compute and d5^{^)/dx 

and used to solve the nonentum equation for • The 

solution continued to loop until consecutive values 

" 8 

differ by less than 10 . Convergence typically occurred 

after three iterations. 

To substantiate this method for determining the vail 
shear stress from the momentum equation, the boundary 
layer data taken by Wieghardt (1951) were used and the 
results compared to those obtained from other empirical 
formula's* The results of this comparison are discussed 
in Chapter VIII* 

In order to determine the wake component, Ti , for t:he 
outer part of the boundary layers, the transcendental 
equation 

1 2 - 
f(K) = i log[{-^^ -65)k] -I- B + ^ 

1 + n 

was solved using a Newton-Ra'phscn iteration r.ethcd. New 
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k dnd C were deterinined from the I09 portion of the 
velocity profile. In the case of the ncn~inanipu-ated 
boundary layers the standard values of k and C, 2.44 and 
5.0 respectively, were found to represent the data well 
and were used in the "law of the wake" solution. In the 
manipulated boundary layers, a least-squares curve fit 
was performed on the u+ vs, I09 y+ data to deter,.ine the 
"best" k and C values. These coefficients, which 
differed from the non-roanipulated case, were then used in 
the solution for n , 

One-dimens i onal Processing 

In order to obtain a statistical cescripticn of the 
digitally acquired turbulent velocity fluctuations r 
one-dimensional time-series digital processing v<as 
performed on the data. The processing ..ncluded 
calculation of the first four nomentSr probability 
density, one-dimensional Fourier transforms, energy 
spectrum, auto correlation and integral length scale, -n 
reaard to studying the wall layer phercmena, 
one-dimensional detection using the average local 
variance as a discriminator is presented along with 
one-dimensional match filtering metnocs for the 
velocity signature pattern recognition. The 
one-dimensional processing was performed on the time 
series velocity data of the single hot-wire sampled at 1^ 
KHz and of the individual rake sensors as well as the 


data from the shear wire. 

In the case of the 10 KHz sampled data, the tannel 
air temperature was not digitized along with the 
anemometer output. In order to compensate for 
temperature differences from the calibration reference 
temperature, the hand-recorded temperatures were used to 
correct the voltage values. Since the change in 
temperature results in an apparent change in the d.c. 
offset and gain in the anemometer voltage, the 
temperature differences were compensated for by modifying 
these two values which were input for the digital 
calibration of each data file. In order to recover the 
actual anemometer voltage value from the digitized value, 
the equation 

Sample value = {E - OFFSET) x GAIK x ADGAIN 
is solved where 


OFFSET = Analog d.c. bias 

GAIN = Analog amplifier gain 

ADGAIN = Analcg-Digital conversion factor. 

In correcting fcr temperature changes which result in 
changes in E, the processed, corrected voltage is £C 
where, 


, E(U,T) 
‘gain X ADGAIN 


OFFSET) 
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Sample value = (EC - OFFSET') GAIN’ x ADGAIN 

where 

OFFSET' = OFFSET/K(T) 

and 

GAIN' = GAIN X K(T) 

Therefore the gain and offset values could be adjusted 
file by file to account for changes in tenperature away 
from the reference. 

The one-dimensional spectra were obtained by taking 
the discrete real valued forward FFT of the 2048 point 
time series and constructing the one-sided power spectral 
density by multiplying the values in the frequency domain 
by their complex conjugate and summing symmetric points. 
The power spectra were normalized by the area under their 
distribution, which is the variance. The 
auto-correlation is obtained by taking the inverse FFT of 
the one-sided spectral function. The integral length 
scale was calculated from the area under the. 
auto-correlation function, where it is defined as the 
time lag intercept of a rectangle of equal area with a 
height equal to the zero lag correlation value, or 
variance. Since each time series constituted a number of 
integral time scales, their correlation functions should 
approach zero for long time delays. As pointed out by » 
Wigeland (1979), very low frequency fluctuations, which 
are not related to the flow characteristics, but rather 
to quasi-steady changes in the laboratory environment. 
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can cause the correlation functions to not asymptote 
toward zero. Since the integral time is computed from 
the area under the curve, this can cause significant 
errors. To remedy this, the method detailed by Wigeland 
(1979) was used to calculate the contribution of these 
very low frequencies and remove these effects from the 
determination of the integral time scale and variance. 

The probability density calculations are actually 
the number density distribution where 

F. = SAMPLE probabilitv that [ d. , < u < d.l 

and d^ specifies the width of a bin. Before sorting tht> 
numbers, the values were scaled by a iriable amount 
corresponding to 10/rms(u) in order to increase their 
significance and to make the bin width a function of the 
rms of the velocity fluctuations. In the distribution 
for the wall shear for each boundary layer case, 200K 
data values were sorted into 300 bins of width 1/10 
rms (r ) . 

"Burst" Signature Detection. In order to determine 
the velocity signature for the wall events, so called 
"bursts", the variable-interval time-averaging (VITA) 
technique introduced by Blackwelder and Kaplan (1976) was 
employed. This technique utilized the localized 
short-time variance defined as 

VAR(x^,t,T) = U^(x^,t,T) - [u(x^,t,T)]^ 

as a discririr.ator, where T is an averaging time which is 
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on the order of the event duration. The detection 
function is defined as 


f 1 if VAR < k.u^ 

D(t) = { 

I 0 otherwise 

where k is a threshold level and u is the 

rms 

long-time-average root mean square of the record of the 
series, namely, 


2 

u 

rms 



VAR 


In essence we are searching for correlated events which 
stand out above the background (Blackwelder , 1?77; 
Antonia, 1972). Except for detection purposes, the 
method has serious flaws. The first of these is that th« 
frequency of occurrence of events detected by this 
technique is dependent on the threshold and local 
averaging time, in addition the technique does not 
always detect the same -class* of events. The "burst- 
event has been classified as a local deceleration 
(ejection) of the streamwisc velocity followed by a rapid 
acceleration (sweep). The large local variance results 
from the large deceleration to acceleration sequence. 
Nishicka et al. (1980) show that although some of the 
latter class of events occur, the majority are of the 
deceleration/acceleration type. Therefore, the VITA 
technique only yields the signature of the most coherent 
event linked to a large local variance. It is with this 
application in mind that the technique was utilized here. 
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In keeping with the results of Blackwelder and 

Kaplan (1976), a threshold value, k, of 1.2 and a 

2 

non-dinensional averaging time, Tu^/v=10 was used. 

Also, the detection was made at a non-dimensional height 
of = 15, 

The detection was performed on the outputs fron the 
rake of hot-wires. In order to determine any Reynolds 
number effect on the frequency of events, the technique 
was also applied to the 10 KHz acquired outputs of the 
traversed single wire. The raw anemometer outputs were 
calibrated on the UNIVAC 11/08 computer and written to 
digital tape. These tapes, containing the instantaneous 
velocity values were then processed on the PDP-11/04 
mini-computer. The processing involved scanning each 
record of 400 points for a local variance greater than 
the threshold. The locations of the detected points were 
subsequently stored on floppy disks for later processing. 
The educted ensemble average of the velocity distribution 
at the time of the event was constructed by averaging 
each velocity time series centered around the detection 
point. In the case of the rake of wires, the pattern was 
educted not only at the detection wire but also for the 
ether sensors. In this manner, the ensenbied signature 
cf the. events for all of the boundary layer cases and 
Reynolds numbers were generated for use as a "pattern" in 
the one-cinensional match filtering which is described in 
the following section. 


One-dimensional Pattern Recognition, Given the most 


coherent velocity signature which was linked to a high 
local variance^ one-dimensional pattern recognition 
processing was done to obtain unbiased statistics for the 
frequency of occurrence of these events. This is 
contrasted to other "pattern recognition" schemes such as 
used by Wallace et al. (1972), which do not have the same 
statistical basis. The pattern recognition was 
accomplished by match filtering the time series of the 
detection wire with the educted velocity signature. The 
match filter has an impulse response equal to the 
time-reverse version ct the signal waveform to which it 
is matched so that the response is identical to an 
autocorrelation functicn. The peak in the 
autocorrelation, marking the match, occurs at a zero lag 
tire and is proportional to the total energy of the 
signal waveforms. It can be shown theoretically (Lynn, 
1973) that the improvement in signal-to-noise ratio 
caused by a match filter is the best possible. 

The match filter was constructed by centering the 
one-dimensional educted time series, obtained from the 
method of the previous section, ir. a 64 point array. In 
order to eliminate the circular effect in the time 
limited autocorrelation function, the tim.e series was 
extended to twice its original length with zeros. The 
real 128 point FFT of the time series was taken and 
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stored in the resulting 128 point complex array. In 
order to -orm the correlation functions, the imaginary 
part of the complex filter function was multiplied by -1 
to form the complex conjugate. 

The velocity time series of each record was broker, 
into six 64 point segments which were each time extended 
with zeros to 128 points. The real I’FT was taken of each 
segment and the real and imaginary parts of the filter 
function and data were multiplied (convolved) in the 
frequency domain. The inverse complex FFT was taken of 
each of the convolved segments to obtain the correlaticn 
function. In order that the zero time lag would 
correspond to the center of the lz8 point correlation 
function, the two 64 point sections were interchanged. _ 
The match filtered output was constructed from the 
filtered segments by the so-called overlap-add method 
which accounts for the fact that the linear convoluticc: 
of each section with the samiple response is larger tham 
the section length. As illustrated in Figure 18, 
adjacent segments were overlapped by an amount equal tr 
the filter width and the autocorrelation values in thir- 
region were sumjned to generate a smooth transition 
between segments. The first and last half filter wicta 
of points in the correlation function are filled with 
zeros and therefore discarded. The original record 
length of points was then retrieved. 

Each record vas analyzed to locate tr.e positive 
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correlation peaks corresponding to a match with the 
pattern for the "burst" event. The match locations for 

• - - i ■ . 

each time series were stored and the difference in time 
between peaks, computed using the sequence sampling rate, 
was used to compute the frequency statistics. 

Two-dimensional Processing 

The two-dimensional space-time processing of the 
outputs from the hot-wire rake was initiated_ to allow a 
comparison between the flow structures deduced from the 
velocity field and those marked by smoke and stored as 
digital im.ages. A number of the processing techniques 
described here are similar to those used in the 
processing of the two-dimensional image data (see Chapter 
V) and, therefore, some overlap will be encountered in 
the discussion. 

If we consider the digital acquisition of the analog 
outputs of the rake c.f hot-wires they can be represented 
by the sample function 

u(y,t) = u(y + qY(y), t + rT) 

where q and r are arbitrary positive integers. Here the 
data are sampled simultaneously in time at the A/D 
acquisition rate with a constant pence, T, and in space 
by the variably spaced, Y(y), sensor wires. We would 
like to reconstruct the discreetly sampled velocity data 
into the low-frequency wave number distribution analogous 
to the digital image. As will be explained, 
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low-frequency is a relative quantity which depends on the 
sampling rate. 

If we consider the sampling theorem for the discrete 
sampling of continuous data, which can be found in any 
number of texts on digital processing, (Lynn (1973)), 
Bendat and Piersol (1971 and 1980), Oppenheim and 
Schafer (1975) , it states that for data sampled at a 
frequency 1/T, the highest frequencies that* can be 
realized in processing is half of that value or 

= 1/2T = l/(j^ 

Here f^^is the well-known Nyquist frequency or folding 
frequency. Furthermore, the sampling theorem guarantees 
that a properly sampled band-limited signal can be 
reconstructed by a properly designed reconstituting 
^^Iter. The reconstituting filter must have the 
frequency response F(tj) such that 

M 1 ^ 

io M > ^ 

or in other words the filter response must guarantee no 
aliasing in the spectrum. In the time domain, this 
function would have an impulse response equal to the 
inverse transform of F(^) or. 


I (t) 


1 

2 - 


^ ( ju,)e- 


d^ 


Kt) 




j-t 
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sin u t 

I(t) = — -t° 
o 

If we convolve this function with our sample function, 

u^Ct) = u^(t) ^ sin cjpt 

CJ 0 t 

it can be shown (Lynn (1973)) that an exact 
reconstruction of the continuous time series is obtained. 
In practice we can reconstruct the continuous data 
function for a range of frequencies below the Nyquist; 
frequency by choosing a low-pass filter which band-limits 
above that frequency. Although a sin(t)/t function was 
used in this analysis, in practice it cannot be applied. 
In this application a finite impulse response filter was 
specifically designed to reconstruct the two-dimensional 
data space. 

Two-dimensional Filter Design . In the case for the 
space-time analysis of the rake of hot-wires, the 
sampling frequencies of interest are the A/D rate and 
wire spacings. If 'we consider first the A/D rate, we can 
choose an integral convection velocity to convert that 
frequency to a spatial passage of a coherent eddy. A 
typical value based on the acquisition rates which was 
used in the experiment is 2.8 mm. /sample. The largest 
spacing between sensors was 11.8 mm. so that the latter 
is the lim.iting value for the design of a two-dimensional 
low-pass spatial filter. 
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The two-dimensional FIR filter was generated from a 
64 point one-dimensional low-pass filter whose impulse 
response was symmetrically expanded into a 64 x 64 
two-dimensional plane. In the 64 point array, the 
maximum sensor spacing corresponded to 12 points so that 
for ,5f/f^ equal to the total 64 points, the maximum 
cutoff frequency for proper reconstruction was O.OSf. 

.The linear phase FIR filter was designed using the 
program written by Mclellan et al. which was one of a 
number of programs* for the processing of digital signals 
(IEEE Press, 1979). This program allows the design of up 
to 128 point, multiple stop-band filters with multiple 
weighting of ideal response characteristics. In this 
case the 64 point filter had a pass band from 0 to 0.05f/f^ 
and a stopband from 0.08 to 0.5f/f^ with an equal 
weighting of 5 for both bands. The weighting and edge 
cutoff frequencies of the bands were designed to give a 
minimum eutiount of ripple and the greatest amount of 
am.plitude cutoff. In this design the amplitude in the 
stopband is reduced by 60dB. 

The program returns the impulse response which was 
subsequently used to generate a two-dimensional symmetric 
response by R(x,y) = R(x) * R(y) . The 64 x 64 array was 
time extended to 128 x 128 by adding zeros to .the 'top and 
fight to ensure linear convolution filtering and the 12H 
point two-dimensional FFT was taken. The frequency • • 

domain representation of the filter was stored for later 
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use in the reconstructive processing. The one- and 
two-dimensional spectral function of the low-pass filter 
is shown in Figures 19 and 20. In the two-dimensional 
representation, only the positive frequency quadrant of 
the digital filter is shown. 

Two— dimes si cnal Discrete Fourier Transform . The 
discrete Fourier transform representation of 
two-dimensional sequences is of considerable 
computational importance in the digital processing of 
two-dimensional signals such as the digital photographs 
and in the reconstruction of the space-time velocity 
features in the boundary layer. In two dimensions, the 
Fourier transform is applied over a finite area 
corresponding to a block of data or part of an image. 
Because of the separable nature of the Fourier transform, 
(se% Oppenheim and Schafer, 1975) the two-dimensional DFT 
can be implemented by using a cne-dimensional DFT on the 
rows and then on the columns or vice versa. 

In its application, as shown in Figure 21, the 64 x 
54 blocks of data were time extended by adding zeros to 
form a 128 x 128 real array of which the two-dimensional 

W' 

FFT was taken. The 12S x 128 real-valued array requires 
32K words of storage and its transform at the worst case 
requires 65K words. Even with the 65K words of . MDP 
display memory, this would only allow one transform to be 
stored in RAX at one time and furthermore leave no rocrr. f 
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in .. eTiTiedi ate and final results of the processing 
sequence. If the- transforms were stored on disk, tie 
time required for disk access would increase the 
processing time by an order of magnitude. Therefore, ve 
utilized the symmetry properties of the Fourier tramsfcrm 
to reduce the storage requirements. The real part c± tie 
transform is conjugate symmetric and the imaginary part 
is conjugate antisymmetric. Therefore, we need only- 
store half of the transform and the other half can fc« 
reconstructed. The redundant transform domain samples 
resulting from the symmetry properties are shown in Fig- 
ure 21. It rust be noted that these symretrv prcpert:ies 
only apply to the Fourier transforms of real functions 
(Oppenheim and Schafer, 1975). In the case of the 
two-dimensional transform the first transform is done on 
a real valued data, the second is done on the result of 
that transform which is complex-valued so that the 
symmetry only applies in the first direction. The 
savings, however, reduce the storage from 65K words to 
32K words so that both the transformed two-dimensional 
filter function and the data block being transformed can 
be stored in the fast access HDP memory, in terms of 
computations, we also reduce the number of FFTs taken in 
the second direction by a factor of two. 

The 128 X 1^8 real array of data to be transform^ed 
initially resided in the upper or lower half of the HZ!P 
display memory’. The irst transform was done in rows by 
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reading the individual "sticks" of data into the MSP 
array processor and. taking the real one-dimensional 128 
point FFTs. The transforms were packed into 64 complex 
values by taking the symmetric half and placing the 
symmetry point into the complex part of the d.c. bin. 

They were subsequently stored in non-bitreversed form in - 
place in the KDP memory. The 64 columns of 128 word long 
complex- values were each transformed by reading in those 
"sticks" of dajta and perfcrming the complex 128 point 
FFTs. The 128 point transforms were written back in 
place in non-bitreversed form in the MDP memory. 

The inverse two-dimensional transform was performed 
by inverting the procedure in the forward transform. The 
64 columns of complex valued numbers were read into the 
MSP array processor, put into bitreversed order, and the 
complex inverse FFTs were taken. The columns were then 
replaced by their bitreversed transform, in place, in the 
MDP graphics memory. The 123 rows of complex values were 
each read from the graphics memory into a one-dimensional 
PDP array. The symmetric half was then reconstructed and 
the 128 point array was complex inverse transformed in 
the MSP. The real values of the result were taken and 
optimally scaled to an integer format whore they were 
written for storage and display purposes in the display 
processor memory. 

In order to reduce execution time, any. processing 
which was repeated a number of tim.es, for example the 
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symmetric reconstruction, was handled by specifically 
written assembly language subroutines. Also, whenever 
possible, the PDP, MD? and MSP processors were allowed to 
run in parallel. The total execution time for a forward 
or reverse 128 x 128 point FFT was approximately 30 
seconds . 

Two-dimensional Reconstruction . As with the 
one— dimens ional analysis of the rake data, the continuous 
time series was broken up into six segments, each 64 
points long. In the space direction the simultaneous 
data from the nine sensors furthest from the wall were 
placed in the array rows 0, 1, 2, 3, 4, 6, 12, 19, and 
31. The spacing? in the array represent their 
proportional spacing in the rake. Any further spacing of 
the sensors to acc'om.modate more of the 64 point space 
dimension was not possible because a suitably lew ripple, 
low-pass filter could not be designed. The array was not 
reducf^d to 32 x 32 because the reduction in points used 
in the design of the impulse response for the digital 
filter degraded the performance. The shear wire and hot 
wire at y+ = 7 were not used in the reconstruction. 
However, 'since the data is over-specified in that region 
of the array for the lew-pass filter cutoff designed by 
the largest spacing, their influence would have been 
minimal. Zero values filled the space between the rows 
where data values were specified. The arrangeifent is 
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shown in Figure 21. 


Since we wished x.o highlight flow structures ,* the 
velocity data for the two-dimensional reconstruction were 
normalized by removing the ensemble mean value and 
dividing by the ensemble rms of each of the sensory. We 
therefore intended no give emphasis to parts of the flow 
which stood out above the background turbulence or 

W ^ * 

"noise". These structures were assumed to remain 
reasonably intact as they were being convected past the 
probe. Furthermore, we expected that these structures 
would be the same as the coherent features visualized by 

smoke and stored as digital images. 

Each segment of data was read into the MDP memory 
and the two-dimensional FFT was taken with the special 
frequency representation being stored in place in memory. 
The convolution was performed by multiplying the real and 
imaginary pares of the column values of the transformed 
data block and filter function. The result then replaced 
the transformed data. The inverse two-dimensional FFT 
was taken of the convolved data and the result was stored 
away. As with the one-dimensional match filter, the 
segmented data were smoothly connected using the 
overlsp^add method* 

The result of the convolution is left in the center 
of the 128 X 128 array. Since these are segments in the 
time series, the overlap-add is done in that direction 
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only. As illustrated in Figure 21, the segments are 
overlapped by an amount equal to the filter width in the 
time direction where the values were added. The top and 
bottom half filter widths of the two-dimensional arrays 
were discarded along with the first and last half filter 
widths of each 400 point long continuous realization. 
Using the match points from the one-dimensional 
analysis, th« two-dimensional representations were 
ensembled about those instants in time to obtain an 
*^avera.ge* view of the "burst" sequence. Also the 
tvo-dimensional.,spectral function was calculated for the 
reconstructed data. In that case the spectra were formed 
by taking the two-dimensional FFT of a reconstructed 
segment of data. The transformed data columns were 
multiplied by their complex conjugate and the real valued 
function formed the positive and negative frequency 
spectra. The spectrum was normalized by the volume under 
the distribution, i,e., its two-dimensional variance. 
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chapt5:r V 

PROCESSING OF IMAGES 

The methods documented in- this chapter describe- the 
two-dimensional processing techniques which were applied 
owa.d systematically extracting information from the 
digitally acquired images of the smoke-wire visualized 
flowfield. These images, whose acquisition Was 
conditioned on a so-called "burst" event, are intended to 
provide visual information about any coherent flow 
structures link.ed to the process. This chapter will 
first discuss the methods used for the visual enhancement 
of the images and the special processing techniques to 
convert the images to a form better suited to machine 
analysis. Later, two-dimensional pattern recognition 
using match filtering will be discussed in order to 

extract and mark coherent features appearing in the 
images ^ 

Enhancer^ and Processing 
The image enhancement process cdnsisted of methods 
that sought to improve the visual appearance of an image 
convert the image to a form better suited to hum.an 
interpretation or machine analysis, in terns of visual 
appearance, one method of enhancement consisted of the 
applicatidn of pseudo-colors to an otherwise gray level’ 
image. Another application of image enhancement was tc 


I 



emphasize salient features of the original image in order 
to simplify the processing task* For example, the images 
could be high-pass filtered in order to emphasize the 
smoke-wire streaklines thereby simplifying streakline 
following routines like those used by Nagib et al* 

(1978) . 

The digital images were stored on a digital magnetic 
tape for post processing after acquisition. In order to 
reduce storage space, each pair of 8 bit pixel intensity 
values were packed into a single double byte PDF word. 
During processing, the images were "unpacked" and the 
pixel values were read into the display memory for 
viewing on the color monitor* The images consisted of 
white smoke on a black background and were recorded in a 
maximum of 256 gray levels* Therefore, an exact black 
and white reconstruction* was accomplished by loading the 
red, green and blue color lockup tables in the display 
processor with identical color intensity ramps from* 0 to 
255. In this way, a pixel value is translated into a 
gray level intensity proportional to its original^ 
luminance. After processing for image detection it was 
often helpful to display the results in peeudo-color 
which were designed to enhance the detectability of 
features in the image. The pseudo-colcr.mapping of the 
original image is. not a color imcage but rather a 
two-dimensional array of values that is converted to a 
color plane. In this application a color table consisting 
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of 36 bands of color hues and saturation levels was 
generated to represent the data. 

Although the digital-camera output values ranged 
over 8 bits, the display processor was capable of 
displaying unsigned integer values of up to 15 bits. 
Therefore we often would rescale the pixel values to 
occupy the full 32K integer range. The easiest way of 
doing this is to compute the maximum and minimum pixel 
values and to do a linear rescaling so that 


P 

F(i ~ niin 

^ J • ' t:* _ IT* 


F - F . 
max min 


X 32767 


This in effect performs a histogram expansion cf the 
original image and thereby increases its contrast. 

Because of the extreme light sensitivity cf the 
digital camera, the images were very sensitive to slight 
differences in background illumination. Differences from 
image to image could come about as a result cf 
non-consistent flash intensities from the strobe light or 
non-uniform vaporization of oil from the smoke-wire. 
Within an image, low-frequency variations in lighting 
arising from shadows or reflections or slight 
non-uniformities in the sensitivity of the Viticon tube, 
were imperceptible in a single realization but they were 
brought out as strongly coherent featu.res in the ensemble 
averaged images. Therefore, before. applying pattern 
recognition techniques to these images it was necessary 
to remove any non-uniformities which were not related to 
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the flow structures. 

In order to remove any stationary variations in the 
image luminance the pixel values were corrected by 
adjusting their values by an amount proportional to the 
"average" background values. The average background 
illumination was obtained by ensemble averaging all of 

I 

the digital images. The random nature of the smoke in 
the images canceled their contribution to the average so 
that only the average stationary background features 
remained. The- .value of each pixel in an image array- was 
changed by an amount corresponding to the average 
background pixel value at that array point. This can be 
written as 


F(j ,k) 


F(jA) 


where K is a coefficient of adjustment used to expand the 
result over a large usable ra.ige. In effect, this 
procedure digitally "dodges" the images in a more 
systematic way than can be done by hand in the dark room. 
This process was performed on all of tne images before 
pattern recognition processing. 

In some cases it was beneficial to remove a dominant 
feature from the images, such as the hot-wire rake, in 
order to better see the smoke around an otherwise hidden 
region. In the case of the sensor rake, this Was 
accomplished by recording an imag'e of the probe alone 
without smoke. It was necessary that this image was 
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perfectly scaled and its position was perfectly 
registered with that in the standard images^ The 
two-dimensional FFT was taken of the rake image alone to 
obtain its spacial frequency representation, A 
corresponding portion of a smoke plus the probe image was 
similarly transformed to the two-dimensional spacial 
frequency domain. The two arrays were then subtracted on 
a point by point basis to remove the frequency content 
associated with the probe. The inverse two-dimensional 
' FFT was taken of the result to leave the smoke image 
without the probe. 

Another method for removing or highlighting image 
features which are oriented with some directional 
coherence was a digital implementation of another dark 
room technique. The pixel values of a digital image are 
linearly rescaled to^utilize the full 15 bits available 
in the display processor and then a "negative" of the 
image was generated by the transfer function: 

F(j,k) = 32767 - F(j.k) 

The negative image was translated in a direction normal 
to the oriented feature of interest and the pixel values 

were averaged and rescaled to obtain a new pixel value 
F(j,k) = F(j,k) F(j + Aj, k + Ak) 

The resulting image brought out the edges of the features 
in the translation direction and everything else not 
coherent in that direction was left as a uniform gray 
background. 
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Two-dimensional high-pass and low-pass filtering of 
the images was performed by convolv4<ig the images with a 
digital filter function. This was done either by a 
discrete convolution of the form 

F{nj^, n^) H + 1, n'2~^2 


or by array multiplication in the spatial frequency 
domain. The criteria for the choice of methods was based 
on the size of the impulse response array, K(i, j) . In 
our case, if the dimensions of the filter array exceeded 
3x3, it became more computationally efficient to do the 
convolution in the frequency domain. The low-pass filter 
masks (taken from Pratt, 197R) that were implemented here 
are 


1 1 


1 . 


H = ^.;i 1 1 


and 
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l! 


^ =1-0 i' 2 r 

These are essentially noise cleaning masks with 
normalized weighting so that they would not introduce a 
brightness bias into the processed image. 

The high-pass filter masks were utilized to 
highlight edges in the images. The masks utilized in the 
discrete convolution were 


f 0 

-1 

0] 

'"1 

5 

-1 

1 0 
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and 
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These possess the property that the sum of their elements 
is unity. The filter masks were implemented in an 
interactive way to enhance visual interpretation of the 
images after processing. 


Another form of edge enhancem.ent that was used to 
highlight features of the match filtered results is 
called statistical differencing (Rosenfeld, 1969) . This 
involves the generation of an image by accounting for its 
statistical deviation from its neighbors. The transfer 
function implemented in the processing was suggested by 


Wallis (1976) 

G(j,k) = 


-• Ao^ 


t.ra 4 (1 -Q| F(j,kj) 




This function enhances the image by forcing a form on the 
first and second order central moments. In this case m. 

Cl 

and c , are the desired mean and standard deviation 

Cl 

factors. The coefficient A is a gain factor which 
prevents overly large values when is snail, and 

is a proportionality factor controlling the edge 
enhancenent. To apply this function, the processed image 
was divided into non-overlapping pixel blocks in which 
the block mean and standard deviation were computed. 
Bilinear interpolation of the four nearest block values 
was used to compute the mean and standard deviation at 
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each point. Along the image bordeir the mirror image of 
the blocks was used in the. interpolation. The block size 
influenced the size of features that were enhanced. In 
the case of the zoomed images in which we viewed 
approximately two boundary layer thicknesses in 128 x 256 
points- we utilized a block size of 16 x 16, 

Ensemble Averaging . Since the images were 

conditioned on the detection of a "burst" event, the 

first processing that was performed was to ensemble 

average the images. In this case not only the average 

image was computed, but also the rms of the differences 

between the images. Regions of low rms signified 

features that were coherent between realizations. The 

same procedure had been done for the phase conditioned 

images of the visualized wake of a cylinder (Nagib et 

al., 1975}. The mean and variance of K related digital 

images was formed in the usual way through the relations 
N 

1 1: [ F^(j,k) - M{j,k)] 2 

n=l 

and ’ 

I I = M(j,k) 

^ n=l ^ 

where is the pixel intensity at point (jfk) in 

realization n. 

As indicated in Chapter IV, the detection method 
used to acquire the images was based on a high local 
variance and w*as not alvrays a reliable indicator. In 
order to make the averages converge more quickly, we 
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included in the ensemble only those images which 
corresponded to an event detected by the more reliable 
match filtering technique applied to the hot-wire sensor 
outputs acquired at the same instant. 

Two-dimensional Spatial Filtering , As in the case 
of the two-dimensional Fourier processing for the 
hot-wire rake signals, spacial filtering was employed in 
the processing of the digital images. Just as the 

individual hot-wires spatially sampled the' velocity field 

. • 

in the boundary layer, the digital image is a 

two-dinensional spatial sampling of points of a physical 
image. In its application with smoke visualization, each 
sample point or pixel is a measure of the image intensity 
only. The pixel intensity is related in some 
non-linear fashion to the density or concentration of 
smoke particles which emanate from a line source and are 
convected downstream. In the boundary layer, the 
transport capability of the turbulence is represented 
(see Tennekes and Lurley, 1972) by the mo.Tientum flux -;uv. 
In a turbulent boundary layer, the three-dimensional 
eddies maintain a good correlation between the u and v 
velocity fluctuations. Indeed, the largest uv 
correlation values are associated with ^the wall "burst" 
events {Lu and Willmarth, 1973). Since the transport cf 
the smoke is related to the relative distribution and 
energy of the turbulent scales, we might expect that the 
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spacial distribution of smoke is an indication of these 
scales. Therefore, rather than looking at the relative 
intensities of the picture elements, we proposed to look 
at their relative groupings to infer information about 
the turbulence structures. This was be done using 
Fourier filtering techniques. 

In the discretely sampled images, the sampling 
frequency is the spacing between pixels. In line with 
the sampling theorem, the smallest scales that could be 
spectorally resolved are two pixels on edge. In the case 
of the zoomed view, the images corresponing to two 
boundary layer thicknesses are made up of 256 x 128 
pixels, so that scales as small as .02: are resolvable. 

In these images, however, we were mainly interested in 
larger scales so that in order to speed the computations, 
the images were bilinearly interpolated to 64 x 128 
pixels thereby cutting our resolution in half. 

In the case of the ensenble averaged images, in 
order to remove the streaklines and enhance the large 
scales, the images were lev-pass filtered. A number of 
64 X 64 Doint two-dimensinal low-pass filters were 
designed from a set of one-dimensional filter impulse 
response functions which were symmetrically developed 
over a two-dimensional rectangular area. The design of 
the impulse response utilized the same computer program 
developed by McClellan and used in the two-dimensional 


The filtering was implemented by dividing the 64 x 
123 images into two 64 x 64 blocks which were stored in 
HDP graphics memory. The area to the tcp and righw of 
the 64 X 64 image blocks was zero filled in order that 
the circular convolution will be identical to the desired 
linear convolution. The 128 x 128 two-dimensional FFT 
was taken of the two image halves and stored in place in 
the display memory. The FFT subroutine was identical to 
that used in the velocity reconstruction so that the 
reader is directed to that section of Chapter IV for a 
detailed description. The filter function which was 
already transformed and stored, was convolved with the 
image halves by complex multiplying columrrs. The 
convolved parts replaced the originals and the inverse 


FFT was performed to return the images to the spatial 
domain. The two image halves were smoothly connected 
using the cverlap-add method described in the previous 

chapter . 

In order to increase the dynamic range of the image 
domain upon reconstruction, the generalized Cepstrum was 
formed. This is a method of nonlinear transform image 
enhancement which consists of taking the logarithm of 


each real tran^forn ccefficienr. 


rhe irr.acinar/ part of 


the transform domain datus points were left untouched so 
as not to destroy the phase information. ' Qualitatively, 
this processing reduces the dynamic range within the 
transform domain and provides a stretching of the 
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resultant dynamic range in the reconstructed image. 

The two-dimensional spectrum of the images was 
computed by taking the FFT of the 64 z 128 images and 
multiplying the Fourier coefficients by their complex 
conjugates. The real valued spectrum functions were 
normalized by their volume, ensemble arranged and plotted 
as two-dimensional level contours with positive and 
negative**spatial frequencies. The volume under the 
spectra corresponded to the statistical variance of 
scales in the image and were saved for comparison. 

Pattern Recognition 

In order to search for structures in the visualized 
flowfield which were associated with rhe wall layer events 
and to obtain the frequency of these events from the 
images viewing 10 boundary layer thicknesses, 
two-dimensional match filtering was performed on these 
zoomed out records. The patterns used for the match 
were the ensemble averaged two-dimensional structure 
shapes obtained from the zoomed in images conditioned on 
the wall layer events. These two sizes of images were 
scaled properly to the same physical size by expanding 
the zoomed back views by padding each pixel with 
neighboring ze;o value points and low-pass filtering 
these images to reconstruct the intermediate! points. As 
in the one-dimensional pattern recognition described in 
the previous chapter, the two-dimensional application ^ 
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utilized a correlation match technique to determine i-f 
and where a match occurred. 

The ccaputaticr.al procedure consisted of 
transforming the 64 x 128 images into wave-number spa^e 
as described in the previous settion on spatial 
filtering. A 64 x 64 representation of the 

r • 

two-dimensional structures was similarly transformed. 

The imaginary part cf the Fourier coefficients of the 
matcl^ing function was r^ultipiied by -1 to form the 
complex conjugate and t:he colunns were conplex multiplied 
with those of the transformed image to form the trans- 
fcrr.ed correlation function. An inverse transforr; was 
taken of the result to obtain the two--dimensicnal spatial 
correlation function. Any peaks in the distribution 
corresponded to location where a match occurred. The 
pattern recognized images were later ensemble averaged to 
obtain the average visual representation of tne flowfield 
approachirfc the point vhere a wall layer event has been 


detected . 




CHAPTER VI 

PAPJ^ETRIC INVESTIGATION OF MANIPULATORS 

In the preliminary stages of this study we chose to 
move toward our objectives in a number of steps. The 
first of these was to verify that the outer intermittent, 
structure cf a turbulent boundary layer could be 
substantially removed. In this aspect we were concerned 
primarily vith altering specific scales of turbulence in 
the boundary layer. Initially, we were not overly 
concerned vith the device drag. However, minimizing the 
drag was part of the overall objective, so we chose a 
device which potentially could provide a minimum 
obstruction, when compared to the screens of Yajnik and 
Acharya (1177) , various types of honeycombs or even 
perforated plates. 

Early cn it was recognized that the initial design 
should allow some fle>:ibility based on a sufficient 
number of parameters sc that one can to match the device 
to the turtulence scales of interest. This is important 
for trying to sort out the complex interaction between 
the inner and outer scales of turbulence in the boundary 
layer. V^e were also interested in how turbulence 
generation near the wall is influenced by the outer 
scales, or lack of them. Therefore, the device had to be 
capable of directly modifying, in a controlled way, a 
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narrow spectrum of scales. An arrangement of stacked 
parallel plates suspended above the wall of the boundary 
layer at different elevations satisfied these design 
criteria. 

Four-Plate Manipulator ; Prel imina ry F.es'Jlts 

Our first attempt at removing the large-scale 
intermittent structures in a turbulent boundary layer was 
reported by Corke et al. (1980) . The turbulence 
manipulator utilized in that study consisted of four 
plateS/ vertically stacked in a parallel arrangement and 
spanning the tunnel test section. All cf the plates were 
0.7 mm. thick and had a streamwise chord of 75 mm. This 
manipulator was designed to remove a vice range of the 
scales in the turbulent boundary layer. Therefore, the 
manipulator plates were spaced successively closer at 
lower y values. The dimensions of this device relative 
to the boundary layer thickness at the streamwise station, 
where it was placed is shown in the top line of the table 
of Figure 22. Here the four-plate manipulator, being the 
first tested, is designated M-1. 

The effects of placing M-1 in the turbulent boundary 
layer was recorded by photographs of the smoke-ware 
visualized flowfield and diditally 'pcguir^ velocity 
profiles taken with a traversing hot;-wire'"probe. The 
photographs revealed the lack of the large-scale 
intermittent structure and a redistributing of fine 
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scales of turbulence in the manipulated boundary la^^ers. 
As will be described later in this chapter ^ these ' 
photographs will be used as a standard to judge the 
effectiveness of other simpler devices in removing large 
scale . 

The velocity surveys of this study recorded an 
approximate 10 percent reduction in* the value of u'/U^ 
below y+ = 100 when the M-1 device was introduced. The 
^streamwise distribution of the fraction coefficient, 
was determined from a momentum balance similar to that 
used in the current study. There were, however, several 
key differences between the ccrputational procedures used 
in the preliminary study and these used in the present 
study. 

As described in Chapter IV, the solution of the 
momentum equation required obtaining functional 
expressions for the streamwise distribution of and 

In the previous study, values cf LV were "fitted" 
with a second order polynomial, its derivative was 
subsequently fit by a linear function which was necessary 
to "smooth" the values of du/dx. The values of ^2(^5 
the preliminary work were fitted with a power-law 
equation as suggested from two-dimensional boundary laver 
analysis (Schlichtihg, 1968). The derivative of the 
pcwer-lcw function, which is used in the momentum balance 
is ever decreasing for increasing streamwise distances. 

In the simplified case of a zero streamwise pressure 
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gradient, the nomentum equation reduces to 


^ ; o < p < 1 


so that the strear.wiEe distribution of t directly 
follows the power-law distribution which the data 
prescribe. In the case of the regular boundary "layers 
the use of the power-law representation for the growth of 
the momentuin thickness provided good agreement with the 
accepted friction coefficients at comparable Reynolds 
numbers, however, for the manipulated boundary layers 
the limitations of the power-law fit resulted in an 
erroneous shear distribution from the momentum balance 
calculati ur. . Tr.e main reason for this error is the local 
influence of the manipulator which results in an 
initially higher momentum thickness. The resulting 
streamwise distribution is not well represented by a 
power- law function in this case , 

In the present results, the functional 
representations of U,^(x) and (x) we’-e obtained by 
fitting the data with a cubic spline w'th weighting and 
da.-ping. The formulation of *-he spline is given '^y 
Greville (1969) . The splin*. 'cximatas a continuous 
and differentiable function over an interval in a 
piecewise fashion using third order polynomials in 
non-overlapping subintervals formed by the data pairs. 

The spline function possesses continuous first and second 
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order derivatives over the full interval. In the limits 
of very high damping the cubic spline will "best fit" the 
data by a straight line. Contrasted with a standard 

I • 

least-squjres fit; which fits the data to a prescribed 
function; the cubic spline formulation fits the function 

to the data. In its application here, the damping was 

♦ 

adjusted to smooth the variations in the streamwise 
values of U..(x) and ^^{x) and their derivatives. The 
smoothed values were then used in the momentum* balance 
computation . 

The results of this preliminary study gave visual 
evidence that the large intermittent turbulent scales in 
the boundary layer wore removed by the four-plate 
manipulator. Near the wall {y+ < 100) the streamwise 
turbulence intensity was reduced by approxim.ately 10 
percent. However, as a result of the large momentum 
deficit produced by the device, no net reduction ir> the 
friction drag was observed. We therefore undertook a 
parametric study of a family of parallel-plate 
manipulators which could cor.parat ively suppress *the large 
scale while exhibiting a lower device drag. 

F am i 1 y of Parallel-Plate Manipulators 

Having documented the effectiveness oi a four-plate 
"turbulence man ipuia tor " for removing the intermittent 
large-scale eddies in a turbulent boundary layer (albeit 
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with an excessive device drag) we documented^ primarily 
through visualization records, the effects of a family of 
such devices that require a fewer number of plates to 
suppress the large scale to the same degree as that 
exhibited by M-1. This parametric investigation was 
aimed at developing a better understanding of the 
mechanisms by which the manipulators work and hence 
schemes for minimizing the device drag. This is an 
essential ingredient in the performance of such 
manipulators since our ultimate goal is to obtain net 
drag reduction for the turbulent boundary layers. The 
summary of results from this investigation is presented 
in the following. These experiments were performed prior 
to the extension of the wind-tunnel test section 
discussed in Chapter II. 

Parametri c Study. The manipulators utilized in this 
study were made up of one, two or four plates in parallel 
or series configurations as shown in Figure 23. In the 
single plate condition, plate angle of attack was also 
taken into consideration. The different parameters for 
each configuration, such as manipulator plate elevation 
in the boundary layer, streamwise spacing, plate chord 
length and angle of attack were adjusted to optimize the 
device effectiveness. This was done while observing the 
smoke-wire visualization in the boundary layer downstream 
of the manipulator. The relevent parameters were 
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systematically varied to obtain a flow which exhibited 
the least amount of large-scale turbulence. The final 
values obtained for this set of devices is tabulated in 
Figure 22. The flow conditions from these optimum 
configurations are recorded by smoke-wire visualization 
photographs and mean velocity and turbulence intensity 
profiles in Figures 24 through 28. It should be pointed 
out that for the sake of clarity, the scales of the 
photographic records and the velocity profiles in the 
conbined figures are not the same. In these figures, the 
distance y on the velocity profiles is approximately 
three times the geometric scale of the photographs. 

The results presented here represent the visually 
observed changes in the large-scale turbulent structure 
of the boundary layer, and the differences in some of the ^ 
mean quantities of the flow caused by the addition of 
each of the manipulators. Except for the cases with 
parametric variation of angle of attack, these results 
are for the devices which have been optimized with 
respect to the downstream boundary layer conditions. 

The results presented in Figures 24 through 26 
demonstrate the effectiveness of each of the devices in 
controlling the intermittant large-scale turbulence. Tc 
aid as a reference, these figures include at the top, a 
typical realization of the visualized boundary layer 
before the addition of a manipulator. The right portion 
of these figures contains a schematic representation of 
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the different device types with manipulator plates shown 
at their relative heights with respect to the mean 
boundary layer profile. The boundary layer profiles 
shown are those obtained in the unmanipulated boundary 
layer at the two positions where the devices are 
installed and where the comparison is carried out 
downstream. The visualized flow resulting from these 
devices is presented in the opposite left part of these 
figures. It should be noted that the streamwise 
dimensions of the device in the schematics of these 
figures are not drawn to scale so as not to obscure the 
profiles . 

In Figure 24, a comparison is made between 
manipulators having one, two or four vertically stacked 
plates. The four plate manipulator is the same device 
used in our early work. The boundary layer within which 
this family of devices was studied, however, had somewhat 
different characteristics from that previously used. 
Device M-1 was found to behave in an identical manner. 

If we compare the results from the two-plated (M-3) 
and single-plate (K-2) devices, in Figure 24, with those 
of M-I , we observe that the introduction of a single 
plate has some effect on the large scale eddies. 

However, at least two plates are necessary to approach 
the effectiveness of M-1. As expected, we found that 
these plates should reside in the outer portion of the 
boundary layer, where the larger scales are being 
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convected. 


In Figure 25 , the result of spacing two single 
plates in the streamwise direction is presented. In this 
case, a quite comparable reduction of the intermittant 
eddies, as with M-1, is observed. This configuration is 
quite appealing since the drag on the downstream plate is 
reduced by being in the wake of its upstream partner. 

Also, in terms of studying the effects of the lack of the ^ 
large-scales on turbulence production and the so called 
"bursting" phenomenon, the placement of this device high 
in the boundary layer limics any direct influence of 
their wakes on the flow near the wall. 

Figure 26 demonstrates the effect of the angle of 
attack of a single plate on the intermittant large-scale. 
Depending on the sign of the angle used, the plate 
generates a wake dominated by eddies which rotate in the 
sam.e sense (negative angle of attack) or opposite sense 
of direction to the large structures in the boundary 
layer. In the former case this tends to augment the 
large scales, instead of suppressing them. Another 
mechanism of large scale suppression is similar to that 
active in honeycombs and documented in earlier studies at 
IIT (Loehrke and Nagib, 1976); i . e . , the inviscid 
inhibition of the transverse velocity fluctuations by the 
solid boundaries. 

As was demonstrated with manipulator H-1, a measured 
reduction in the streamwise turbulence intensity in the 
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boundary layer is found below y+ = 100. In Figure 27, 
the turbulence intensity profiles for devices M-1, -3 and 
-4 exhibit the same characteristic reduction in that 
quantity near the wall. Also included in this figure are 
the mean velocity profiles for the same cases. The 
velocity defects in the profiles with the devices give an 
indication of the device drag when compared to the 
profile of the unmanipulated boundary layer. Both 
manipulators with only two plates, either in parallel or''" 
in series, 'result in less device drag compared to the 
case of M~l, 

In order to extend our measurements to farther 
downstream distances from the manipulator in the same 
fetch of wind tunnel test section, we generated a 
turbulent boundary layer approximately one-third as thick 
as that used in previous cases. Using a properly scaled 
two-plate manipulator of the M-2 design, we observed the 
flew to approximately 7C boundary layer thicknesses 
downstream of the device. The result of that experiment 
are summarized in Figure 28. As far as we can judge witih 
the visualization and the u* profiles, the effectiveness 
of the manipulator seems to persist for the entire 
downstream fetch. 

As a result of this simple parametric study of 
various flat-plate manipulator configurations, we have 
verified that an optimum exists for the placement of 
single or multiple plates in parallel or series 
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arrangements. The optimum is based on minimum device 
drag for the comparable effectiveness in suppressing the 
large-scale turbulent structures in the boundary layer. 
The results suggest that two plates spaced in the 
streamwise direction^ i.e.f in tandem/ give a reduction 
in large-scale intermittant eddies in a turbulent 
boundary layer which is comparable to a stacked, 
four-plate manipulator with a much less cost in device 
drag. With the newly extended wind-tunnel test section, 
we utilized these conclusions to establish the decree of. 
net drag reduction that may be achieved by such passive 
devices in turbulent boundary layers as well as 
documented some changes in the mechanisms of turbulence 
production in the wall region. 
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MEAN CHARACTERISTICS OF TURBULENT B0UN:3ARY LAYERS 

In order to demonstrate the changes which occurred 
in the boundary layers as a result of the suppression cf 
the intermittent outer scales, the time-meBn and unsteady 
characteristics of the naturally developed turbulent 
boundary layers, into which the manipulator was placed, 
were* first documented. These measurements not only serve 
for comparison to the manipulated boundary layers, but 
also act as a bench mark for comparisons to the classical 
literature and as a reference for the reader to judge the 
quality of the present measurements. 

This chapter depicts the mean and turbulent velocity 
distributions of these boundary layers for the range of 
momentum thickness Reynolds numbers from 2000 to 5500. 

This range was achieved by operating at two 
different free-stream speeds of 6,7 and 10.4 m/s. The 
"burst" conditioned m.easurements for these layers will be 
presented in Chapter IX where a more direct ccmpariscr. 
with the "manipulated" boundary layers is aiade. For the 
figures to be referenced in the following chapters, we 
used a Cartesian coordinate system, that has its origin at 
the leading edge of the boundary layer plate on the 
center plane of the test section. The x, y and z axis 
are in the strear.v ise , vertical and transverse directicr.s, 
respectively . The data were taken at 10 screamwise , 
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stations along the plate centerline {z equals zero) as 
depicted in the schematic of Figure 3. The figures in 
this chapter will typically show data from alternating 
stations for clarity purposes^ although the processed 
data exists for all of the measurement locations. 
Furthermore, all of the digitized raw data tire series as 
well as the piccessed results are stored on digital 
magnetic tape to be used for any further processing. 

The non-dinensionalized mean velocity and turbulence 
intensity profiles fot the boundary layers at different 
streamvise stations for the two free-stream speeds is 
shown in Figures 29-32. The vertical, y, distribution of 
the mean velocity in Figures 29 and 3C is 
non-dir.ensionalized vrith the local free-stream velocity, 

, and the momentum thickness V7e chose the outer 

variable -2 non-dirr.ens ionalize the vertical distance 
rather than the boundary layer thickness, because the 
former is an easily determined physical characteristic of 
the flew unlike the somewhat vague definition for the 
boundary layer thickness. Whenever possible, 
non-dimensionalization with outer variables will be in 
terms of the momentum thickness. Exceptions wrill be 
those cases where comparisons are made to the results of 
other investigators which involve the boundary layer 
thickness as a parameter. For the readers reference, 
the ratio of the boundary layer thickness to momentum 
thickness is approximately 10. The rms of the 
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fluctuating velocity distribution for the boundary layers, 
presented in Figures 31 and 32, is non-dimensionalized 
with the local wall-friction velocity . The values 
for were determined from the moraentuir. balance and will 
be discussed later in this chapter. The friction 
velocity was chosen to non-dimensionalize the fluctuating 
streanjwise velocity distributions because it provided a 
collapse for the full vertical extent of the boundary 
layer. Non-dimensionalizing this quantity with the outer 
free-stream velocity would only provide a collapse of the 
downstream development of the fluctuating velocity at the 
outer porticfn of the vertical distribution. This was 
expected since the turbulence stresses are linked to the 
production at the wall and proportional to the friction 
velocity u^. 

Utilizing these non-dimensional forms, the 
streamwise development of the mean velocity and 
turbulence intensity for the two ranges of Reynolds 
numbers show a high degree of collapse that demonstrates 
their self-similar behavior. Included for reference in 
Figures 31 and 32 is the turbulence intensity 
distribution taken from the results of Klebanoff (1954) 
for a Reynolds number of 7500. Comparing the turbulence 
distributions in these two figures we observe that the 
higher Reynolds number flew case results in a slightly 
greater distributicn of turbulent stress fluctuations 
away from the wall and thereby a better comparison to the 
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higher Reynolds number flow case of Klebanoff (1954) . 

The values for the wall-shear velocity used in the 
previous figures were obtained by a balance of the 
momentum flux along the boundary plate. As a result of 
the extremely mild streamwise pressure gradient^ which is 
shown in non-dimensional form in Figure 33 f the 
functional form for the momentum balance included the 
streamwise rate of change of the total pressure. The 
pressure gradient is mainly a result of the pressure 
field generated by the perforated plate at the end of the 
test section which r .aintained the sligr.tly r.igher pressure 
on the measurement side of the plate and ensured parallel 
flow at the leading edge. In this respect the perforated 
plate acted like the flap often used to control the 
circulation on test boundary layer plates. The reduced 
streamwise growth exhibited in the manipulated boundary 
layers had no measurable effect on the streamwise 
pressure gradient for the two operating free-stream 
velocities. The consistent agreement and calculation 
approach allowed us to maintain comparable accuracy in 
using the momentum, balance for the determination of sk.cn 
friction in regular and manipulated boundary layers. 

In order to test the software program for 
calculating the values of wall-shear stress from the 
momentum, balance/ the data of V.'iegnardt (1944)/ as 
presented by Coles and Hirst (1968) / were utilized. 

Since in the latter reference a num±>er of the parameters 
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for this data were calculated, including the wall-shear 
stress, this offered a good comparison of the results 
from different conputational techniques. In particular, 
for the determination of the friction velocity, the 
values obtained from the empirical log-relation and the 
Ludwig-Tillmann relation were compared with the results 
from the mcnentum balance,. The number of profiles and 
their streamwise spacing was chosen from the wiegHardt 
data so as to be representative of the two parameters in 
this study. In the momentum balance the damping or the 
spline fits to the distributions of dV^/dx and do^/dx 
v;ere adjusted tc obtain the "best" correlation of rricticn 
coefficient values with those provided by Coles and Hirst 
(1968). The point of best agreement was also found tc 
coincide with the best collapse of the mean velocity 
profile values in the leg-linear region. Varying the 
spline dampings also provided us with a m.easure for t.te 
dependence of the results on these parameters. 

The results for the calculation of the fricticn, 
coefficients for V'ieghardt's data (1944) are shown in 
Figure 34. The distribution of the friction coefficient 
values for the range of Reynolds numbers in this figure 
show a good collapse of the results from the 
Ludwieg-Tillmann relation and the momentum calculcticn. 
Below a Reynolds number of 6C00 the comparison is not as 
satisfactory. However less emphasis was placed on 
results from the the fit of the log portion of the 
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profile because the sitial!^ nuntber of data values taken in 
that region (typically 2 to 3) added to the uncertainty 
of the fit. Furthermore, since the data of Wieghardt 
(1944) was used in part to obtain the Ludwieg-Tillman 
form for the friction velocity, the degree of agreement 
between their relation and the momentur. balance was taken 
as a more satisfactory indication of the performance of 
the program. In addition to the friction velocity, the 
integral properties of the boundary layer of Wieghardt 
were com, pared to the results published by Coles and Hirst 
(1968) and found to be in close agreement. 

After satisfactorily exercising the momentum, 
calculation with the Wieghardt data (1944), the profile 
data from, the present experim.ent were analyzed. The 
results for the distribution of. the local friction 
coefficient versus Reynolds num.ber are presented in 
Figure 35. The values for presented in this figure 
were obtained in three manners-, a least squares curve 
fit of the log portion of the profile from. y+ equal to 50 
to 200, the Ludwi eg— Til Imann relation, and the solution 
of the rconentum, equation. Also shown on this figure the 
corresponding distribution from. Wieghardt’ s (1944) data 
ever a comparable Reynolds num.ber range. The data is 
presented with individual symbols corresponding to each 
of the 10 streamwise stations for the two ranges of 
Reynolds numbers. 

An im.portant feature of the data is that three of 
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the stations have overlapping Reynolds "number values 
based on the momentum thickness. In the lower 
free-stream velocity, this corresponds to stations at the 
downstream end of the plate and in the higher free— stream 
velocity case to points at the upstream end. Any 
property of the flew which exhibits only a Reynolds 
number dependence should present a smooth distribution 
through this overlapping range In the distribution of 
Figure 35, all three methoc's gave values vithin 5 
percent of a mean distribution drawn through the points. 
As a result of the scatter of points in this figure that 
average was chosen to represent t.'e results. 

The values of the friction velccity taken the 
results i'n Figure 35 were used in the 


non-dimensi .alization with inner variables of the mean 
velocity distributions presented in Figures 36 and 2? for 
the two Reynolds number ranges. Drawn for comparison in 
these figures is the log— law relation using the 
coefficient values suggested by Coles (1963) . In both 
figures, the collapse of the data in the log region, 
50<y'*'<200 , is quite good. It shouid be pointed out that in 
this approach we have determined the value . of the . 

von Karman constant and the log-linear intercept constant 
from the moment'um balance. This in contrast to the often 
used approach of determining the skin friction velocity from, 
these' two constants, i.e., by a fit or tne log-linear 
profile. T n i s result is signi — — can^. i** u..cxw » 
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substantiates our application of a careful momentum 
balance for determinincr the wall shear stress. The 
agreement with the law-of-the-wall relation is, however, 
not very satisfactory in these figures. The lack of 
similarity in this range is most likely a result of the 
accuracy of positioning the traversing probe in the thin 
sublayers of these boundary layer flows in air. 

The wall-unit non-dimensionalization of the 
fluctuating velocity distributions is shown in Figures 38 
and 39 for the two ranges of Reynolds numbers. As 
indicated by these figures the approximate similarity in 
wall units of the rms values of the streamwise velocity 
fluctuations is maintained only out to an y+ equal to 15. 
A similar trend with Reynolds number was presented by 
Purtell, Klebanoff and Buckley (1981). They suggest that 
the decrease in the u'/u, values at lower Reynolds numbers 
is a result of the suppression of all but the largest 
scales of turbulence. The maximum in the distribution 
occurs at y+ values between 15 and 20. It was at y+ of 
15 that Blackwelder and Kaplan (1976) observed the 
largest effect on their "burst" conditioned averages. 

This region of the flow field has been pinpointed by 
Corino and Brodkey (1969) as the site for the ejection of 
low-speed fluid near the wall during a "burst" event, and 
as contrihut ir.vT as much as 60 percent of the total uv 
Reynolds stress (Willmarth and Lu, 1972) in the boundary 
layer . 
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We focus on the wake portion of the boundary layer 
profile in Figures 40 and 41. in this figure the wake 
component 

W(y/5)= ^ u - (i log y"*^ *+B) 

is plotted versus the non-dimensional boundary thickness 
y/i . The respective values of k and B of 0.41 and 5.C 
presented in the log distributions of Figures 36 and 37, 
were used in this formulation. For reference purposes 
wake distribution suggested by Coles (1956) is included 
in these figures. At the lower Reynolds number range the 
agreement with Coles' wake distribution is very 
satisfactory. At higher Reynolds numbers in Figure 41, 
the agreement is not as good below half the boundary 
layer thickness. Close inspection of the wake function, 
w(y/5), reveals that it is not entirely antisymmetric 
^^bh respect to y/c== 0.5. Actually we expect (see Hinze, 
1975) , deviations from the sin^ distribution at the top 
and bottom of the profiles. In the representation of 
Figures 40 and 41, the deviation from the antisymmetric 
distribution appears to display some Reynolds number 
dependence . 

Normalized one— dimensional spectra of the streamwise 
velocity component at different self-similar heights 
along alternating measurement stations is shown in 
Figures 42 to 47 for both Reynolds number ranges. The 
power spectral energy in these distributions is 
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non-dimensionalized by the area under them which is equal 
to the .variance of the instantaneous velocity 
fluctuations. The frequency axis is normalized by the 
local mean velocity at each respective height in the 
boundary layer to obtain the wave number, k. Also 
• included in these figures for reference purposes are the 
lines with slopes of -1 and -5/3 corresponding to the 
viscous-convective and inertial subranges, respectively, 
for one-cimensinal velocity spectra. The spectral 
distributions exist for the non-dimensional heights 
between those shown on the-fee figures but were omitted for 
clarity purposes. However, they did reinforce the trends 

exhibited in these figures. 

At any of the streamwise stations, the changes in 
the energy spectrum with non-dimensional height reveals 
that the contribution to the turbulence energy in the 
low-wave-number range, which is made up by larger eddies, 
decreases as the wall is approached. However the spectra 
indicate that the large scale, low-wave-number energy 
decreases only slightly down to y+ equal to 30, thereby 
indicating the influence of the large eddy structure very 
near to the wall. In a compatible way, as the wall is 
approached, the contribution to the turbulence energy m 
the high-wave-number range. is increased. 

In the outer part of the boundary layer, y+ greater 
than 200, the -5/3 law is followed rather closely over a 
discernable wave-number range. In the case for tne 
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lowest Reynolds number in Figure 42 the wave-number range 
having a -5/3 distribution correspondance extends 

from approximately K equal to 50 to 200 m With 
increasing Reynolds number we observe a growth in the 
extent of the inertial subrange to encompass a wave 
nut±ier range between 5 and 200 m Nearer to the wall, 
y+ less than 200, the spectrum shows a range where E(k) 
varies nearly according to k ^ , indicating ’a strong 
interaction between the mean and turbulent flow. It is 
at these heigtits in the boundary layer where we also have 
strong production of turbulent energy. Comparing the 
extent of the viscous-convective subrange in these 
spectra indicates no discernable Reynolds number 
dependence . 

The vertical distribution of the integral length 
scale at alternating streamwise stations for the range of 
Reynolds numbers is show’n in Figures 48 and 49. The 
length scales were determined by multiplying the integral 
time scales, obtained from the one-dimensional 
autocorrelations, by the local mean velocity at each 
respective streamwise location and height in the boundary 
layer. That is, obtaining these from a local application 
of Taylor's hypothesis. In Figures 48 and 49 the 
integral scales and vertical y positions are 
non-dimensionali zed with the outer characteristic length, 
5^ . These distributions show that the maximum length 
scales occur at a height approximately 20 percent of the 
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boundary layer thickness. This suggests the presence of 
turbulence structures which are elongated in the 
streantwise direction and which reside fairly rear the 
wall. This interpretation is consistent with the 
space-tiir.e correlation measurements of Grant (1958) and 
Blackwelder (1970). Near the wall, the streamwise extent 
of the turbulence structures is observed to lengthen with 
increasing Reynolds numbers. 

The integral length scales are not observed to 
behave in a self-similar manner with Reynolds number when 
non-c’ir-ens.ional ize:. vitli outer- variables as in Figures 48 
and 49. To evaluate the growth with Reynolds number of 
the scales at different self similar heights in the 
boundary layers, the data are presented in the manner of 
Figure 50. In this figure the integral length scales, 
which were non-dir.ensionalized with the momentum 
thickness, are plotted versus Reynolds number along 
constant y+ lines. The actual data were "best fit" by a 
linear distribution and for clarity purposes only these 
lines are presented on the figure to show the trends. 

The typical correlation for the fit was approximately .6 
- .7, where 1.0 represents a perfectly correlated 
numerical fit of the function to the data. The data as 

presented in Figure 50 can be interpreted as the 

« 

streamwise growth cf the average scales along self 
similar paths in the boundary layer. 

With increasing Reynolds number, near the wall below 
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y+ equal to 30, the turbulent scales are observed to 
decrease with streamwise position. Between y+ values of 
50 to 200 the scales are not growing appreciably with 
increasing Reynolds number. This corresponds to the 
log-region in the mean velocity distribution. At y-^ 
equal to 400 we see a rapid growth of the streamwise 
scales with Reynolds number. This most likely reflects 
the growth of the constant stress layer. At higher 
elevations, y+ = 800 , there is a continued growth of the 
streamwise scale although at a less rapid rate than in 

ff 

the initial wake region of the mean profile. 

To summarize this chapter, these boundary layers, 
developed over two ranges of Reynolds numbers, exhibit 
mean and turbulent characteristics which are comparable 
to similar flows reported in the literature. The 
streamwise development of these quantities also 
demonstrates the desired self-similar behavior. 
Importantly, they display a full spectrum of turbulent 
scales which follow the expected wave number 
distributions. Other statistical quantities providing 
further documentation of these turbulent boundary layers 
will be presented in Chapter VIII where they will be 
directly compared to the results fromi the m>anipulateld 
boundary layers. 
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MEAN CHARACTERISTICS OF MANIPULATED BOUNDARY LAYES^S 

The time-mean velocity measurements in the 
manipulated boundary layers are presented in this cb-apter 
so as to document the effects of suppressing the 
intermittent outer structures on the mean characteristics 
of the flow. This chapter depicts their mean and 
turbulent velocity distributions for the range of 
momentum thickness Reynolds numbers from 2400 to 49CD 
corresponding to the two free-streem speeds of 6.7 and 
10.4 m/s. In addition to the Cartesian coordinate system 
utilized in the results presented in Chapter VII, 
additional laboratory coordinates, which are related to 
the position of the manipulator, are introduced in this 
chapter. The coordinate length is the streanwise 
distance between the manipulator and the leading edce of 
the boundary-layer plate and denotes the downstre^ 
distance of a measurement station from the trailing edge 
of the manipulator. The data were taken at nine 
stream.wise, x, stations along the plate centerline, z 
equals zero. This is in contrast to the data runs 
without the manipulator in which measurements were 
recorded at 10 streamwise locations. For the results in 
the boundary layer with the turbulence manipulator, the 
first measurement station falls between stations 2 3Jid 3 
of the data runs without the device. In the manipulated 
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layers data staV.ion 2 corresponds to station 3 cf the 
other case and the remaining streamwise stations 
correspond to their N + 1 equivalents. As in the 
previous chapter, the data will typically be presented at 
alternating stations although the processed data exist 
for all of the measurement locations. 

The turbulence manipulator utilized in this final 
set of measurements was oriented in the series 
arrangement, M-4 of Figure 23, and optimized for the 
boundary layer characteristics depicted in Chapter VII. 

As described in Chapter VI, the optimization involved 
varying the different parameters to obtain the maximum 
suppression of the intermittent outer scales in the 
boundary layer. The degree of large-scale removal was 
judged visually with the aid of the smoke-wire technique. 
In this final set of results, the manipulator was placed 
at the streamwise distance, x^=2.25 m., corresponding tc 
the x-location of measurement station 1 in the regular 
boundary layer case. The manipulator plates were made 
from 0.2 mm. thick brass shim stock and each was placed 
at a height from the wall of 80 percent of the boundary 
thickness in the regular boundary layer at station 1. 

The streamwise chord length was equal to their distance 


from the wall which w; 


0 . 8 . 


i.ie strea. wise spacina cf 
the manipulator plates was 8 ; at that location. The 

arrangement of the manipulator plates is represented in 
the schematic of Figure 3 and shown in the photegraph of 
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Figure 4 • • 

The boundary layer characteristics r visualized using 

the STnoke-“wi re f before and after use of this series 

manipulator are shown in the photographs in Figure 51, 

The top two photographs correspond to the regular 

boundary layer and the bottom two correspond to the 

boundary layer with the final series manipulator* The 

f r ee^st r earn velocity at which these photographs were 

taken corresponds to the low Reynolds number range with 

Re ranging from approximately 28CO to 37Q0 over their 
6 c 

streamwise extent. Although these photographs represent 
one instant in the turbulent flowfields, their overall 
characteristics are representative of features of the 
boundary layer with and without the manipulator. 
Quantitative analysis of digitized image records of the 
visualize'^ flow for the two types or bounoary layers will 
be pres( ed in Chapter IX. The photographs in Figure 51 
allow corrpariscn to the visualized flows presented in 
Chapter VI, which were obtained in a boundary layer with 
different character isti*cs than that utilized in this 
final set of results. The photographs of Figure 51, as 
well as those discussed in Chapter VI indicate that if 
the large scales are being reformed downstream of the 
manipulator, their return process is extrer.ely slow. In 
Ficure 51 no strong evidence of their rebirth is found as 
far as far as more than 40 boundary layer thic'xnesses 
downstream. 
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Qualitative examination of the photographs in Figure 
51 clearly reveals the suppression of the larger outer 
bulges in the turbulent/non-turbulent interface with the 
addition of the series manipulator. Close inspection of 
the top two photographs in this figure and other related 
photographs reveals regions of potential flow, visible as 
intact streaklines, extending well into the bulk of the 
boundary layer and marking the intermittent non-turbulent 
portions of the flow. This mechanism is clearly 
inhibited by the removal of the outer scales as evidenced 
in the bottom two photographs. As will be demonstrated 
later in this chapter, this effect reduces the 
entrainment of ’ rrotational fluid into the boundary 
layer . 

The non-dimensionalized streamwise mean velocity 
distributions over a range of momentum thickness Reynolds 
numbers from 2400 to 4900 are presented for the 
manipulated layers in Figures 52 and 53. For comparison, 
a line corresponding to the non-dimensionalized velocity 
distributions of the boundary layer without the 
manipulator is also included in these figures. Also, the 
manipulator height, normalized with the characteristic 
length of the boundary layer at the first . streamwise 
station, is shown for reference in the figures. 

The wake of the manipulator plates is' evident in the 
profiles approximately 36 momentum thicknesses downstream 
of the device in Figures 52 and 53. At the next 
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downstream station, 5 /O 2 approximately equal to 136, the 
turbulent velocity fluctuations have smeared out the r^iean 

wake profile of the dev’ice. 

In the lower 50 percent of the mean profile taken 
closest to the device, the non-dimensional distribution 
looks like that for the regular boundary layer. In the 
outer part of the distribution, > 6, the lack of the 

larger scale is immediately evident by 36 momentum 
thicknesses downstream,, as indicated by the higher 
velocities at the edge of the layer with the device. The 
lack of the large scales is felt in the lower 50 percent 
of the boundary layer by 136 momentum thicknesses 
downstream where we observe a self-similar distribution 
having lower mean velocities than in the regular boundary 

layer. 

The distribution of the rms of the streamwise 
component of the velocity fluctuation, non-dimensional ized 
with the friction velocity u^ , is shown in Figures 54 and 
55 for the r.ianipuiatec. layers. The streanwj.se 
distribution for u. was obtained from the momentum 
balance in the manner described in the previous chapter. 
That distribution will be presented later in this 
chapter. Included. in these figures is the distribution 
taken from Klebanoff (1954) which was shown, in Figures 
31 and 32, to represent quite v.ell tne fluctuating 
velocity profiles in the regular boundary layer. 

Clearly, the streamwise turbulence fluctuations do not 
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display the self similar behavior exhibited by the 
regular boundary layer. Nearer to the deviccr up to 136 
momentum thicknesses downstream, we observe a reduction 
in the streamwise velocity fluctuations except in the 
vicinity of the wake where higher fluctuations, 
suggestive of a vertex street, result in locally higher 
rms values. Further downstream, higher turbulent 
stresses are observed to be filling the boundary layer up 
to approxirately 6C percent of its thickness. This 
redistribution of turbulence stresses away from the wall 
is reminiscent of the Reynolds number dependence observed 
in the low Reynolds number range of the regular boundary 
layer. In that case, the distributions in Figure 31 
demonstrated lower u*/u^ values in the constant stress 
layer portion of the profile when compared to the higher 
Reynolds number case of Klebanoff (1954). Near the 
floor, at the edge of the sublayer, the maximum 
values are approximately 16 percent higher in the 
manipulated boundary layer than in the regular bcondary 
layer. The maximum u'/Loo values in the manipulated, 
boundary layers are, however, reduced as presented in the 
preliminary work on these manipulated boundary layers ty 
Corke et al. (198C) . This is also demonstrated in the 
results from a parallel study which is presented in 
Appendix B of this manuscript. 

A comparison of the streamwise growth of the 
integral thicknesses, ^^9^2 ^3 t^etween the regular 
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and manipulated boundary layers for both Reyholds number 
ranges is presented in Figures 56 and 57. Included for 
the momentum thickness distribution, f^Cx), are the lines 
corresponding to the spline function fit used in the 
momentu , calculation for the friction velocity, u, , The 
figures demonstrate an initial thickening of the boundary 
layer with the manipulator. However, the streamwise rate 
o’f change of those quantities is less than in the regular 
boundary layer case. After approximately 25*0 momentum 
thicknesses, based on 5^ station 1 which corresponds 
to the leading edge of the manipulator, the momentum and 
energy thickness distributions cross their counterparts 
in ‘the regular boundary layer. The displacement 
thickness requires somewhat farther distance for the 
distributions to cross. 

Tne rar’o value? of the integral lengths are shown 
in Figures 56 and 55. In. Vic-ire 5:, she ratio of the 
aisplacem.er.s tr rrcr.’entun’. thicknesses, is presented 
for a range of roiienturr. thic.kr.ess Reynolds numbers 
in the regular and manipulated boundari layers. 

As demonstrated by Purtell et al. (1981', the shape 
factor for the regular boundary layer is solely dependent 
on Reynolds number, .Re^. . This, however, is* not the case 
for the boundary layer with the manipulator. As 
evidenced by the distributions in Figures '56 and 57, 
the displacem.ent thickness is disproportionately larger 
t; an the mcnentum thickness in the manipulated layers. 
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This is opposite to the effect of adding large surface 
roughness which resulcs in a reduction in the 
displacerr.ent thickness (see Hinze, 1980, Figures 7-15). 
As will be further supported in the following parts of 
this chapter, this has led us to conceptualize these 
manipulated layers as boundary layers with ’negative" 
roughness . 

The ratio of^the energy thickness to momentum 
thickness, ^ 2 ^' 2 ' plotted versus the shape factor in 
Figure 59. Included are the data from Wieghardt with the 
integral values given by Coles and Hirst (1968) and 
obtained fror. the present ccnputational scheme. Also 
included in Figure 59 are the lines cor respondina to the 
analytic expression relating with and given by 
Wieghardt and Tillmann (1941) 


which was obtained fr oir. the ve 1 oc i ty law U/C = ( v / - } 

■ ■ These distributions are presented for the ccnstar.rs. given 

i by Wieghardt and Tillmann (1941), which they had adjusted 

I agreement to their data, as well as constants 

i obtained by a ’best fit" of the relation to Kieghardt's 

as well as the present data. The shape factors clearly 
follow a unique relationship of the form obtained from 
the analysis, including the man-ipulated boundary layers 
However, over this range of values, new constants were 
I necessary to best represent the present data as veil as 

I 

I 
I 
1 


129 



* .> 


ORIGINAL FACS IS 
OF POOR QUALITY 

■c 

Wieghardt's data which was, in part, the basis for the 
Wieghardt-Tillmann {1941) constants. 

The Reynolds number distribution of the wall 
friction coefficients, r for the manipulated boundary 
layers is shown in Figure fiO. The line in the figure 
designates the cor respo'hding distribution for the regular 
boundary layer. The open symbols in the figure indicate 
the values obtained from the momentum balance and the 
half-open symbols indicate the values of obtained from, 
the log-linear and Ludwi g-Ti 1 Imann relations. In the 
latter two cases, the same coefficients as those 
presented for the analysis of the regular boundary layer, 
were used to calculate the streamwisc distribution of the 
wal 1- shear stress . 

It is clear from. Figure 60 that the shear stress 
distributions obtained from, the analytic expressions are 
not valid since neither demonstrates the necessary smooth 
Reynolds number dependence. This behavior was, however, 
anticipated and prompted the use of the momentum, 
calculation to determine the wall shear stress. As was 
stated earlier, there is no reason to expect that either 
of the m.ethods used in Figure 60 other than the mcnentum 

balance v:ork in the raninu la tod layer. They have been 
included here only for purpose of comparison and 
illustrations. In the manipulated case, the last point 
from the low Reynolds number range is the furthest 
streamwise station from the manipulator while its 
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datum point corresponds to the streamwise 
station closest to the manipulator in- the high 
free-stream velocity experiment. 

Comparing the local C£(x) distributions between the 
regular and manipulated boundary layers in Figure 60 we 
observe an approximate 30 percent reduction in the local 


drag on the wall as a result of the suppression of the 


large intermittent outer scales. This reduction in drag 
is evident in the momentum thickness growth which was 
shown in Figures 56 and 57. Here we observed a decrease 
in the rate of change of 6^(x) which, in absence of a 
pressure gradient, is proportional to the local 
values. That decrease in the slope is approximately 30 
percent, which compares well to the reduction in the 
magnitudes of obtained from the momentum balance which 
had included the contribution of the extremely mild 
pressure gradient. 

The net reduction in the drag must include the 
device drag of the manipulator. This can be obtained by 
calculating the net loss of momentum which, because of 
our near zero pressure gradient, is 


Percent Reduction 


^^2,9^1 - <^2.1Q^R 

^^2,10 ■ *^2,1^R 


By this definition, we obtain a net reduction in the wall 
drag of 21 percent. The 52^^) distributions, however, 
are still diverging at the last streamwise station so 
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that greater net reductions nay be achievable at larger x 
dimensions. 

The values of (x) were obtained from the momentum 
calculation and used in the non-dimensionalizaticn with 
inner variables of the mean velocity profiles which are 
shown in Figures 61 and 62 for both Reynolds number 
ranges in the manipulated boundary layers. It is evident 
from these figures that a log-linear region exists in 
these profiles. However^ the coefficients of the 
analytic expression have changed. Drawn for com.parison 
is the solid line in the figure which corresponds to the 
log-linear relation with the coefficients, 1/K = 2.44 and 
B = 5,0, which were given by Coles (1968) and found to 
represent the boundary layer profiles without the 
manipulator. The dashed line is the "best fit" of ti]e 
log-linear relation to the data over the range 50£y"*^£200. 
This fit yields the coefficient values, 

1/k =2.87 
B = 6.53 and 7 . 56 

over the range of Reynolds numbers for the manipulated 
layers. 

If we draw on the analogy of the increasing wall 
roughness which results in a downward shift in the log 
linear distribution cor responding to a reduction in the 
value of the coefficient B, the boundary layers without 
tne large outer scales em.ulate the qualities of a 
hypothetical "negative" roughness flow, Rotta (1950) 
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connected the value of B with the thickness of the 
viscous sublayer for large values o'f u^y/v to obtain the 
form 


1 

B = ^ (leg 4k -1) + — I— £ 


where i, is the sublayer- thickness. Therefore, an 
increase in the value of B can be associated with a 
thickening of the sublayer. Using the values for k, B 
and u, for the regular and manipulated boundary layers 
the sublayer thickness is estimated to have increased by 
approximately 17 percent in the manipulated boundary 
layers. The value of y+ corresponding to 6^ is 6,S, 

The thickening of the sublayer has also been 
observed in boundary layer flows with drag reducing 
polymers (Donohue et al., 1972; Achia and Thompson, 1977; 
Hanratty et al., 1977), The mechanisms, however, are 
clearly different since the drag reducing polym.ers work 
by inhibiting momentum transport through viscoelastic 
action near the wall. 

To account for the change in the value of the 
von-Karman constant in the manipulated boundary layers we 
consider the Reynolds number dependence of that quantity 
which was examined by Tennekes (1966) in turbuient pipe 
flow. In that reference, Tennekes obtains a second-order 
friction law and a Reynolds number dependent correction 
to the "universal" slope of the logarithmic velocity 
profile. This correction takes the form 
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1/K = 2.95 - 4.4 

If we take R* to be iu /u we obtain values of 1/K of 2.55 
and 2.87 for the middle ranges of Reynolds numbers in the 
regular and ma^nipulated boundary layers respectively. 
Although this theory has generally been unaccepted by 
others (see Yaglcm, 1979) it predicts the observed change 
in the value for k in this case with surprising accuracy. 

It is important to realize that the log-linear 
portion of the velocity profile and its coefficients were 
verified and presented in all cases with values of 
taken from a momentum balance. Therefore, no a priori 
bias was applied to the results except that momentum must 
be conserved. Thus the logarithmic region was found to 
exist in the boundary layers both with and in the absence 
of large outer turbulent structures. Based on the 
persistence of the logarithmic region down to m.omentum. 
thickness Reynolds numbers of at least 465, where the 
larger scales of the turbulent fluctuations become 
increasingly dom.inant, Purtell, Klebanoff and Buckley 
(1981) suggest seme form, of Reynolds number independence 
of the large-scale structures. The results of the 
present investigation do not appear to support this 
conclusion since the manipulated boundary layers lack 
these large intermittent scales. Another explanation is 
that the large scales are a result of the transition 
process to turbulent flow and are actually composed of 
other sm.aller scales such as inclined hairpin eddies 
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(Klebanoff, 1962) which are generated through a wall 
instability process and which are sustained even in the 
fully turbulent flow. The logarithmic region vrculd 
reflect the character of these eddies and, because they 
are related to a wall instability, would display some 
Reynolds number dependence. Further evidence toward 
this concept will be presented in results to follow’. 

The wall-unit non-dimensionalization of the 
fluctuating velocity distributions for the manipulated 
boundary layers is shown in Figures 63 and 64. Comparing 
these figures to the distributions for the regular 
boundary layer in Figures 38 and 39 indicates that the 
peak in the fluctuations in the streamw^ise velocity is at 
a lower w^all unit in the manipulated layers. In the outer 
part of the profile there also appears to be a much 
sharper cutoff of the streamwise fluctuations near the 
edge of the boundary layer in Figures 63 and 64, The 
wake function for these manipulated layers is presented 
in Figures 65 and 66. For this caser the ”best fit” 
values of k and B taken from log-linear distributions in 

Figures 61 and 62 were used in the determination of v:. 

2 

Included in these figures is the sin distribution 
suggested by Coles (1960) . Although the wake strength/ 

7., was determined by prescribing that W(y/5)=2 at y=:, 
the peak in those values is observed to occur at y/: 
approximately equal to 0.8. This would indicate that the 
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wake in the manipulated boundary layers is associated 
with structures which scale with an_ intermediate Iengt:h ^ 
smaller than the boundary layer thickness. 

The Reynolds number dependence of the wake strength 
parameter, II , is shown in Figure 67 for regular and 
manipulated boundary layers. Besides the dependence on 

Re^^ displayed in the fiaure, it is also a function cf th 
streamwise pressure gradient. It was proposed by Coles 
(1968) that the strength cf the wake component is 
and that it might be an indicator of the local rate cf 
production of large-scale turbulent energy. This 
quantity, non-dimensionalized with the free-strean 
velocity, is plotted versus momentum thickness Reynolds 
number in Figure 63. In this figure we observe a 
reduction in this quantity in the manipulated boundary 
♦layers. There is, however, a . discontinuity in the 
distribution between the two Reynolds num*ber ranges in 
the manipulated case indicating a dependence on the 
initial conditions imposed by the device on the flew. 

Another indication of the strength of the large 
outer scales is to consider the velocity of entrained 
fluid into the boundary layer. In a self-preserving 
boundary layer, the velocity of propogation of turbulent 
fluid into the f ree-st ream , namely the entrainment 
velocity, Vp, is 


'■’p = h 

which was shown by Bradshaw in his 1967 paper. Ir was 
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also found empirically by Bradshaw et al. (1967), th<=- 

/ r-2 

the quantity, Vp/Uco is very nearly proportional to 
both in equilibrium and non-equilibrium boundary layers. 

In that reference, they also imply that the ratio of the 
large-eddy velocity fluctuation to the small-scale velocity 
fluctuations in the central part of the boundary layer is 
proDortional to Furthermore, the ratio of 

advection or diffusion to production in the outer part cf 
the layer is shown to be proportional to 6V^/~ and tn.erec\ 
proportional, by the previous form, to 

In Figure 69, the ratio, Vp/Uoc is plotted versus He^ ^ 
for the regular and manipulated boundary layers. The tcp 
most dashed line represents the distribution of the 
friction coefficient, = i/h-L' « the regular 

boundary layer. This distribution was fit to the data to 
.obtain the proportionality constant between the entrainnent 
velocity and the non-dimensional shear at the wall. In 
this figure 


^ X 10^ = 0.86 (Cj X 10^) -1.3 


is found to provide a smooth distribution through the 
data values. 

In the manipulated boundary layers, the suppres.=»ic_ 
of the large outer scales results in lover propcrticnal 
entrainment velocities. The solid line in this figure 
corresponds to thi distribution for the manipulaued 
layers when scaled with t.he same proportionality ccustan 
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as in the regular layers. In this case the line does not 
represent a best fit to the V^/u^, values indicating the 
proportionality constant is not exactly the same for both* 
cases. In the case of the manipulated boundary layer, 
the relation 

^~P X 10^ = 1.16 (C.:XlO^) -1.6 

IV "- 

best represenis the data values. The average 
difference between the latter distribution and that 
obtained using the coefficient values for the regular 
boundary layer is approxinmately 20 percent, '“he average 
reduction in the entrainnent velocity in tne manipulated 
boundary layers is observed to be approxir.ately 40 
percent. 

The results of Figure 6 9 show sor.e interesting 

trends. The first of these is that the entrainment 

velocity is becoming smaller with increasing Reynolds 

number. This would indicate that the strength of the 

large-outer scale turbulence structures is decreasing 

with increasing . The ratio of diffusion to 

production in the outer part of the layer is decreasing 
2 

like T/stpc = 2C^ in the figure. This is a result of the 
decrease in the ratio of the large-eddy velocity 
fluctuations to the small-scale velocity fluctuations 

I- 

proportional to (2C^) Thus the outer part of the 
boundary layer is filling with smaller scale turbulence 
at higher Reynbolds numbers that are more quickly 
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dissipated by viscosity. By removing the large outer 
scales in the boundary layer ve have hastened this 
process. This is evident in the visualization records of 
the manipulated boundary layers which appear to be filled 
with finer scales of turbulence than in the corresponding 
regular boundary layer cases. This gives some evidence 
that the manipulated boundary layers emulate higher 
J^eynolds number layers. 

The wave^number distribution of energy in the 
manipulated boundary layers is shown in Figures 70 to 78. 
In Figures 70 and 71, the development of the frequency 
associated with the vortex street generated by the 
manipulator plates is presented. At the closest 
strearwise station, a peak corresponding to a strouhal 
frequency of approximately ft/U, = 0.03 is evident for 
both Reynolds number ranges. The scale of eddies 
associated with this peak is approximately 0.5 cm. The 
peak is observed tc be less dominant at the higher 
Reynolds number. Ey the next measurement station, 
approximately 17 plate heights or approximately lOf 
downstream of the device, we observe no spectral peak, 
corresponding to a regular vortex pattern, in the 
one-dimensional spectra. Other differences in the 
spectral distributions are a result of taking the spectra 
at a constant y-value equal tc the device height rather 
than along self-similar y-values in the boundary layer. 

The normalized one-di ’.ensional spectra of the 
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streamwise velocity component at different self-similar 
heights along alternating measurement stations is shown 
in Figures 72-76. Included in these figures are the 
lines with slopes of -1 and -5/3 corresponding to the 
viscous convecti’.'e and inertial subranges, respectively, 
for onc-diTTiensicr.al velocity spectra. The correspcncino 
spectral distributions for the regular boundary layers 
were^prcsented in Figures 42 to 47. 

For both Reynolds number ranges in the manipulated 
boundary layer ve observe a better- collapse of the 
spectral distributions taken at different heights at the 
station just downstream of the device (Figures 72 and 
75). This would indicate morG even distribution of 
turbulence scales across the boundary layer. Furtherncre 
we observe that each of the profiles nas a range of wave 
numbers over which they follow a -1 distribution. In the 
regular boundary layer y only the spectra taken very near 
to the wall (below yw = 30) exhibited = -I distribution. 
This would indicate that the vi scous-ccnvective 
interaction, and therefore greater turbulence production,- 
is occurring ever a greater portion of the boundary' lay^er 
close to the manipulator. Indication of the inertial 
subrange is only' evident in these figures at y £ 

In the regular boundary layer the -5/3 portion in the 
spectral distribution could be observed down to y+ = 2CC. 

At the farthest downstream station from the 
manipulator for the low Reynolds number range in Figure 
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74, we still observe a fairly good collapse of t±e 
spectral distributions for wave numbers greate'r the- IOC 
at y+ values less than 200. At a comparable Reynolds 
number in the regular boundary layer, the inertial 
subrange was evident in the «spectra taken at y+' values 
greater than 200, where it ranged between wave nonbers cf 
50 to 200 m ^ In Figure 74, the spectral distributions 
demonstrate a -5/3 slope down to y+,= 75, The wave 
number rang.' corresponding to the inertial subrange in 
this case is a much lower 5 to 50 r"* A k"^ distribaticn 
is observed in this figure for a range of wave numbers 
between approximately 8 to 100 at y-^ = 30. 

At the highest Reynolds number, shown in Figure 77, 
the inertial subrange is observed to expand to encoimass 
wave numbers between 5 and 100 Tnis is still an 

order of magnitude lower in wave numbers comparec to the 
regular boundary layer at a comparable Reynolds number. 
This would indicate that the constant stress layer has 
extended over a larger portion of the boundary layer. 
There is still evidence cf the visccas-convective 
interaction near the wall, below y+ of 75 at the bicher 
Reynolcs number. This would indicate that the events 
occurring near the wall in regard to turbulence 
production are reasonably unchanged in the boundary 
layers without the large outer scales. 

The one-dimensional spectra of the wall shear, t.. , 
for the regular and manipulated boundary layers at tvo 
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Reynolds numbers are shown in Figure 78. The 
wave-numbers are computed in this figure by using tfte 
average value of the shear velocity, u. , for each case. 
These distributions show a small range of wave numbers 
over which a energy charge is evident. The emaining 

distribution exhibits a behavior indicating strong 
viscous dissipation of turbulent energy. In the 
manipulated boundary layer, which corresponds to the 
lowest Reynolds number in Figure 11 , the wave-number 
distribution of the normalized spectra follows more 
closely the high Reynolds number distribution. This is 
especially more evident for lower wave-numbers, less than 
1000 m ^ in this figure. 

A ccmparison of the normalized probability 

distributions of the instantaneous shear stress for the 

« 

manipulated and regular boundary layers is shown in 
Figures 79 and 80 for both Reynolds number ranges. In 
these figures, the distributions are normalized by 
removing the sample mean and dividing the distribution 
axis by 0.1 times the sample rms. Therefore every 10 
bins correspond to one rms increment from the mean. This, 
normalization provides an excellent collapse of the 
distributions even* though the mean shear at the wall and 
the rms of its fluctuations are significantly reduced in 
the manipulated boundary layers. As is evident from 
these figures, the distribution of the surface shear 
fluctuations is skewed so that the most probable value is 
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less that the mean, a similarly skewed shear stress 
distribution was presented by Sandborn (1981) for a fully 
developed open channel flow. The probabilrty 
distribution is, however, Reynolds numbe Jependent, as ■ 
shown in Figure 80 where the skewness ha decreased by 
approximately 50 percent for an increase in Re^ of 
approximately 30 percent. Also with increasing Reynolds 
number, the t'/t shear stress intensity is observed to 
decrease by approximately 17 percent. In the manipulated 
boundary layers at a comparable Reynolds number, - ’ 'Y is 
higher by 22 percent over the regular boundary layer 
case. This was similarly true for the u'/u. values near 
the wall in the distributions in Figure 54. However, it 
should be pointed out that the absolute value of the 
wall s-.e^r iluctuations is reduced when the boundarv 


la\er is manipulated. As the Reynolds number is 
increased in the manipulated layers we observe a 25 
percent decrease in this value, which is about half of 
the decrea-se- demonstrated for the regular boundary layer. 

The shape of the distribution, indicated by the 
kurtosis, is also Reynolds number dependent, in the case 
of the regular boundary layer, the kurtosis is observed 
to decrease by approximately 20 percent. More 
importantly, it has moved from a value highef than that 
for gaussian to one less than for gaussian noise. 
Therefore the shear-stress distribution is becoming 
narrower at higher Reynolds numbers indicating a more 


143 


spiked time series f^r the shear. There is no significant 
change in the kurtosis values for the manipulated 
boundary layers. 

The St reamwise-velocity skewness and kurtosis 
distributions across the boundary layer for manipulated 
and regular boundary layers are shown in Figures 81 and 
82. Th^^skowness distribution in Figure 81 shows a 
maximum positive value between y+ of 1 and 7. This 
corresponds to most probably velocities less than the 
distribution mean, which would be associated with fluid 
moving away from the wall. Between y+ of 50 to 200 , the 
skewness is approximately constant and near zero 
(gaussian) . This corresponds to the constant stress 
region associated with the log-linear velocity profile. 

In the outer part of the boundary layer, y+ > 200, the 

t 

skewness becomes largely negative. This indicates that 
the most probable velocities are greater than the mean. 
Similar statistics were presentee by Gupta and Kaplan 
(1972) . 

Near the wall in the kurtosis distribution of Figure 
82, we observe a peak kurtosis value in the sublayer 
between y+ of 1 and 7. In this region the kurtosis is 
greater than three indicating a probability distribution 
which is narrower than gaussian and a velocity time 
series which is more spiked in nature. This could be 
categorized as indicating rapid parcels of low velocity 
(S greater than 0) fluid moving away from the wall. This 
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type of motion has been classified as the ejection phase 

of the bursting process. 

In the buffer region between y+ of 20 to 50 we 
observe the minimum kurtosis values in the velocity 
distribution. The values correspond to a broader 
distribution than in a gaussian function. These 
statistics indicate that relatively large (less spiked 
signal) parcels of fluid are moving through this region. 
The near zero average skewness indicates that these 
parcels are an even distribution of locally low and high 
velocity regions. 

Between y+ values of 50 to 200 we observe a nearly 
constant average kurtosis with a magnitude slightly less 
than for a gaussian signal. Again this corresponds to 
the region in the velocity distribution of constant 
stress. In the outer part of the boundary layer, the 
kurtosis distribution acoin attains values araater than 
that for a oaussian signal. This is associated with the 
spiked signal indicative of rapid, streamwise narrov 
regions of higher velocity fluid approaching from the 
outer regions of the boundary layer. 

The greatest differences that occur in these 
statistics between the regular and manipulated layers is 
in the regions below y+ of 50 and above y+ of 200. In 
other words, statistic similarities in the stream.wise 
velocity fluctuations are only unaltered in the 
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log-linear region of the velocity profile® 

The vertical distribution of the integral length 
scales in the manipulated boundary layer at different 
downstream stations is shown in Figures 83 and 84 for 
both Reynolds number ranges. These results can be 
compared to the integral length scale distribution for 
the regular boundary layer shown in Figures 48 and 49. 

r ^ « 

Comparing the figures reveals a growth in the streamwise 
scales near the wall _in_ the manipulated layers. The same 
characteristic lengthening of the streamwise scales with 
increasing Reynolds number is exemplified in the 
manipulated layers in the lower half of the boundary 
layer. At the first two measurement stations in the 

outer part of the manipulated boundary layer, the lower 

* 

integral scales are indicative of the vortex street 
emanating from the device. 

The distribution of the integral length scales with 
Reynclns number is shown in Figure S5 for the manipulated 
boundary layer. This is comparable to the results for 
the regular boundary layer in Figure 50. Near the wall, 
'.iOlrv y+ = 30, the streamwise scale is observed to remain 
nearly constant with Reynolds number. In the regular 
b^’^nfary layer, the non-dimensional length scale 
deerrased by a factor of two^ In the log-linear region 
of the manipulated bor.ndary layer in Figure 65, the 
ncr.-dim.en.sicr.al streamwise scale is observed to increase 
with Reynolds number. In t^e regular boundary layer in 
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Figure sO, the streamwise scales demonstrated almost no 
growth in the log-linear part of the velocity profile. 

At y+ = 400, the ncn-dimensional integral length scale is 
increasing approxinately 30 percent faster in the 
manipulated layer than in the regular case. In both 
cases, the non-dimensional streamwise scale at y+ = 800 
increases at a slower rate with increasing Reynolds 
number, tTiaii those at y+ = 400. In the manipulated 
boundary layer, the growth rate of the scale at that 
ncn-dir-ensional height is approximately twice that in the 
regular boundary layer. 



CHAPTER IX 


UNSTEADY CHARACTERISTICS OF WALL LAYER 

The results presented in this chapter correspcnd to 
the analysis of the "burst" conditioned velocity 
measureirents obtained from the rake of hot-wires and the 
simultaneous digiti;;ed smoke-wire images of the flew. 

The techniques used for the acquisition and processing of 
these data are detailed in Chapters III, IV and V. The 
rake of wires, including a surface mounted shear wire, 
was located at measurement station 9 or the equivalent 
station 8 in the regular or manipulated layer, 
respectively. The mean characteristics of the boundary 
layers at those stations have been documented by the 
results in Chapters vii and VIII. The m-' n 
characteristics of the boundary layers s streamw’ise 

station are also summarized in Table 2 this chapter. 
The results will induce one-dimensional (time) analvsis 
and "burst" detection, two-dimensinal (space/time) 
reconstruction of the rake outputs and analysis of the 
digital images. In all cases, the results from the 
regular and manipulated boundary layers will be compared 
to draw light on the physical processes occurring ir. the 
near wall region of these turbulent boundary layers. 
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Table 2. Mean Boundary Layer Characteristics at Rake 
Location r 


BL 

(m/s) 


, : r 3 ( cm> ) 

(m/s) 

C 2 

C ^(xlO^) 

1 e 


R 

6-7 

1.231, 

.881, 1.548 

.265 

3819 

3.026 

.0628 

. 67- 

R 

10-4 

1,118, 

.811; 1.430 

.378 

5473 

2,659 

.0741 

.91' 

M 

6.7 

1.117, 

.783, 1.371 

.218 

3364 

2.136 

.0628 

.73: 

M 

10.4 

1.026, 

.732, 1.285 

.316 

4933 

1.852 

.0741 

.794 


One-Dimensional Analysis of Velocity; Time-series Data 

The outputs from the 11 hot-wire sensors can be 
treated either as 11 individual time series or as a 
two-dimensional space-time array of velocity data. In 
this section we will first analyze the rake outputs in 
the former manner. The latter case of the results of a 
space-time analysis will be presented in the following 
section. 

Figures 86-89 display a time-series record of the 
digitized outputs from the 11 hot-wire rake sensors for 
the low and high Reynolds number ranges of the regular 
and manipulated boundary layers. The time axes of these 
figures are noQ-dim.ensionalized by the outer variables 
and * . The spacing between sensors and their distance 
from the wall are denoted by their self similar 
coordinates based on inner-wall characteristic quantities 
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and V. The tine-series values are presented as the 
instantaneous fluctuating quantities honnalized by the 
record average rns. This was equivalent to the long tine 
average mean since the record lengths were chosen so that 
the mean statistics were ergodic. The central moments 
within a data record correspond to the long-time average 
moments, at least to the first four orders. Other than 
for the sensor at y+ = 0, the quantity plotted in these 
figures is 

(Q(y) - C(y) )/Q(y) Q = [ ^ : y"^ > O 

rir.s 1 ^ ^ . y+ ^ 0 

For the remaining time series at the wall, the velocity 
in the above expression is replaced by the wall shear 
stress, This method of ijon-dimensionalization is 

necessary for presenting these time series in this form 
because their mean values are increasing while the rms of 
their fluctuations are decreasing at greater 
ncn-dimensional distances from, the wall. Also, since we 
are interested in focusing our attention on the coherent 
'’str uctures" in the turbulence, this quantity represents 
a measure of the fluctuating intensities which stand out 
above the background turbulence fluctuations. In 
particular we were looking for the coherent "signatures" 
of these structures sc as to obtain statistics about 
their distribution in the boundary layers. 

The time series in Figures 86-89 correspond to the 
first 315 discrete data values in one record of continuous 
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data. As described in Chapter Illf the data records in 
the low Reynolds number range were acquired 
simultaneously with the digitized images of the 
smoke-wire visualized flow. The instant in time at which 
the images were taken corresponds to point 50 of the 315 
points displayed in Figures 85 and*88 and is denctec by 
an arrow on the time axis of these figures. The 
triggering was based on a locally high variance in the 
velocity signal for the wire at y+ = 7. Inspection of 
the time series at that non-dimensional distance frcm the 
wall in Figures 86 and 88 reveals a locally high 
excursion in the signal which would resi:lt in a high 
variance. Inspection of the time series at other 
distances from the wall at that instant in time reveals 
similar features on a time scale of apprcximateli* 2-35/u^. 
These features appear coherent^ over an acrrerlabie 
vertical extent of the boundary layer, including the 
time-distribution of the wall shear stress. As a first 
step in understanding the physical process depicted by 
these time series, analysis of the velociry distribution 
near the wall is first presented. 

This phase of the experiment was designed to make 
comparisons between the regular and manipulated boundary 
layers over a range of Reynolds numbers. In particular 
we were interested in determining what effect, if any, 
the suppression of the large outer scale had on the wall 
"burst" process. As discussed in Chapter IV, the "burst* 
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events were also (detected using a match filtering 
technique. In order to generate the velocity signature 
to be used as the impulse response filter function 
associated with the "burst" process, the VITA technique 
of Blackwelder and Kaplan (1976) was utilized. In 
keeping with the method presentee in that paper, a 

2 

threshold value of 1.2 and an averaging time of Tu^/v 
approximately equal to 10 were used. The conditionally 
averaged time series were obtained from the data by 
surrjring up the velocity values about the detection point 
and dividing that total by the number of detections. 
Therefore the conditional average of the quantity Q is 
defined as 

N 

< Q(x^,T) > y" = I i: Q(x., + T) 

3 = 1 . 

where x denotes the streamwise position at which the 
sampling was made rnd y-f denotes the vertical location of 
the detection wire. In the case of the rake, x was a 
fixed value. However, to investigate any Reynolds number 
effects on the results, the time series data taken from 
the traversable single wire probe, which was used to 
obtain one-dimensional velocity spectra, was analyzed for 
bursts. In this case, five streamwise pcsitiens at y =15 
were used for the four bcundary layer cases and tvo 

free-stream velocities. For the one-dimensicnal 
analysis, in all cases, the detection was made on the 
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signal from the hot-wire at y-?- = 15. The quantities t^ 

in the relation are those points in the time series at 

which the detection cccurrec. and the quantity T is a 

positive or negative time delay used to determine the 

temporal behavior of Q before or after detection 

occurred. Because cf the band limited nature of the 

data, events detected within less than half a wave length 

of the beginning or end of a record were not included in 

the ensemble. The visriable T. in the expression is the 

total number of detections included in the ensemble. 

The ensemble average of the fluctuating velocity, 

non-dimensionalized by the local average rms at the 11 y+ 

locations, which were conditioned on a VITA detection 

2 

with a threshold of i.2 and time averaging Tu^/v = ic, 
are presentee in Figure 90 for the regular boundary layer 

with Re, = 3819. The horizontal lines in the figure 

^ 2 

correspond to the record average m.ean value at each 
vertical height in the boundary layer. The spacing 
between these lines corresponds to 0.54 where Q is 

the stream.wise fluctuating wall shear, for y+ = 0. 

For the other non-dirrensional heights, Q is the 
strodmvise velocitVf Approxiinat6ly 100 enseintles are 

included in the average distributions, which were 
conditioned on the VITA detection in Figure 90* 

If we compare this figure to the conditioned 
streamwise velocity distributions obtained by Blackwelcer 
and Kaplan (1976) we observe the same velocity signature 
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which is associated with the "burst" process. If we 
focus on the normalized fluctuating velocity at y+ = 15, 
the detection height, we observe the characteristic local 
deceleration of the flow, very rapid acceleration and the 
subsequent slower deceleration. At this Reynolds number, 
the wavelength of the process is approximately 4c/u^ and 
the maximum peak- to-peak amplitude - is ‘approxiriately 2,2 
u(rm.s). This compares extremely well to the 
characteristics of the velocity signature at the same 
non-dimensional height in the boundary layer of 
Blackwelder and Kaplan (1976) for a somewhat lower 
momentum thickness Reynolds number of 2550. At other 
distances from the wall we observe a characteristic 
pattern which extends as far as the sensor at y+ = 500, 
which is just at the lower edge of the wake region in the 
boundary layer, Blackwelder and Kaplan (1976) 
investigated the conditioned velocities out to y+ = 100. 
In the velocity signatures at y+ values from 200 to 500 
we observe a characteristic higher velocity bulge which 
is upstream in time when conditioned on the wall events. 
This would indicate a structure which is inclined in the 
streamvise direction. This- of course, is not surprising 
since any flow structures in the boundary layers are 
subjected to a mean shear. Inclined coherent structures 
originating near the wall and extending out into the bulk 
of turbulent boundary layers was first proposed by Grant 
(1957) from space-time velocity correlation measurements. 
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At the wall/ the normalized instantaneous shear st-ess is 
significantly changed du'ting an event. The ensemfcj.e 
average of in this figure demonstrates a local 

decrease in the wall shear followed by a rapid increase 
which is associated with the rapidly accelerating _lUid 
near the wall. The time scale for the event in the shear 
distribution appears to be slightly sjnaller than that in 
the scale of the coherent velocity fluctuations. 

As was indicated in Chapter IV, reliable and 
objective "burst" frequencies cannot be obtained ty the 
VITA technique. This was consistently demonstrated by 
Offen and Kline (19731 who, in addition to showing that 
the number of detections was non-linearly related to the 
threshold value for the VITA technique, also showed that 
the other commonly used methods exhibited the same 
problem. In addition to the dependence of the detection 
frequency on the threshold, the statistical distribution 
of events is non-linearly related to that threshcj.d. 

This is demonstrated in Figure 91 in which the normalized 
probability distribution of burst events detected with 
two different thresholds are presented. As in the 
probability distributions presented in the results of 
Chapter VII for the instantaneous wall shear stress, the 
distributions in Figure 91 are plotted as number 
densities at rms increments away from their respective 
means. This type of normalizaiton should collapse 

In the distributions of Figure 91 


related distributions. 


for only a +10 percent deviation away from the 1.2 

■ • 

threshold used by Blackwelder and Kaplan (1976) , we 
observe differences in the burst statistics which are 
related to the detection threshold value. 

In spite of this criticisin of the VITA technique for 
measuring the burst frequency,’ it is valid in the sense 
that it does generate the nest coherent vel*ccity 

1 

signature that is associated with a large locally high 
rms in the streair.wise velocity near the wall. Through 
such visualization studies as that of Offen and Kline 
(1973) , who found a majority of VITA detections 
corresponded to the visual observance of a "burst" event, 
we can relate that average signature to .that for an 
average "burst" cycle. Although Offen and Kline (1973) 
found, for a low Reynolds numfc-er boundary layer that a 
VITA detection did detect burstSf they observed that the 
technique only detected 53 percent of those that were 
observed to visually occur in the region ,cf the probe. A 
similar factor of two was observed in that study for the 
other techniques that they investigated. This sane trend 
was also observed by Bogard and Tiederman (1981) • 

In light of these findings, the velocity signature 
obtained through the VITA technique was only used as a 
pattern in the match filter detection which was used to 
obtain the frequency statistics for the occurrence of 
wall events. In all the cases considered, we utilized a 
pattern that was generated from the data for that case. 
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As with the VITA detection, the patterns and match filter 
detection were generated from the velocity time series at 
y+ = 15. 

In Figures 92 to 95 we have displayed typical cases 
of the raw data for velocity time series and the 
corresponding match filter correlation for regular and 
manipulated boundary layers at each of the^ respective 
Reynolds numbers. The top time series is the 
instantaneous velocity non- dimensionalized by the average 
rms. The bottom time series in these figures is the 
cross correlation between the pattern and the velocity 
series. The dashed lines correspond to zero level in 
both af the traces in these figures. As a reference, the 
tim.e a/.is is presented non— dimensionalized with outer 
flow parameters. 

The peaks in the correlation correspond to' locations 
where a match occurred between the average burst 
signature and the time series. Inspection of points of 
positive oorrelatin peaks in these figures reveals that 
they correspond to the characteristic signature of the 
local deceleration of the flow followed by a rapid 
acceleration. Thece correlations are normalized by the 
output for a perfect correlation that corresponds to 
correlating the match filter function with itself. It is 
interesting to note that with this normalization it is 
possible to have cases where the correlation is greater 
than one. These cases are associated with detected 
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bursts that are spaced less **han an average wave length 
apart where the positive correlations can sum up. The 
negative correlations that occur in these figures 
correspond to the cases pointed out by Nishioka et al* 
(1980), namely the deceleration of the same magnitude and 
time scale as that for a burst. Without any 
specification on the sign of the velocity change, the 
VITA technique, as it was proposed originally, would 
associate these occurrences with bursts. 

The burst locations were taken as the locations of 
the significant positive correlation peaks in the match 
filter outputs. We chose a normalized correlation of 30 
percent to be the minimum for a relevant match in the 
velocity series, i.e., minimum acceptable peak amplitude. 
This value corresponded to one end of a region of 
insensitivity of the detected bursts to the correlation 
level between 0.3 and 0.8. 

Using the detection points obtained by the match 
filter, the ensemble averaged time series for the burst 
events at the 11 heights in the boundary layers were 
constructed and presented in Figures 96 - 99. The scale 
for the abscissa in these figures is the same as that for 
the ensembled time series obtained from the VITA 
detection in Figure 90. By comparing Figures 90 and 96 
ve observe that the main features of the ensem»ble 
averaged velocity signatures are basically the same. 
However, we observe that in the case of the match 
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detections, the velocity signature is coherent over a 
larger wave length than in the vita detection ensembles. 
This can be .attributed to the inclusion of non~events in 
the VITA ensembles. This leads us to conclude that the 
VITA technique does not only miss some of the wall events 
but it also includes some non-events. 

At the higher Reynolds number range in Figure 97, we 
observe a decrease in the peak— to— peak amplitude of ^he 
event signature below y+ = J.hO. Mcst of that decrease is 
a result of the much smaller level of deceleration in the 
beginning phase of the event at the higher Reynolds 
number. After the rapid acceleration in the 
instantaneous velocity, the velocity signatures exhibit a 
Diateau which was not present at the lower Reynolds 
num.ber. It is important to point out here that since the 
average signature at the higher Reynolds number exhibits 
an amplitude which is a smaller percentage of the 
turbulence fluctuating velocities, the threshold val..e in 
the VITA detection must be reduced in order to include 
all classes of the events. Inspection of the ensembles 
for the manipulated boundary layers in Figures 98 and 99 
shows them to resemble the higher Reynolds number case in 
Figure 97. One most notable difference occurs in the 
shear stress signature. In the case of the manipulated 
layers we observe a much reduced decrease in the local 
shear which is associated with the decelerating fluid 
near the wail. 
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Another noticable feature at the plateau is a 
consistent pattern of regular oscillations. This may 
very well reflect the instability of the shear layer 
developed by the "burst" cycle and observed by 
visualization, i.e., by Kline et al. (1967) and during 
the largest stages of transition, i.e., by Nishioka et 

al. (I960). 

The probability distributions of the period between 
burst events is shown in Figures 100 to 104*. A 
comparison between the probability distributions for the 
events detected by the VITA and match detection is 
presented in Figure 100, As with the previously 
presented probability distributions, these are normalized 
by subtracting the mean period and dividing by the rms. 
Comparing the statistics in this figure shows that with 
the threshold value recommended by Blackwelder and Kaplan 
(1975) , the match detection makes approximately twice as 
manv detections. This result is consistent with the 
results of Offen and Kline (1976) who found that the VITA 
technique, as well as other methods which employ a 
threshold, detected only half of the events that were 
observed to occur with visualization. Comparing the 
statistical distribution of events as reflected in the 
skewness and kurtcsis yields further results. The match 
detected events are skewed approximately twice as much as 
these of the VITA technique. This indicates that the 
VITA technique with the given threshold is not detecting 
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the events occurring with longer periods between them,. 
This would indicate that these events were less intense 
than those occurring nearer to the mean period. The 
kurtosis of the distribution for the VITA method 
indicates a distribution which is narrower than gaussian. 
In contrast, the kurtosis for the match detection is 
approximately three times greater, indicating a much 
broader distribution of events. 

Blackwelder and Kaplan (1976) , using the. VITA 
technique, estimated that the burst events occurred 
approximately 25 percent of the time in the life of the 
boundary layer. The statistics presented here for the 
match detection indicate that these events are occurring 
over a much greater percentage of time. The 
twc-dimensicnal reconstruction of the velocity from the 
hot-wire rake, to be presented later in this chapter, 
will give more insight into this aspect. Hov;ever, it is 
clear from this figure that a technique which relies on 
detection of events through locally large statistical 
quantities may be accepting only part of the total 
distribution of these events. In particular, setting a 
local threshold only includes events of a particular 
class that are sufficiently energetic or otherwise pass 
some biased criterion. Other events of the same class 
may exist but are being ignored because they do not fall 
within a narrow sampling criteria. This is not said to 
criticize this method of signal detection, but rather to 
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point out that results based on such an approach nust be 
qualified in a manner v?hich is appropriate with the 
selection parameters. 

The distributions are compared for the regular and 
manipulated boundary layers at both ranges of Reynolds 
numbers in Figures 101 and 102. These results indicate 
that the distributions become less skewed and less broad 
with increasing Reynolds number. This trend is 
consistantly demonstrated in the manipulated layers which 
appear to emulate higher Reynolds number boundary layers. 
A direct comparison of the regular and manipulated 
boundary layers at the same free-stream velocity is 
presented in Figures 103 and 104. These are the same 
distributions as in the previous two figures but 
replotted to allow a direct comparison. In all cases, 
the peak in the distribution occurs at approximately one 
rms period less than the mean. The addition of the • 
manipulator, which we observed earlier to result in lower 
friction at the wall and a thicker sublayer, leads to an 
increased absolute period, at both free-stream speeds. 

Two-dimensional Analysis of Velocity: Space/Time Data 

As described in Chapter III, space-time analysis of 
the output of the rake of 11 hot-wires was done by 
performing reconstructive two-dimensional low-pass 
filtering.' The one-dimensional filter function was shown 
in Figure 19 and the generated two-dimensional filter 
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function was shown in Figure 20. Typical records of 
two-dimensional reconstructed velocities are presented in 
Figures 105 and 106 for all four boundary layer cases. 

The velocity data are presented as contours of constant 
normalized velocity, U(y) - U(y)/u(y)^^. The dashed 
line contours correspond to "structures” having negative 
normalized velocities, or velocities less than the mean. 
The solid line contours correspond to structures 
containing velocities greater than the mean. . 

Structures corresponding to the wall bursts are 
designated by their one-dimensional signature to be the 
deceleration of the fluid near the wall followed by the 
rapid acceleration. By this definition, the burst cycle 
will correspond to a dashed line structure followed by a 
solid-lined structure for the notation of Figures 105 and 
105, In the contours of Figure 105, the acquisition was 
initiated by the detection of a burst event which is 
visible as the dashed-solid pair centered about TU /v 
approximately equal to 1 in the time sequences in that 
figure. The high velocity data were not acquired based 
on a "burst" detection and the same comparison cannot be 
made at the non-dimensionalized time of 1. Inspection of 
Figures 105 and 105 strongly suggests that the velocity 
field near the wall is in a constant state of 
acceleration and deceleration. In fact it almost looks 
like an bistable mode. This is reminiscent of the 
behavior of the wall layer proposed by Black (1966) , 
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Einstein and Li (1956) and Hanratty (1956) which modeled 
the flow to be in a constant state of dissipation and 
regeneration. Another way to interpret the 
quasi-periodic low and high velocity structures near the 
wall would be in terms of a bulk instability of the wall 
layer. This aspect will be explored more, later in this 
chapter . 

These two-dimensional velocity reconstructions also 
bring out another startling result. In the bottom part 
of Figures 105 and 106, there appears to be soce coherent 
eddy structure reminiscent of the vertex street emanating 
from the manipulator plates. The presence of these is 
clearly missing in the one-dimensional velocity spectra 
taken at this height in the boundary layer. However, the 
cv;o-di— ensional data indicate that the eddies are 
sparsely coherent, thereby resultinc in a broad banc 
distribution cf their one-cir.ensiona' spectral energy. 

The ensemble averaged reconstructed velocity 
profiles conditioned at the time-instant when the digital 
images were acquired are shown in Figures 1C7 and 1C8 for 
the regular and manipulated boundary layers, 
respectively. These indicate that most cf the digital 
images were triggered on a burst event, although the 
detection is observed to be more sensitive to the high 
velocity sweep of fluid near the wall. This information 
will be used in the interpretation of the enserbled 
digital images. 
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The ensemble averaged reconstructed velocity 
profiles conditioned at the match filter detected time cf 
an event are presented in Figures 109 - 112 for the 
regular and manipulated boundary layers at the two 
free-stream velocities. The range of contour levels in 
these figures is the same as that for the previous tvc 
figures. The features depicted by the contours about the 
detection point, TUa=/6 = 0, correspond to the average 
"burst" structure. If we compare the average structure 
between Figures 109 and IID we observe a streamwise 
extension with increasing Reynolds number. In the 
manipulated layers in Figures 111 and 112 we observe that 
the structure has a longer streamwise extent compared to 
the regular boundary layer. This is consistent with the 
idea of the manipulated layers emulating higher Reynolds 
number boundary layers. The fact that some coherence is 
evident away from the detection point is indicating cf 
the quasi-per iodic nature of the process. 

In order to better understand the mechanisms of the 
wall layer events, the average structure associated with 
burst periods greater than 3c away from the mean period 
was investigated. This was done so as to pinpoint what 
differences, if any, existed in these structures for 
which the nearest event occurred 3c time scales before. 
These are presented in Figures 113 and 114 for the 
regular boundary layer at the low and high 
free-stream speeds, respectively. The non-matching cf 
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some contours is a result of having to form the enseables 
in separate 64 point segments so as ^o maximize the 
number of realizations in the ensemble average for these 
doubly restrictive conditions in the band-limited time 
sequences. The 3c burst time scale corresponds to 
approximately one-third of the continuous sampling time 
and as seen in Figure 99 only approximately cne percent 
of the total events occur at periods greater than 3c away ' 
from the mean. As a result of these limitaticns/ the 
ensembles in Figures 113 and 114 contain only 
approxim.ately six realizations. However, certain 
conclusions can be drawn from the results. In particular 
these figures indicate that when an event has not 
happened for a previously "long" time, the wall structure 
remains at a low velocity state, compared to the mean 
flew, which is observed to extend vrell up into the 
constant stress layer. A high velocity shear layer is 
observed to form above the low vclccity structure prior 
to the fermation of the high velocity vail structure. 

This is evident in the ensemble averages shown in beth of 
the figures. 

In order to lock at the wave-numiber distribution of 
the structures near the wall, two-dimensicnal spectra 
were computed of the two-dimensional reconstructed 
velocity distributions of the type shewn in Figures 105 
and 1C6 . The two-dimensional spectra are presented in 
Figures 115 - 118 for the regular and manipulated 
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boundary layers at the two free-stream speeds. The 
spectra are nornalized by the volume under their 
distributions which corresponds to the two~dimens\onal 
variance and plotted as contours of constant spectral 
levels. In addition, the two-dimensional spectra are 
plotted as all four quadrants of positive and negative 
wave numbers. 

At the low free-strear> speed for the regular 
boundary layer in Figure 115, we observe narrow band 
peaks at equal to approximately 2 and 0.5, along a 

spatial frequency 5cy/6 equal to approximately 0.005. The 
higher frequency corresponds to c non-dimensional period 
of 0.5, which is evident as the smaller scales observed 
at the lower Reynolds number in the top part of Figure 
105. The other peak corresponds to an approximate 
non-dimensional period of 2 Uc/d , which is the 
approximate mioan burst period determined from the 
one-dimensicnal match filter analysis. The peak in the 
spatial frequency coordinate corresponds to a vertical 
scale of approximately 0.2f which is the approximate 
scale of the structures with velocities greater than the 
mean in Figure 105. The other spatial peak corresponds 
to a vertical scale of 0.C5J which is the approximate 
scale of the wall structures having velocities less than 
the m.ean. 

At the higher Reynolds numiber, the two-dimensional 
spectrum in Figure 116 does not shew a peak at fi/U^ = 2 
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and the other peaks are observed at lower wavenumbers. 

The peak at approximately - 0.25 corresponds to a 

structure period of approximately /).02S seconds which is 
approximately the burst period for the regular boundary 
layer at the higher free-straam speed. These structures 
have a spatial wave number of approximately ky/c=0.004 
which corresponds to vertical scales of approximately 
0.165. Thus the wall structures are observed to increase 
their stream.wise extent and decrease their vertical extent 
with increasing TReynolds number. 

The two-dimensional spectral distribut icns for the 
manipulated boundary layers exhibit the same trends as 
brought out by the match filter analysis, in that v:e see 
an increase in the streamwise extent of the wall structures 
with the addition of the manipulator- In addition, the 
spectra also indicate that the vertical scale of the wall 
structures is increasing. This is consistent with the 
observed growth of the sublayer in tr.e mr.nipilated boundary 
layers . 

T^wc- dimensional Image Analysis 

The rake of hot-w’ires provided information about the 
spatial features of the boundary layer in a region 
primarily below the mid-point of the boundary layer 
thickness. In order to view the coherent structures 
across the whole bulk of the boundary layer, the 
dicitized images cf the smoke-’wire visualized flowfield 
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were analyzed. Grey-level representations of typical 
digitized images are shown in the top half of Figures 11? 
and 120 for the regular and manipulated boundary layers. 
These are photographs taken from the color television 
monitor as output from the KDP display processor with a 
black and white color table which gives an increasing 
grey intensity corresponding to the digitized luminance 
levels in the smoke image. In these representations, 
black corresponds to the lowest intensity and white 
corresponds to the maximum intensity value. The images 
in these figures correspond to the zoomed-in view of the 
flow in the near vicinity (a^proxirately 2j) of the 
hot-wire rake, which is visible in the left bottom corner 
of the photograph. 

In order to distinguish the coherent features in the 
visualized flow, the images were match filtered utilizing 
the same two-dimensional filter used in the velocity 
reconstructions presented in the previous figures. This 
was appropriate since the physical scalings were 
identical in both of the processes. The details of the 
match filtering for the digitized images were presented 
in Chapter V. The results of that process yields 
enhanced images consisting of features whose spatial 
variations in luminance correlated well with the spatial 
characteristics of the two-dimensional filter response. 
We therefore are searching for particular scales of 
features in these images which are consistent with those 
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brought out in the velocity distributions of the 
two-dimensional reconstructions. 

The bottom photographs in Figure 119 and 120 are the 
corresponding match filter reconstructions of the 
respective digitized images in those figures. In these 
cases, the black/white color table was loaded with a step 
function sc as to distinguish distinct levels and form 
contours of constant match correlation. Black 
corresponds to a negative correlation and white 
corresponds to a maximum positive match. Comparing the 
digital images and their corresponding match-filter 
representations points out the features that this 
processing highlights. The hot-wire rake, which yields a 
strong correlation, was left in the processed images to 
aid as a reference to the photographs. 

Comparison of the images for the regular and 
manipulated layers in Figures 119 and 120 reveals some 
differences. It is important to remember that the 
acquisition of these images was initiated by the 
detection of a "burst" event. This was indicated by a 
large local variance in the velocity signal at y =7 in 
the boundary layer. At that instant :n the regular 
boundary layer, the match-filter enhancement in Figure 
119 reveals some coherent features near the wall and in 
the outer part of the flow that extend toward the wall. 
In the manipulated layer in Figure 120, we observe a 
larger feature near the wail and what appears to be a 
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regular “ddy pattern in the outer part of the flow. 

These, hcwever, correspond to one realization only of the 
"liiursf event. In Figure 121, the ensanble-averaged 
enhanced-iraages of these events are presented. As v;as 
indicated earlier in this chapter, the detection of wall 
events based on a locally high 'variance is not always 
reliable. Therefore, the image realizations Included in 
the ensemble were those which corresponded to a 
match-filter detection at that instant in the velocity 
signal, which was acquired simultaneously. The 
acceptance rate in this case was approximately one out of 
th'ee . 

The ensemble average of the images without any 
enhancement tended to smear the details, but 
conditionally averaged match filtered images which 
are presented in Figure 121 revealed seme features. 

If we focus cn the result for the regular boundary 
layer, in the tep of the figure, we observe what 
looks to be a relatively- large slender feature which is 
inclined in the flew direction. This shape is 
reminiscent of the Townsend (1970) inclined eddies or the 
Brown and Thomas (1977) description of the outer inclined 
structures they associated with the wall ■bursting". The 
loc^.l maximum at the outer tip of the feature indicates 
that the smoke was rolled-up in that region suggesting a 
possible eddy structure. This may very well be a visual 
ensemble of the hairpin-like vortices. 
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In the ensemble averaged j.mages for the manipulated 
boundary layer in the bottom half of Figure 121, we see a 
similar feature. However we find it residing in a region 
below approximately 60 percent of the boundary layer 
thickness. Above that we observe the alternating light 
and dark regions that are indicative of the ed'^y pattern 
in a vortex street, most likely emanating frcm the 
manipulator plates. If we compare these images to the 
ensem.ble-averaged reconstructed velocity distributions in 
Figures 107 and 108 we see some similarities. The 
velocity distributions were normalized so as to bring out 
features that stood out against the oackground 
non-coherent turbulence. For the regular boundary layer 
the normalized velocity contours in Figure 107 reveal a 
lifted structure which is inclined at an angle similar to 
that of the feature in the image. In the manipulated 
layer, the structure implied from the velocity contours 
in Figure 108 is observed to exhibit less of a vertical 
extent than in the regular layer. This was also evident 
in the ensemble image for the manipulated layers. 

In order to view a greater streamwise portion of the 
boundary layer, the images viewing approximately 10c in 
the flow were analyzed. These are shown in Figures 122 
and 123 for the regular and manipulated boundary layers. 
Again the match filter enhanced images were generated for 
these realizations which also were acquired upon the 
detection of a wall event. These are comparable. 
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although not the exact same realizations, to the 
two-dimensianal reconstructed velocity tin^e series in 
Figure 105. There are clear and identifiable differences 
in the enhanced irnages for the two cases. The regular 
boundary layer exhibits long coherent features that 
extend across the bulk of the boundary layer and which 
are inclined along the mean shear. In the manipulated 
boundary layer these features are not evident. Rather, 
the outer part of the boundary layer flow appears to be 
dominated by a stream of alternating vortex patterns. 

This feature was pointed out in the discussion of the 
features for the two-dimensional reconstructions of 
earlier in this chapter. 

Considerably more analysis is needed to interpret 
these images. However, they do bring out coherent 
features which are linked with the burst events and which 
are consistent with the picture we obtained from the 
simultaneous velocity surveys. Wo hope that they will 
also provide us with the time sequence of events from an 
objective and quantitative assessment of flow 
visualization records . 
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CHAPTER X 
DISCUSSION 

In the previous two chapters , the consequences of 
suppressing the large-scale intermittent outer structures 
in turbulent boundary layers were presented.^ In this 
chapter these results will be tied together so as to 
pinpoint the physical mechanisms involved in these 
manipulated boundary layers. In doing this, emphasis 
will be placed on the role of the intermittent outer 
corrugations on the boundary layer growth as well as on 
the production and distribution of the turbulence scales 
of motion. Next, the cuter flow processes in the boundary 
layer will be tied to the mechanism for turbulence 
generation in the wall layer through their interaction 
with the coherent motions near the wall. This analysis 
will attempt to lead to a unified phenomenological 
r eor esen t a t i on for the mutual interaction of inner and 
outer fluid layers in boundary layer flows. Whenever 
possible, the results of other investigators will be 
drawn upon for comparison and to reinforce the 
conclusions arriv'ed at in this investigation. 

Meehan isrr.s G overn i n g Cpe ration 
of Parallel-Plate Manipulators 

The manipulated boundary layers were observed to 
have marked differences in their mean and unsteady 
characteristics when compared to the regular turbulent 
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boundary layers over a comparable Reynolds number range. 
The major differences occurring in the manipulated 
boundary layers are characterized by the following; 

1) a visual evidence of the suppression of large 
outer turbulent scales that leads to a lack 
of intermittent e^ccursions of high speed 
potential fluid into the boundary layer, 

?.) a decrease in the streamwise growth of the 
boundary layer thickness and integral 
thicknesses 6^, 6^, and 6^, 

3) a disproportionately larger shape factor ^-^ 2 ' 

4) a 30 percent decrease in the local skin-fricti.on 
coefficient over the regular boundary layer at 
comparable Reynolds numbers, 

5 ) a growth in the sublayer by approximately 17 
percent which is refleced in an upward 
translation of the log-linear velocity 
distribution , 


6) a 17 percent decrease in the von Earroan 
p constant reflected in a more sloped log-linear 

velocity distribution for the manipulated 
layers, 


a decrease in the entrainment velocity which is 

} 
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also observed to decrease with increasing 
Reynolds numbers, 

8) a reduction in the wake strength coefficient, 
n, in the manipulated boioncary layers at 
comparable Reynolds numbers, 

9) a lower level of u'/Uno in the manipulated 

boundary layers below y. = 100 and an increase 
in u'/u^ in the wall layer up through the 
constant stress layer, ^ 

10) growth of the constant stress layer and an 
increase in the rate of growth of the 
streamwise length scale of turbulent eddies, 
and 

11) a reduction in the absolute frequency of 
occurrence of the so-called "burst" events in 
the manipulated boundary layers. 

We look to these documented changes in the 
manipulated boundary layers in order to ascertain the 
mechanisms governing the operation of the parallel plate 
manipulators in a turbulent boundary layer. It should be 
kept in mind that these cited effects appear to persist 
for appreciable distances (at least 70 boundary layer 
thicknesses) downstream of the manipulator. 

If we consider the results of the parametric study of 
this class of devices in Chapter VI, we note that an 
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optimum exists for the different parameters in the 
different manipulator plate configurations. The table in 
Figure 22, for which the optimum parameters are scaled 
with the characteristic boundary layer lengths, indicates 
that the maximum reduction in the outer intermittent 
scales is produced by devices which have dimensions and 
are placed at heights that are comparable to the average 
si 2 e of the eddies which are to be suppressed. The fact 
that two-plate parallel and series manipulators worked as 
well as the four plate manipulator suggests that the 
lower plates in H-1 had little effect on the large scale 
turbulence. What is extremely interesting in the 
parametric study is that the optimum spacing between the 
streamwise oriented plates of M-4 is 6 to 8 boundary 
layer thicknesses which correspond to the average life of 
a large scale eddy estimated from the velocity 
correlation measurements of Blackwelder and Kovasznay 
(1972). 

There appears to be four key mechanisms for the 
observed suppression of the outer scales in the turbulent 
boundary layer. Any one or all of these mechanisms may 
be acting at once although seme are raost likely more 
dominant than others. The first explanation for the 
effectiveness' of these simple devices is that the 
manipulator plates restrict the vertical velocity 
components associated with the large scale motions 
(Kovasznay, 197 in very much the same way as a 


177 



ORIGT'JAL TACF. iS 
OF POOR QUAUTY 

honeycomb supresses the lateral velocity components of 
grid turbulence within eight mesh lengths (Loehrke and 
Nagib, 1976) . The fact that these scales dissipate so 
rapidly may be an indication that they are relatively 
inactive motions. 

The unsteady circulation generated by the flat plate 
in the fluctuating turbulent flow leads to the second 
mechanism. A flat plate at zero incidence in an unsteady 
shear flow similar to the boundary layer will generate 
circulation in the same and opposite sense to that of the 
large scale structures. The exercise carried out with a 
one-plate manipulator at various positive and negative 
angles of attack in Chapter VI was aimed at demonstrating 
this mechanism. In those experiments we observed both 
suppression as well as augmentation of the large scales 
depending on the angle of attack being negative (or zero) 
or positive , respectively . The augmentation of the large 
scale motions is an appealing idea for improved heat 
transfer applications. 

The third mechanism is based on generation of nev: 
scales in the outer part of the turbulent boundary layer 
as a result of the vortices being formed in the wakes of 

the plates. Although the one-dimensional velocity 

« 

spectra taken at the height of the manipulator plates at 
downstream stations indicate that these vortices rapidly 
lose their identity through the turbulence velocity 
fluctuations, the two-dimensional velocity 
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reconstructions and imago processing records indicate 
that some coherent remnants may exist in intermittent 
packs far downstream of the device. These eddies, which 
are of a smaller scale than the larger outer bulges in 
the boundary layer may act to redistribute small scale 
turbulence. Through visualization of the outer flow of a 
turbulent boundary ]ayer, Taneda (1981) and Taneda et al. 
(1981) have demonstrated that a roost amplified mode may 
exist for the passage of large outer flow structures. 

The generation of a vortex street with passage 
frequencies away from the most amplified case may disrupt 
this process. 

The last mechanism, which may have a role in 
suppressing the large outer scales and in redistributing 
the fine scale turbulence through a greater part of the 
boundary layer, results from having a wake profile 
embedded in the mild mean velocity shear in the outer 
part of the boundary layer. The turbulent energy 
equation for a two-dimensional incompressible mean flow, 
outside the viscous sublayer (Townsend, 1956) takes the 
following form 


1 

P 


^-2 3-2 

(U + V ^ ) 


9x 


3y 


advection 


- puv § pq^v) + pe = o 

production diffusion dissipation 


Normally in the outer wake part of the boundary layer ^ 
is small, as well as the uv correlation, thereby 
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r 6 sultin 9 in a near zero turbulence production. In the 
wake of the manipulator plates, we have a locally large 
dD/dy. This is evident just downstream of the device as 
shown in Figure 52. On the wall side of the manipulator 
plate we might expect that dO/dy would be negative but 
because of the mean shear the U{y) distribution directly 
below the plate is nearly constant as shown in Figure 52. 
On the potential flow side of the manipulator plate, we 
observe a substantially larger positive velocity 
gradient. As a result of the unsteady vorticity 
generated by the manipulator plates, we also might expect 
that the uv correlation would increase. Therefore, the 
vake of the plates is a likely site for turbulence 
production in the manipulated boundary layers. Additional 
evidence and discussion of this is presented in Appendix 
B. 

The Ir.st mechanism may explain why the manipulated 
boundary layers are filled with fine-scale turbulence as 
decicted in the photographs in Figure 51. Furthermore 
this would explain why the u'/u^ distributions in the 
manipulated layers in Figures 54 and 55 have higher 
values than their regular boundary layer counterparts. 

In the regular case, the streamwise velocity stress was 
found to be related at all heights in the boundary layer 
to the wall shear velocity, u^ , as was shown in Figures 
31 and 32. This was not the case in the manipulated 
boundary layers. However, since u. is a measure of the 
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turbulence production at the wall» it would not account 
for other sources which contribute to the turbulent 
stresses away froF. the wall* With the added turbulence 
production in the wake of the manipulator plater we would 
expect a thickening of the constant stress layer which 
indeed was observed in the one-dimensional velocity 
spectra. 

Although this concept is very appealing for 
explaining one of the turbulence control mechanisniSr the 
mean velocity distributions indicate a fairly rapid decay 
of the wake of the plates. Therefore, the turbulence 
production must be localized to the near wake of these 
dates. However, once the fine scale turbulence is 
generated, it might be convected for long streamwise 
distances • 

Each of these mechanisms that may govern the 
effectiveness of these large scale turbulence 
manipulators in the' boundary layer have justifying 
ascects. Most likely, all of them, are act.ng to 
suppress the large-scale outer structures and to 
redistribute turbulence throughout the boundary layer. 
The fact that an optinu-- exist? for this process 
indicates that these devices are acting on some 
fundamental properties of the turbulence in the boundary 

layer. 
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Drag Reduction and Reynolds Number 
Effects on Bou ndary Layer 
Large Scales 

In Chapter IX we documented a 30 percent reduction 
in the local skin friction in the manipulated boundary 
layers. Accounting for the drag of the device, this 
resulted a net reduction in the friction coefficient of 
21 percent by the last streamwise station at the end of 
the tunnel test section, i,e., about 55 boundary layer 
thicknesses downstream of the manipulator. This part of 
the chapter will address the physics behind this behavior 
and make some judgments as to how long this effect may 
persist. 

In the photographs of Figure 51, along with those in 
Figures 24 and 25 from the parametric study, we observed 
a lack of the intermittent excursions of high speed 
potential fluid rnto the recesses in the boundary layer 
wita '-.he manipulator. As a result of this we measured a 
decriase :.n the strearwise growth of the boundary layer 
and all of its integral thicknesses as evidenced by 
Figures 56 and 57. By impeding the high speed fluid from 
reaching close to the wall we would expect to reduce the 
high local instantaneous wall shear stresses and thereby 
the time averaged value. The reduction in the growth of 
the boundary layer and reduction in the local average 
wall shear stress are a result of the lower rates of 
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entrained potential fluid into the boundary layer. This 
vas demonstrated in Figure 69 which presented the mean • 
rate of propagation of turbulent fluid into the 
free-stream or the entrainment velocity, as a 
function of Reynolds number in the regular and 
manipulated boundary layers. In this figure we observe 
an approximately 40 percent reduction in the entrainment 
velocity with the large scale manipulators. Bradshaw 
(1967) found that the quantity Vp/u® was very nearly 
proportional to i/pU^ in equilibrium and 
non-equilibrium boundary layers. Furthermore that ratio 
of the large-eddy velocity fluctuations to the 
small-scale velocity fluctuations in the central part of 

2 k 

the boundary layer is proportional to ('/pU ) . 

Moreover, the ratio of advecticn or diffusion to 
production in the outer part of the layer is proportional 
to C^. 

Figure 69 demonstrates that the Reynolds number 
distribution of friction coefficient is related, within a 
proportionality constant, to the ratio of the entrainment 
velocity to free-stream velocity. Therefore, removing 
the large outer scales leads to a reduction in the 
entrainment rate in the manipalateld boundary layers 
which then results in a subsequent reduction in the local 
friction coefficient. 

Considering the regular boundary layer case in this 
figure we find that the entrainment velocity is 
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decreasing with increasing Reynolds nuraber in a manner 

proportional to C^(Rej^). This indicates that the large 

outer scales, defined by the regions of interinittency in 

the turbulent boundary layer, are becoming less coherent. 

In a similar fashion the ratio of diffusion to production 

in the outer part of the layer is decreasing like 2C£(x). 

This is a result o- the decrease in the ratic of the 

large-eddy velocity fluctuations to the small-scale 

velocity fluctuations proportional to {2C^) . Thus the 

cuter part of the bound: y layer is filling with smaller 

_ '/> 

scale turbulence like (He. )~ . The process has been 
observed through visualization records of turbulent 
boundary layers over a range of Reynolds numbers 
including transitional values by Nagib et al. (1978). 

In the manipulated boundary layers we are hastening 
this process by removing the intermittent outer scales 
and thereby "aging" the boundary layer in a Reynolds 
number sense. These manipulated boundary layers 
therefore emulate higher Reynolds number flows. By this 
same argument, however, '-.'e are implying that the large 
outer scales are associated with "young" boundary layers. 
Extending this argument one step further leads to the 
hypothesis that these large outer scales are the remnants 
of the turbulence transition process. This reasoning is 
reinforced by the fact that we have not observed a 
recurrence of these large outer bulges in the manipulated 
boundary layers within the 70 bounday layer thickness 
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available to us downstream of the device. 

Given the relationship between the entrainment and 
friction coefficients, we examined the Reynolds number 
dependence of for the regular and manipulated boundary 
layers in Figure 124, This figure demonstrates that a 
power law relation exists between the friction 
coefficient and Reynolds number. In particular, the 
values for the regular boundary layer are observed to 
decrease like the -1/4 power of This is predicted 

by the Ludwig-Til Imann relation for a constant shape 
factor n.anipulated boundary layers, the 

dependence on Reynolds number fails on a higher slope of 
-3/7. The fact that the C^(Re^ ) distribution is not only 
shifted but also diverging from the regular case 
indicates that the boundary layer is not only "aged" by 
the manipulator, but that it is "aging" faster with 
increasing Reynolds numbers. 

We can estimate the equivalent Reynolds number 
corresponding to the manipulated boundary layers by 
translating the distribution in Figure 124 so that the 
values fall on the -1/4 distribu-tion of the regular 
layer. This is equivalent to multiplying those values by 
a coefficient which corresponds to a translation on the 
log axis. The coefficient value, however, increases with 
increasing Reynolds numbers in the manipulated boundary 
layer to account for the greater slope. This process 
yields a range of equivalent Reynolds numbers for the 
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manipulated boundary layers from 7,000 to 22,000. 

Based on the above arguments and, in particular, the 
hypothesis that the large scales are remnants of the 
transition process, one may expect that the effects of 
the manipulator on the wall friction would persist for 
very long distances. Only if the large scales would be 
regenerated by a truly turbulent boundary layer, one may 
expect a return toward the regular boundary layer 
correlation. Additional evidence on this matter is 
presented in the next section. However, if the return of 
the l^rge scales is not likely or is very slow, the 
manipulators investigated here may present a very 
attractive tool for reducing drag over bodies with 
turbulent boundary layers. 

Wall Layer Instability 

The one- and two-dimensional analysis of the 
coherent features in the wall region in turbulent 
boundary layers indicated that the so-called wall 
"bursts" which were determined by other investigators to 
occur approximately 25 percent of the time (Blackwelder 
and Kaplan, 1976) actually occur in a continuous 
quasi-periodic process. The random nature associated 
with these events is a result of their extremely broad 
distribution of wavelengths. It is this behavior which 
has prompted us to associate these events with a bulk 
instability mechanism in the wall layer. In thir section 
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of the chapter we will address this issue and relate the 
frequency of occurrence of these events to the boundary 
layer parameters so as to identify the physical process 
governing these important features in the wall layer. 

In the results in Chapter IX, we demonstrated that 
the one-dimensional match filter detection of the wall 
events gave the same results as those obtained from 
two-dinensional spectral analysis of the reconstructed 
velocity distributions. Given the consistency of both 
approaches, the results from the traversing single wire 
probe and stationary rake of wires were used to establish 
the Reynolds number dependence for the frequency of the 
wall events. In line with an instability approach we 
formed a Strouhal number based on a wall layer velocity, 
U-, ar.c. length, , to non-dimensionalize the wall event 
frequency. That Strouhal number is 



L- , 


This can be put in terms of quantities characteristic of 
the sublayer, where we take an arbitrary value for y+ 
which is near the edge of the sublayer to obtain 

St^ = fy^v/u^ = f 
y u_ 

This yields the traditional inner variable non- 

d iper.sional izat icn of the burst period or 
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The quantity 1/st^ is plotted versus Reynolds number 
in Figure 125 for the regular and manipulated boundary 
layers. In this form the event frequencies show some 
very interesting trends. If we concentrate on the 
frequencies obtained from the rake of wires, shown as 
half filled symbols, for which the change in 
obtained by changing we find the trend denoted by the 

dashed line in the figure. The slope of this line, both 
for the regular and manipulated boundary layers, is 1.5. 
If, however, we measure the burst frequency at one 
free-stream velocity and increase , by increasing the 
downstream x location, we obtain a different trend as 
denoted by the solid lines in Figure 125. In the regular 
boundary layer for both free-stream velocities, the slope 
of the line is -3/8. The manipulated boundary layers 
’show a slightly steeper slope, however. This is a result 
of the differences in the distributions, which were 
shown in Figure 124, and will be accounted for in the 
next figure. It is clear that both trends must intersect 
at the X station for the hot-wire rake. The two data 
values at that station are an indication of the 
repeatability, albeit from two different probes and two 
different experimental runs. 

The results in Figure 125 clearly indicate that in 
order to collapse the nondimensionalized burst period we 
need to introduce two parameters, namely, and 

Without introducing this second parameter we cannot 
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collapse the data which occur over an overlapping range 
of momentum thickness Reynolds numbers. Given this 
insight as to this two parameter system, we can reexamine 
the data taken by other investigators and look for 
similar trends. 

Figure 126 presents the data from this experiment 
along with the "burst" frequency data of other 
investigators who employed a number of detection 
techniques. The solid symbols, corresponding to the 
manipulated boundary layers, have been shifted out to 
their equivalent higher Reynolds number values in this 
figure (see last section for more details) . The scatter 
of values in this figure is typical of what has been 
observed in the past. However, if we focus on the data 
taken with flow visualization by Kim et al. (1971) 
corresponding to the flagged circles and by Schraub and 
Kline (l?6o) corresponding to the flagged squares, in the 
lower left corner of the figure, we observe consistent 
trends. The cited data of Kim et al. (1971)was taken 
with increasing free-stream velocities at a constant x 
station and it demonstrates the expected upward trend. 

The data of Schraub and Kline (1965) was obtained at a 
constant free-stream speed and increasing x locations and 
it demonstrates the expected downward trend. The consist- 
ent trend is observed in the visually detected burst 
frequency results of Corino and Brodkey (1969) in a 
turbulent pipe flow. There they obtained an approximate 
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1.5 power dependence on Reynolds number with increasing 
velocity. They, however, could not explain the opposite 
trend exhibited by Schraub and Kline (1965) . It is clear 
from this that early investigators were compounding the 
probleir. by indiscriminately changing both of these 
parameters so as to obtain a range in one of them. 

Based on these results, we used our data to obtain 
the dependence on Re^^ required to collapse the curves for 
the regular boundary layer. That dependence is plotted 
in the log-log distribution in Figure 127 versus Re^^ , 

In collapsing the data for the regular boundary layer we 
have also succeeded in collapsing all of the data values 
in the previous figure on the sam.e line with a slope of 
5. From this we conclude that 

+ X ^ 

St (x ) = CONSTANT 

The implications of this result are that in a 
boundary layer flow which is self— similar by all 
para-etric definitions, the wall layer bulk instability is 

dependent not only on the characteristics of the layer (<5;) 
but also on how you reached these characteristics (x) . 

The relatively large power (5) on these parameters is a 
result of the near parallel nature of their growth in the 
boundary layer. In that regard it may be thought of as a 
non-parallel flew correction to the instability process. 
The fact that this correlates the results taken over a 
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range of Reynolds numbers where the boundary layer is 
fu^y turbulent and dovn to = 600, where the flow is 

transitional, indicates that the mechanisms are the same 

in both of these flow regimes. 

It is only simple to conclude at this point that the 
large scale structures, which have been studied in great 
detail over the last 20 to 30 years *n laboratory 
boundary layers, are remnants of the transition process. 

The turbulent spots observed and documented again in 
detail during the transition of boundary layer and 
channel flows must persist well into the ti:rbuler.t 
regime. The hairpin-like vortices observed within them 
must be the same ones that have been captured recently by 
various techniques in laboratory turbulent boundary 
layers. However, there is no strong reason to expect 
these spots to oe regenerated in a turbrlent boundary 
layer if they are suppressed by a manipulator. 

In presence of these spots or large scales, the wall 
layer reacts to them by undergoing some violent wall 
events or "bursts". However, this reaction to the 
■external" forcing is only part of the sensitivity of the 
viscous sublayer. The instability of the same viscous 
sublayer appears to be a continually ongoing process 
which scales with the inner variables and is dependent on 
as well as Re - . Re^ dependence is not surprisina 

' X : * 

when one considers tne instability of this layer 
dominated by viscous effects from the wall. 
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The aanipuJatrrs may then be viewed as a way of 
"aging" the boundary layer past its post-transition 
state. This of course leads to the attractive lower 
values of high Reynolds number boundary layers. If the 
wall layer instability would lead to ?arge scales like 
those of the "spots"^ the return to the original bebavicr 
nay be forthcoming. Even so, this return may be very 
slow and net gain for technological devices may be 
feasible. If the wall generated scales would not lead to 
such large scales, better performance may be even 
available to us. 
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CONCLUSION 

The conclusions address two separate but 
interrelated subjects that were experimentally 
investigated. First, the main results dealing with the 
management of turbulent boundary layers with the 
objective of reducing their drag are summarized. In the 
second part, the main conclusions describing the n:^w view 
of turbulent boundary layers which emerged from the 
detailed measurements are put forth. An extensive 
discussion of the findings of the experiments is 
available to the reader in Chapter X and the details of 
facilities, techniques, flowfields and results are 
contained in earlier chapters. 

The following section will address the various 
conclusions in detail. However, the major findings can 
be summarized by the three statements; (i) Turbulent 
boundary layers can be manipulated by suppressing their 
large scales and net drag reduction can be achieved, (ii) 
These large scales, which have been the subject of many 
recent laboratory investigations, turn out to be remnants 
of the transition process, i.e., slowly decaying 
turbulent spots that are imbedded in the boundary layer, 
and (iii) The wall layer events leading to the production 
near the surface, which have been termed "bursts", appear 
to be a result of an instability of the viscous sublayer. 
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It was then not surprising to find that they are 
correlated by the two combined parameters Re^, and Re^ , 
rather than by the first one alone, as has been assumed 
for many years. 

Conclusions 

The experiment verified that large-scale outer 
structufes' defined by the intermittent excursions of 
high-speed potential fluid into the boundary layer can be 
suporessed by simple parallel— plate "manipulators . A 
parametric investigation of the different types of 
parallel plate devices verified that optimum designs 
existed for the different configurations of stacked and 
tandem plate arrangements. A final design consisting of 
twc single plates residing at a height of 0.85 and spaced 
8; apart in the streamwise direction was chosen for the 
final set of measurements. This manipulator was selected 
because it provided the minimum device drag among the 
arrangements considered and because the plates were well 
away from the wall layer. The latter criterion was 
important so that the effect of suppressing the large 
outer scales on the wall layer events can be examined 
without the added influence of the plate wakes. Several 
mechanisms have been identified to explain the operation 
of these manipulators and are discussed in Chapter X. 

As a result of suppressing the outer scales in the 
turbulent boundary layer with the manipulator, over the 
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range of frcs 2000 to 5800, we obtained a decrease 

in the streamwise growth of the layer and its integral 
thicknesses. The entrainment velocity was decreased by 
approxiinatlely 40 percent in the manipulated boundary 
layers. Consistent with the reduction in the entrainment 
was 3 30 percent decrease in the local friction 
coefficient.. After accounting for the drag on the 
manipulator plates, the net drag reduction in the 
manipulated boundary layers reached a value of 
approximately 20 percent within 55 boundary layer 
thicknesses downstream of the device. Further net drag 
reduction is likely at farther downstream distances since 
we see no evidence for recurrence of the large scale at 
the farthest downstream distance in the test section.- 

In the absence of the intermittent outer scales, the 
manipulated layers exhibit a 17 percent thickening of the 
wall sublayer. Consistent with the this was a vertical 
translation of the log-linear distribution which has also 
been observed in boundary layers with drag reducing 
additives. In addition, the frequency of bursts at the 
same x location is reduced by a comparable amount. The 
manipulated layers consistently exhibit a 19% lower u co 
turbulence intensity in the region below y+ = 100 
although they display a higher u'/u ^ distribution over a 
bulk of the layer, primarily because of their low skin 
friction. 

Various features of the manipulated layers are 
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consistent with the idea that they emulate higher 
Reynolds number boundary layets. The wall friction 

for the manipulated layers is equivalent to 
that for regular boundary layers in the Reynolds number 
range from 7000 to 22000 based on the momentum thickness. 
This concept suggests that by removing the large-scale 
turbulence, we hasten the "aging" of the flow in a 
Reynolds number sense and that at higher Reynolds 
numbers, the large intermittent outer scales are 
diminished. 

With the aid of the manipulators we were able to 
isolate and focus on boundary layers with comparable 
oUi_er characteristics and a lower wall fricticn, i.e., 
smaller wall variables. This led to a better insight 
iSwO the mechanisms of the wall layer. Introducing new 
techniques for objective detection and description of the 
two-dimensional velocity field near the wall, these 
mechanisms were documented. These techniques included 
match filtering, two-dimensional velocity reconstruction 
as well as digital image processing and pattern . 
recognition. 

Correlation of the results established that the 
frequency of occurrence of coherent wall structures is 
simultaneously dependent on two parameters, namely on 
Re., and Re^ . By forming a Strouhal number, St"^ , based 
on a sublayer thickness and velocity we obtained the 
familiar inner scaling nondimensionalization of the event 
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frequency. A correlation between this Strouhal number 
and the two parameters was formed and used to collapse 
the results of this investigation for both regular and 
manipulated boundary layers. For the first time, results 
of other investigators, over a range of momentum 
thickness Reynolds numbers from 600 to 38000, are 
correlated by the resulting relation, i.e.. 



The ratio multiplying the Strouhal number can be thought 
of as a non— parallel flow correction for the wall layer 
bulk instability frequency, it further points out that 
the mechanism governing the coherent motions near the 
wall must include information about the history or 
initial conditions of the boundary layer. The fact that 
this relation correlates the data over the Reynolds 
number range for fully turbulent flow, and down to 
transition Reynolds numbers below 1000, confirms that the 
physical processes near the vail are similar for both of 
these flow regimes. 

Based on the above conclusions, there is little 
reason to expect these spot-like large scale structures 
to be regenerated from a turbulent boundary layer after 
they are suppressed by a manipulator. As pointed out 
earlier during the discussion, one may view the 
manipulator as means of "aging* the boundary layer past 
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its post-transition state. If the wall-layer 
instabilities would lead to large scales like these 
"spots", a return to the original trends, including the 
skin friction coefficient, may be forthcoming. Even so, 
this return is likely to be very slow and net gain for 
technclogical devices may be within our reach. Better 
performance may materialize if wall generated scales in 
an "aged" turbulent boundary layer would not lead to «5uch 
large structures. While in free shear layers, wakes and 
jets, the mean velocity profiles may provide continuously 
unstable modes, the outer boundary layer velocity profile 
is most likely not to do the same. It is conjecture', 
that the main source of instability leading to 
exponential growth and rapid production of turbulence 
would continue to be the wall layer, i.e., smaller 
scales. • 

Recommendations 

The recommendations for further work must include 
surveys with multiple sensor rakes oriented in the 
transverse direction and placed close to the wall in 
order to study the so called "streaky" structure in 
regular and manipulated boundary lavr'^s. Two-dimensional 
velocity reconstructions and pattei ' jognition would 
give an altogether new view of these f.ow structures and 
yield quantitative statistics cf their spacing with 
changing wall layer parameters. In particular we can 
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determine if they scale in accordance with the results 
for the wall events. 


The experiment should also be extended to higher 
Reynolds numbers and greater distances downstream of the 
manipulator. Preferably, the Reynolds numbers should 
fall in the range or exceed the proposed equivalent 
Reynolds numbers of the manipulated layer so that direct 
comparisons can be made. In this work, additional 
velocity^ surveys farther out in the boundary layer should 
be conducted in order to calculate the intermittency 
distribution and determine its Reynolds number 
dependence. In the present investigation, that data was 
not available. Pattern recognition may be quite useful 
in detecting the turbulent/non-turbulent interface at the 
edge of the boundary layer. It may also be very 
interesting to subject the manipulated boundary layer 
flows to strongly adverse pressure gradients in order to 
see if the large scales reemerge. 

In terms of the image processing, a great deal more 
analysis can be done on the present images. This 
includes two-dimensional pattern recognition of bursts 
from the long-view images. In future measurements, 
however, images more analagous to the Eulerian velocity 
acquisition of the hot-wire rake can be obtained through 
the illumination of the smoke along a line by a laser 
beam and scanning the image across the digital camera 
with a rotating mirror. This would give the needed time 
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aspect to the digital images and should yield some very 
interesting results after digital enhancement and 
frequency analysis. 
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INSTRUMENTATION DESIGN CONFIGURATION 


AS indicated in Chapter II, a cake of 11 hot-wires 
wao used to spatially sample the velocity field in the 
turbulent boundary layers. It was cost and space 
prohibitive to utilize commercially available constant 
temperature anemometers, to operate the hot-wires or 
amplifiers and low-pass filters, to pre-condition the 
analog signals before digitization and acquisition. We 
therefore designed and built constant temperature 
anemometers and analog signal conditioners which were low 
in cost and compact in size. The latter criterion was 
important because we wanted to place the anemometer 
circuitry relatively close to the sensors in order to 
minimize electronic noise and provide frequency 
compensation for the cables which connected the 


anemoneters. * 

The circuit diagram for the CTA was shown in Figure 
The etch pattern for that circuit is shown on the top 
of Figure 128. This was designed for single-sided, 
copper-clad printed circuit board. The photo negatives 
of this pattern were placed together, in line, in groups 
of three so that the parallel lines formed a continuous 
connection. These lines carried the *15 volt supply to 
the individual circuits. This arrangement was used as 
the pattern on a 6 x 22 cm photo-negative printed circuit 
board from which the copper was etched to leave the 
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conducting lines indicated by the black regions in this 
figure. Two of these boards, containing six anemometers, 
were stacked in an 8 x 8 x 28 cm component box which 
contained the required BNC connectors and switches, A 
total of 3 complete sets, corresponding to 18 channels, 
were built during this experiment. 

Included with each group of six anemometers was a 
square wave generator based on a 7405, Hex— inverter . The 
etch pattern for that simple circuit is shown in the 
middle of Figure 128. This was sucessively connected to 
each anemometer bridge through a rotary switch, which 
also connected the anemometer output to a common BNC for 
monitoring the response during feedback adjustments to 
the anemometer bridge. The square wave generator was 
turned off after adjustments were completed. 

The pre-conditioning of the analog outputs from the 
anemometers was done by the bias/gain/filter circuit 
shown in Figure 7. The etch pattern for each channel of 
that circuit is shown on the bottom of Figure 128. In 
this case, six channels were compactly fit into a 15x8x25 
cm. component box. Each channel utilized external 
connectors to allow for different feedback resistors to 
set the gain of each amplifier stage. Also, each channel 
incorporated a single rotary switch which simultaneously 
connected two sets of matched precision resistors to 
change the cut-off frequency of the low-pass filter. The 
phase and frequency characteristics of all of the filters 
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designed in this experiment were adjusted to give the 
identical response to that shown in Figure 7. External 
connectors were also available on the component box to 
allow for the use of different resistor combinations for 
different frequency cut-offs. 

In all three circuits, special attention was placed 
on keeping electronic noise to a minimum. In this 
regard, all power inputs t« the generator amplifiers were 
DC coupled with .01 ;.f capacitors. Also, all external 
leads were properly shielded so as not to be sources of 
power line noise. 
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SUMMARY OF RESULTS FROM A PARALLEL EXPERIMENT 

A parallel experiment on the control of the large- 
scale intermittency in a turbulent boundary layer was 
performed by H. Nagib, in collaboration with R. 
v:eatphall, during the period from December 1980 to 
Ferruary 1961, while on sabbatical leave at the HTT?* 
Laroratory of the Mechanical Engineering Department of 
Stanford University. In that work, a turbulent boundary 
lever was developed on one wall of a 2.4 m long 
wird-tunnel test section which had been adjusted to give 
a tear zero streamwise pressure gradient. The 
crtss-sectional dimensions of the tunnel were 0.2 on by 
0.5 m with the measurements being performed on the longer 
dimensioned side. The point of transition to turbulent 
flow was fixed at the entrance by sandpaper in a manner 
similar to that used in this inves^gation. The 
torbulent manipulator consisted of two 0.25 mm thick 
plates which were suspended in the tunnel in an 
arrangement of the M-4 type,i.e., very similar 
configuration as that of the final results of the present 
investigation. The streamwise chord of the plates was 
aoproximately equal to the boundary layer thickness, 5, 
as their distance from the wall was approximately 0.5:. 
Toe streamwise spacing between the plates corresponded to 
t'- as in the present investigation. Therefore, the nain 
cofference in the design of the manipulator in the two 
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investigations is the elevation of the two tandem plates 
from the wall. Here we used 0.86 while Nagib used 0.56. 

The manipulator was placed in the boundary layer at 
the streamwise location x = 1.3 m where the Reynolds 
number based on the momentum thickness was approximately 
1700. The free-stream velocity was 10.6 m/s and the 
turbulence intensity at that speed was approximately 0.5 
percent. The free-stream velocity fluctuations were 
mainly a result of low frequency contributions. 

The measurements consisted of mean velocity and 
turbulence intensity profiles which were analyzed to 
provide the mean properties of the flow including the 
integral thicknesses, shape factors, and friction 
coefficients. These were obtained for 9 stations along 
the streamwise centerline as well as 2 off center 
stations, in addition to the mean prof iles,- pulsed-wi re 
probe measurements were made near the wall at each of 
these locations to evaluate the mean wall-friction 
coefficients, their rms and probability density 
functions. The hot-wire and pulsed-wire data were 
acquired digitally and stored on floppy disks. The data 
presented here are based only on the hot-wire data and 
were processed at IIT using the same routines developed 
for this investication . Various problems were encountered 
with the acquisition and interpretation of the 
pulsed-wire data and, therefore, they are not summarized 
here. 
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Smoke-wire visualization was also employed to give 
visual evidence of the performance of the turbulence 
manipulator and to aid in optimizing the device. The 
results to follow are intended to summarize the work of 
tnat study. Additional details may be published later on 
tnis parallel effort. In particular, the evidence from 
tr.e smoke wire visualization in the form of photographs 
is very similar to that presented here. 

Figure 129 documents the mean velocity and 
turbulence intensity profiles just upstream of the 
streamwise location of the device. The dashed and solid 
lines in the figures correspond to the distributions with 
and without the m.anipulator in place, respectively. In 
all these figures, the mean and rms velocities are 
ncn-d imens ional i z ed with the local free— stream velocity. 
This independence of the initial condition^ from the 
presence cf the manipulator is amply verified by this 
figure. 

The distributions just downstream of the device 
(x/h=2.3) are shown in Figure 130. These clearly show 
the wake cf the manipulator plates. Also evident is the 
lever u'/Go values in the manipulated boundary layers, 
which was also observed by Corke et al. (1980) for the 
fcur-plate manipulator. The velocity defect in the 
manipulator wakes was larger in this work by Nagib 
because of the somewhat thicker plates and larger chord 
lengths. As will be seen, this will result in a somewhat 
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smaller net drag reduction in this investigation for a 
distance from the device comparable to that for the 
present work. This added drag is also evident in the 
larger increase in Re^^for the manipulated boundary 
layer at the same free-stream speed, i.e., larger 
increase in momentum thickness due to the device. 

Farther downstream of the device in Figures 131 and 
132, the wake of the device be dissipated by the 
turbulence in the boundary layer. By 21 device heights 
downstream of the manipulator in Figure 131, the nomentum 
thickness in the manipulated layer is still somewhat 
larger than in the regular case but it clearly is not 
growing as quickly. By the last station, 90 device 
heights downstream, the mean velocity profile 
demonstrates the same characteristic differences as snown 
in Figures 25 and 26 of the present work. The lower 
turbulence intensity in the outer part of the boundary 
layer indicates a reduction in the velocity fluctuations 
associated with the large scale excursions of turbulent 
fluid away from the wall. By this station, the momentum 
thickness is less than that exhibited at the same 
x-location in the regular layer; i.e., without the 
device. 

The Reynolds number distribution of the local 
wall-friction coefficients is shown in Figure 133 for the 
regular and manipulated boundary layers. The results of 
Nagib are plotted along with the C,, distributions 
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obtained from the present investi^^ation. The values for 
for Nagib’s data were obtained from the same momentum 
balance software written for the present Investigation 
and summarized in Chapter IV- The values for the regular 
boundary layer fall in close agreement with average 
values obtained for the present results. The error bars 
for those data points correspond to the range of values 
obtained from the momentum balance, the Ci^au^er log-fit 
and the Ludwieg-Tillmann correlation. 

The results in the manipulated boundary layer 
demonstrate a comparable reduction in the local friction 
coefficient to the results of the present work. The 
error bars in these data points are a result of 
uncertainty in measuring dP/dx along the test section. 

In one limit, the distribution was calculated 

utilizing the momentum balance with the streamwise 

* 

pressure difference of zero set to an accuracy of 0.03 mn 
of water by Ilagib over the entir® length of test section. 
The other linit of the error bars was obtained for the 
du/dx distribution which was returned from the analysis 
and is based on the hot-wire readings in the free stream. 
Both values show approximately 30 percent reduction in 
the local friction coefficient. If we i 'dude the added 
drag of the device, based on its momentum defect, the 
results of Nagib indicate a net drag reduction of 
apprcximatcly 15 percent with a fetch of less than 60 
boundary layer thicknesses. 
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The results of this experiment by Nagib demonstrate, 

in a separate facility, the effectiveness of the series 

type manipulator in reducing the drag in turbulent 

boundary layers. They further substantiate the findings 

of the present work and document the effects of altering 

the intermittent scales in a turbulent boundary layer. 

Finally, they clearly identify one of the mechanisms of 

this type. The reduction in turbulence level beneath and 

• % 

in the lower part of the plate wake, depicted in Figure 
130 , is very important to the operation of these devices. 
Here we can see that the transfer of energy to the mean 
flow in that part of the wake is partly balanced by a 
transfer back into the turbulence from,the mean flow 
above the plate elevation. The change in the sign of 
across the perturbation introduced to the equilibrium 
boundary layer flow in the form of the wake is a very 
important ingredient to this two-way transfe?. While the 
reduction of turbulence near the wall is the important 
ingredient, this two-way transfer may make its 
implementation very efficient. 
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Figure 6. Circuit Diagram of DC Power Supply Used for 12-Channel C.T.A, 
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Fiqure 7. Circuit Diagram of Variable DC Bias Amplifiers 
P^eproceas Analog Signals Before Digitization 


and Low-Pass Filters Used to 


oK!g;.\al page fs 

OF POOR QUALITY 




220 



I'iyuro 0. f'roquoncy Rosponae Curvoa for Socond-Ordcr Duttervorth Low-Pass Filters Shown 
in Fiqure 7 


ORIGINAL PAGE 13 
OF POOR QUALITY 



221 


I 18 
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Figure 9. Schematic of PDP-11 Based Fast Acquisition System 
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DAPS PROCESSING AND DEVELOPMENT SYSTEM (PDS) 


Figure 10. Schematic o 


f PDP-11 Based Software Development and Data Processing System 
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DAPS SPECIAL PURPOSE MODULE (SPM) 


Figure 11. Schematic of Special Purpose Module Containing Microprocessor-Controlled 
Array and Display Processors 
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Figure 13 . Schematic of Instrumentation Used in Mean Profile Measurements 
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Fxgure 23. Plate Configurations for Boundary Laver 
Manipulators ^ 
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Manipulators 
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Figure 27. Comparison of Mean Velocity and Turbulence Intensity Pr . files for Boundary 
Layer Alone and with Optimized Manipulators 
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Figure 28. Visualized Flow Without and With the Addition 
of Two-Plate Manipulator, Showing Streamwise Extent (7.05) 
of Large-Scale Turbulence Suppression 



2*^0 





STATION, X(m) 

§ 

0 

1, 2.25 

V 

□ 

3, 2.46 

6 

< 

5, 3.37 


V 

7. 3.98 

i 

o 

9, 4.59 

Q 

Ra » 

2005-3319 








O 0.3 0.4 0.5 0.9 1.0 

Figure 29. Non-dimensional Mean Velocity Profiles at Variou.'? Streamwiso Stations in 
Regular Boundary Layer for Low Reynolds Number Range 


ORIGINAL PAGE »3 
OF POOR QUALITY 










STATION, 

x(m.) 

V 

0 1 , 

2.25 

8 

□ 3 , 

2.46 


< 5 , 

3.37 

% 

^ 7. 

3.98 

\ 

O 9 . 

4.59 


Rec = 2005-3819 
®2 


/ FROM KLEBANOFF 




ro riles at 
: Range 


ORlGiriAL rAGi 
OF POOR QOAUn^ 






.6 2.0 2.4 2.8 


ofiloa nt VarloviP fltrcnmwiao fitatlona 
ir Ranqe 




x(m) 


l-’iyuru JJ. 


Utrortinwibo Udvelupniont of Nun-dlmuiialon«l |"i oo-UL turtni I'l'oatMuo Uirroiomjo 




Fiqiire 3'’i. Compnrison of WaU-Friction Coofficiont Values Obtaincrl hy Throe Methocls for 
I Ml MimmmImi y t<ayMi 



248 


II I III Hill 


Hill I II II I lll■llllll 






ORIGINAL 
OF POOR quality 


tfZ 



r 


ORIGINAL PAGE 13 
OF POOR QUALITY 



ORIGINAL PAGE IS 
OF POOR QUALITY 



sbuBH -laqumN spTOuAan «oq Jtoj jaAcq Ajtepunoe 

snoT^BA uoxqnqiaqsTa AqxooiaA BuxqBnqonu pazxxBuoxsuauixp-uoN 

♦ X 


suoxqBqs asxMUiBajqS 

^T“n-IT®M *8C a-inBxj 



o 

m 

CM 


uVU 



*A 


D009 


0001 009 


h 

90M-S662- '‘•M 


69> 

869 


^At>a 

o 


910 ‘ C I 

9Z'Z ' I ^ 

(ui)¥ • Noime 




252 



Piaure 40. wake Function Distribution at Various Streamwlse stations in fiegular Boundary 
Layer for Ixjw Reynolds Number Range 





Figurf! 41. Wake Function Distribution at Varioi 
Layer for High Reynolds Number Range 
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Figure 42. One-dimensional S^ctra of the Streamwise 
Velocity Component at Different Self-Similar Heights 
at Station 2 in the Regular Boundary Layer for Low 
Reynolds Number Range 
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Figure 43. One-dimensional Spectra of the Strecunwise 
Velocity Component at Different Self-Similar Heights 
at Station 5 in the Regular Boundary Layer for Low 
Reynolds Number Range 
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Figure 45. Cne-dimensional Spectra of the Streamwise 
Velocity Component at Different Self-Similar Heights 
at Station 2 in the Regular Boundary Layer for High 
Reynolds Number Range 
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Fiqure 46. One-dimensional Spe9tra of the Stre^ise 
Velocity Component at Different Self-Similar Heights 
at Station 5 in the Regular Boundary Layer for High 
Reynolds Nuirlser Range 
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Figure 47. One-dimensional Spectra of the Streamwise 
Velocity Component at Different Self-Similar Heights 
at Station 9 in the Regular Boundary Layer for High 
Reynolds Number Range 
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Figure 50. Reynolds Number Dependence of Non-dimensional Streamwise Integral Length Scale 
nt Uifferwnt J!«l C-JJ Iml Inr Ualghtn In t.'H iu?gular nmiridary Layer 




Figure 51. Two Realizations of Smoke-wire Visualization 
Photographs of Flowfield in Regular Boundary Layer (Top) 
and Manipulated Boundary Layer (Bottom) Representative 
of Final Data Runs 
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Fiquro 52. Non-dimensional Mean Velocity Profiles at Various Streamwise Stations in 
Manipulated Boundary Layer for Low Reytiolds Number Range 
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Fiqure 56. Comparison Hetween Streamwise Distribution of Integral Thicknefises in 
Regular and nampulated Boundary Layers for Low Reynolds Number Range 
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Figure 57. Comparison Between Streamwisa Distribution of Integral Thicknesses in 
Regular and Manipulated Boundary Layers for High Reynolds Number Range 
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Figure 58 . Reynolds Number Dependence of Shape Factor for Regular and Manipulated 
Boundary I^ayers 
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Figure <>0. Comparison of Wall-Friction Coefficient Values Obtained by Three Methods 
in Manipulated Boundary Layer and Average Distribution in Regular Boundary Layer 
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Figure 64, Wall-Unit Non-dimensionalized Fluctuating Velocity Distribution at Various 
Streamwise Stations in Manipulated Doundary Layer for High Reynolds Number Range 
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Figure 65. Wake Function Distribution ai. Various Streamwise Stations in Manipulated 
Boundary Layer for Low Reynolds Number Range 
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Figure 66. Wake Function Distribution at Various Streamwise Stations in Manipulated 
Boundary Layer for High Reynolds Number Range 
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Figure 69. Reynolds Number Dependence of Non- 
Correlation to Wall-Fiiction Coefficient 
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Ficure 70. One-diir.ensional Spectra of Streamwise Velocity 
Component at Various Stations Behind Manipulator Plate 
for Low Reynolds Nuir±>er Range 



EU/VflR(U) 





EU/VRRIU) 


ORIGINAL PAGE IS 
OF POOR QUALITY 


HRNIPULflTED BL. 
STRTION 1 
LO 


K (l/M) 

Fiaure 72. One-dimensional Spectra of the Streamwise 
velocity Component at Different Self-Similar Heights 
at Station 1 in the Manipulated Boundary Layer for Low 
Reynolds Number Range 
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Figure 73. One-dimensional Spectra of the Streamwise 
Velocity Component at Different Self-Similar Heights 
at Station 4 in the Manipulated Botandary Layer for Low 
Reynolds Number Range 
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Figure 75. One-dimensional Spectra of the Streamwise 
Velocity Component at Different Self-Similar Heights 
at Station 1 in the Manipulated Boundary Layer for High 
Reynolds Number Range 
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Figure 77. One-dimensional Spectra of the Streamwise 
Velocity Component at Different Self-Similar Heights 
at Station 8 in the Manipulated Boundary Layer for High 
Reynolds Number Range 
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Figure 78. One-dimensional Spectra of Wail-Shear Stress 
for Regular and Manipulated Bo’undary Layers at Both Free- 
Stream Speeds 
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Fiqure 81 . Vertical Distribution of Skewness of Streainwise 
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Figure 82. Vertical Distribution of Kurtosis of Streanwise 
Velocity in Regular and Manipulated Boundary Layers 
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Figure 83. Vertical Distribution of Non-dir.ensional 

Streanwise Integral Length Scale at Various Downstream 
Stations in Manipulated Boundary Layer for Low Peynolds 
Number Range 
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Fiqure 34. Vertical Distribution of Non-dimensional 

Stroamwise Integral Length Sqale at Various Downstream 
Stations in Manipulated Boundary Layer for High Reynolds 
Number Range 
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Figure 37. Kormalized Time Series Outputs From Kot-wire 
^ke and Surface Mounted Shear Wire for Regular Boundary 
Layer at High Reynolds Number 
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Figure 90, Normalized Burst-Conditioned Ensemble-Average 
Time Series Obtained from VITA Detection in Regular 
Boundary Layer at Low Reynolds Number 
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Figure 91. Normalized Probability Distributions of Burst 
Periods from VITA Detection with Different Thresholds in 
Regular Boundary Layer at Low Reynolds Ikumber 
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Ficure 54. Normalized Velocity Time Series From Serscr at 
y+ = 15 and Corresponding Match Filter Output i^ 
Manipulated Boundary’' Layei for Low Reynolds Numh-er 


306 




OF POOR QUALITY 




.54QfTO (Q(y) -Q(y))/Qrms(y^ 


Dar.t IS 



TU^/5 


Fiaure 96. Normalized Burst-Conditioned 
Time Series Obtained from Match Filter 
Regular Boundary Layer at Low Reynosos 
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99. Nornijlizod Burst-Conditioned Ensemble Average 
Time Series Obtained from Match Filter Detection in 
Manipulated Boundary Layer at High Reynolds Number 
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Ficure 104. Normalized Probability Distributions of Burst 
Periods from Match Filter Detection in Regular and 
Manipulated Boundary Layers at High Reynolds Number 
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Figure 109. Ensemble Averaged Two-Dimensional Boun^ 

VolocUy for Time Instant of ni^ital Image Acquisition in Manipuiacea ao y 
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Figure 112. Burst-Conditioned Ensemble Average of Normalized Two-Dimensional Velocity ^ 
Reconstructions Obtained From Match Filter Detection in Manipulated Boundary Layer 
at High Reynolds Number 
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Fiqure 113. Ensemble Average of Normalized Two-Dimensional Velocity Reconstructions 
for Events with Periods Greater than Three Standard Dev.\.ations from Distribution 
Mean at Low Reynolds Number 
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Picpj-r-e 115. Normalized Two-Dimensional Spectrum of 
Reconstructed Streamwise Velocity Distributions fo 
Regular Boundary Layer at Low Rei-nolds Number 



igure 116. Normalized Two-Dimensional Spectrum of 
Reconstructed Streamwise Velocity Distributions for 
Regular Boundary Layer at High Reynolds Number 
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Figure 118. Normalized Two-Dimensional Spectrum of 
Reconstructed Streamwise Velocity Distributions for 
Manipulated Boundary Layer at High Reynold.*; Number 
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'iqure 119 . Digitized Realization of Visualised Flow 
Conditioned on Wall Event (Top) and Corresponding 
Twc-Dinensional Match-Filter Enhancement in --^egu-ar 
Boundar-/ Laver 


ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 




Figure 120. Digitized Realization of Visualized Flow 
Conditioned on Wall Event (Top) and Corresponding 
Two-Dimensional Match-Filter Enhancement in 
Manipulated Boundary Layer 
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Figure 121. EnseirJjie Average of Enhanced Images 

Conditioned on Match Detection in Regular (Top) and 
Manipulated (Bottom) Boundary Layers 


334 


ii SiiiiiilBliiiiliiiiii 


i ii III II II! nil mil II 


i: miisiBi 


SI I! !l 


lll!l!l III I I lilt H9I 



Figure 122. Digitized Realization of Visualized Flow conditioned on Wall Event (Top) 
and Corresponding Two-Dimensional Match-Filter Enhancement (Bottom) Viewing 106 
Regular Boundary Layer 
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Figure 123. Digitized Realization of Visualized Flow Conditioned on Wall Event (Top) 
and Corresponding Two-Dimensional Match-Filter Enhancement (Bottom) Viewing 106 
in Manipulated Boundary Layer 
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Figure 124. Roynolda Number Dependence of Wall-Friction Coofficionta in Regular and 
Manipulated Doundary Layers 
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Figure 125. Reyr.olcs y"® 
Strouhal Frequency in .-^eg 
Layers 


Decendence of Wall-Layer 
ar'and Manipulated Boundary 
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Figure 129. Mean Velocity and Turbulence Intensity- 
Prof lies in Regular and Manipulated Soundary Layers 
Upstrean of Device Position 
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Figure 132. Mean Velocity apd Turbulence Intensity 
Profiles in Regular and Manipulated Boundary Layers 
at Location Approximately 90 Device Heights Downstream: 
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Figure 131. Mean Velocity and Turbulence Intensity 
Profiles in Regular and Manipulated Boundary Layers 
at Location Approximately 21 Device Heights Downstream 
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Figure 130. Mean Velocity and Turbulence Intensity 
Profiles in Regular and Manipulated Boundary Layers 
at Location Approximately 2 Device Heights Dcwnstrear: 
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